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Vorwort des Herausgebers 


Die Fahrzeugtechnik ist kontinuierlich Veränderungen unterworfen. 
Klimawandel, die Verknappung einiger für Fahrzeugbau und -betrieb 
benötigter Rohstoffe, globaler Wettbewerb, gesellschaftlicher Wandel 
und das rapide Wachstum großer Städte erfordern neue Mobilitätslö- 
sungen, die vielfach eine Neudefinition des Fahrzeugs erforderlich ma- 
chen. Die Forderungen nach Steigerung der Energieeffizienz, Emissi- 
onsreduktion, erhöhter Fahr- und Arbeitssicherheit, Benutzerfreund- 
lichkeit und angemessenen Kosten sowie die Möglichkeiten der 
Digitalisierung und Vernetzung finden ihre Antworten nicht in der sin- 
gulären Verbesserung einzelner technischer Elemente, sondern benöti- 
gen Systemverständnis und eine domänenübergreifende Optimierung 
der Lösungen. 


Hierzu will die Karlsruher Schriftenreihe für Fahrzeugsystemtechnik 
einen Beitrag leisten. Für die Fahrzeuggattungen Pkw, Nfz, Mobile Ar- 
beitsmaschinen und Bahnfahrzeuge werden Forschungsarbeiten vorge- 
stellt, die Fahrzeugsystemtechnik auf vier Ebenen beleuchten: das Fahr- 
zeug als komplexes, digitalisiertes mechatronisches System, die 
Mensch-Fahrzeug-Interaktion, das Fahrzeug in Verkehr und Infrastruk- 
tur sowie das Fahrzeug in Gesellschaft und Umwelt. 


E-Motoren sind als Antriebsmaschinen in Personenkraftwagen gut ge- 
eignet. Wie auch Verbrennungskraftmaschinen tragen sie zur Ge- 
räuschsituation im Fahrzeug bei und beeinflussen so Komfort und Qua- 
litätseindruck entscheidend mit. Es ist daher wichtig, über ein entspre- 
chendes Instrumentarium zu verfügen, um die Geräuschentwicklung 
von E-Maschinen passend zum Charakter des Fahrzeugs gestalten zu 
können. 


Vorwort des Herausgebers 


Hier setzt die vorliegende Arbeit an. Sie befasst sich nicht nur mit den 
vielfältigen Ursachen von Wechselkräften, die durch eine E-Maschine 
erzeugt werden und zu Geräuschen führen, sondern nutzt die Möglich- 
keit, durch die Ansteuerung der Maschine diesen Schwingungen entge- 
genzuwirken. Dazu erarbeitet Herr Stretz schnell rechnende Methoden 
für die Modellierung der E-Maschine und der Schwingungsweiterlei- 
tung über Strukturbauteile sowie zur Bestimmung von Betrag und 
Phase von Überlagerungsströmen zur Schwingungskompensation. Die 
Methodik erlaubt auch, Geräuschursachen zu ermitteln. Ihre hohe 
Wirksamkeit wird in breit angelegten Versuchsreihen belegt. 


Frank Gauterin 


Karlsruhe, im Juli 2022 
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Kurzfassung 


Der Einsatz elektrischer Radnabenantriebe als Traktionsantrieb bringt 
viele Vorteile mit sich. Einer der Wichtigsten ist der Gewinn zusätzli- 
chen Bauraums innerhalb der Fahrzeugkarosserie, da dieser durch die 
Integration des elektrischen Motors und der Leistungselektronik im Rad 
frei wird. Zusätzlich ergibt sich die Möglichkeit von sehr großen oder 
auch radindividuellen Winkeleinschlägen durch den Wegfall der An- 
triebswellen. Zentrale Traktionsantriebe werden aus vibroakustischen 
Gründen über Entkopplungselemente an der Fahrzeugkarosserie befes- 
tigt, um die Übertragung von im Motor entstehenden Vibrationen in den 
Fahrzeuginnenraum zu minimieren. Radnabenantriebe hingegen wer- 
den direkt am Fahrwerk beziehungsweise der Radnabe montiert. Der 
Einsatz von erweiterten Entkopplungsmechanismen ist weder zwischen 
Antrieb und Fahrwerk noch zwischen Fahrwerk und Karosserie mög- 
lich, ohne das Fahrverhalten zu beeinflussen. Denn neben den Antriebs- 
kräften müssen über diesen Weg auch alle anderen fahrdynamischen 
Kräfte aufgenommen werden. Vor allem höherfrequente Anregungen, 
welche durch die Magnetkräfte des Motors entstehen, werden gut über 
das Fahrwerk in den Fahrzeuginnenraum übertragen, da bisherige Fahr- 
konzepte auf eine derartige Anregungsquelle nicht ausgelegt sind. Dies 
führt schließlich zu störenden Geräuschen und Vibrationen in der Fahr- 
zeugkabine. 


Das im Projekt „E-Wheel-Drive“ verwendete Fahrwerkskonzept einer 
Verbundlenkerachse mit Radnabenantrieb besitzt diese vibroakustisch 
negative Eigenschaften und wird aus diesem Grund zu einer näheren 
Untersuchung herangezogen. 


Im Rahmen dieser Arbeit werden die Entstehungsmechanismen dieser 
störenden Geräusche und Vibrationen untersucht und in einem multi- 
physikalischen Modell nachgebildet. Aus diesen Erkenntnissen werden 
Gegenmaßnahmen abgeleitet und bewertet. Die Gegenmaßnahme mit 
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dem größten Potential, ein vibroakustisch optimierter Steuerungsan- 
satz, wird weiterentwickelt und auf das zu untersuchende System so- 
wohl simulativ als auch experimentell angewandt. Des Weiteren wer- 
den Grenzen und das Zusammenwirken des optimierten Steuerungsan- 
satzes mit der Motorregelung näher beleuchtet. Aus den gewonnenen 
Ergebnissen wird eine Betriebsstrategie für die spätere Nutzung abge- 
leitet. 


Der in der Arbeit dargestellte vibroakustisch optimierte Steuerungsan- 
satz ist ebenso auf andere elektrische Maschinen übertragbar. Mit der 
entwickelten Methodik können die dafür notwendigen Parameter in den 
verschiedenen Betriebspunkten durch nur wenige Messungen schnell 
bestimmt werden. 
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Abstract 


The use of electric wheel hub drives as a traction motor has many ad- 
vantages. One of the most important is the gain of additional space 
within the vehicle body, which becomes available through integration 
ofthe electric motor and the power electronics in the wheel. In addition, 
very large or wheel-individual wheel steering angles can be realized, 
since no drive shafts are needed anymore. However, in contrast to 
powertrain concepts with a central engine, all forces and moments of 
the wheel hub drive to the body act at the connection point to the chas- 
sis. Thus, the resulting forces from the electromagnetic fields of the 
electric machine are routed through the chassis into the passenger cabin, 
where they can lead to disturbing noise and vibrations. 


Within the scope of this work the root causes of these disturbing noises 
and vibrations are investigated and simulated using a Multiphysics 
model. Based on the insights obtained through these simulations, coun- 
termeasures are derived and evaluated. The countermeasure with the 
greatest potential, a vibroacoustically optimized control strategy, is de- 
veloped and applied to the wheel hub drive system wheel hub drive in 
simulation and experiment. The limitations and the interaction of the 
optimized control strategy with the engine control are analyzed as well. 
From these results, an operating strategy for later use is derived. 


The vibroacoustically optimized control strategy presented in this thesis 
can also be applied to other electric motors. With the developed meth- 
odology, the required parameters can be determined quickly at various 
operating points by only a few measurements. 
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1 Einleitung 


1.1 Motivation 


Mit dem Einsatz von Elektromotoren als Antriebsmaschine in Fahrzeu- 
gen sind neue, innovative Fahrzeugkonzepte denkbar. Allerdings erge- 
ben sich hieraus auch neue Herausforderungen wie der geänderte vibro- 
akustische Charakter eines Elektrofahrzeuges gegenüber eines mit Ver- 
brennungsmotor angetriebenen Fahrzeuges. 


Eine Ausprägung dieser neuen Fahrzeugkonzepte bilden Fahrzeuge mit 
elektrischen Radnabentrieb. Hierbei sind alle für den Antriebsstrang 
notwendigen Komponenten wie Leistungselektronik, Motor und 
Bremse in der Felge integriert. Dieses Konzept bietet neben der Mög- 
lichkeit von erhöhten und unabhängigen Radwinkeleinschlag [Tin09] 
und dem direkten Kraftfluss durch den sehr kurzen Übertragungsweg 
zur Straße vor allem die Möglichkeit der Vergrößerung des Nutzraums 
innerhalb der Karosserie [Dep16]. 


Im Gegensatz zu Fahrzeugen mit elektrischen Antriebsstrang sind die 
bei konventionellen Antrieben entstehenden Wechselkräfte, welche zu 
hör- und spürbaren Vibrationen im Innenraum führen, eher niederfre- 
quent. Der Mensch hat sich an diese niederfrequenten Vibrationen und 
Geräusche der Verbrennungskraftmaschine, die in der Fahrzeugkabine 
wahrgenommen werden, gewöhnt. Oft werden diese Effekte vom Fahr- 
zeughersteller sogar gezielt genutzt, um ein Fahrzeug beispielsweise 
besonders sportlich wirken zu lassen. Die ungewohnten höherfrequen- 
ten Geräusche und Vibrationen, die aus der elektrischen Antriebsma- 
schine eines E-Fahrzeuges stammen, wirken hingegen auf den Men- 
schen eher störend [Suk14], [Mar12]. Besonders unangenehm für den 
Fahrer sind diese Effekte bei niedrigen Geschwindigkeiten, da hier 
noch die Maskierungseffekte der vorbeiströmenden Luft und die Roll- 
geräusche des Reifens sehr gering sind. Somit treten die Geräusche und 
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Vibrationen der elektrischen Maschine dort besonders stark in den Vor- 
dergrund [Alb14]. Bei höheren Geschwindigkeiten dominieren dann 
Geräuschanteile, die von Wind und Reifen verursacht werden und die 
bei niedrigen Geschwindigkeiten als störend wahrgenommenen Geräu- 
sche der elektrischen Maschine rücken in den Hintergrund. 


Da der Komfort eine der wichtigsten Eigenschaften eines Fahrzeuges 
ist und dieser stark durch die in der Fahrgastzelle wahrgenommen Ge- 
räusche und Vibrationen geprägt wird, trägt dieser zur Kaufentschei- 
dung eines Fahrzeuges wesentlich bei, wie Abbildung 1.1 zeigt. 
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Abbildung 1.1: Komfort als einer der Hauptgründe für die 
Kaufentscheidung [Ara15] 


Somit ist es vor allem bei Elektrofahrzeugen von großer Bedeutung, 
störende Geräusche und Vibrationen zu minimieren, um zu einem mög- 
lichst hohen Komforteindruck in den Fahrzeugen beizutragen. 
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Ein effektiver Ansatz, die Übertragung der durch die elektrische Ma- 
schine entstehenden Wechselkräfte auf die Karosserie und in den In- 
nenraum zu reduzieren, ist eine Entkopplung durch geeignete Lagerung 
der Antriebsmaschine [van16], wie es bereits auch bei Fahrzeugen mit 
Verbrennungsmotor der Fall ist. Aufgrund der sehr ähnlichen Topolo- 
gie mit zentraler Antriebseinheit kann diese effektive Möglichkeit bei 
Fahrzeugen mit zentral oder radnah angeordneten Elektromotoren an- 
gewandt werden. Im Falle des Radnabenantriebs ist dies jedoch auf- 
grund der notwendigen steifen Anbindung des elektrischen Motors an 
das Fahrwerk nicht möglich, da über die Anbindung auch die komplet- 
ten Antriebs- und Seitenführungskräfte des Fahrzeugs übertragen wer- 
den müssen. 


Somit ist es besonders für dieses neuartige Fahrzeugkonzept notwen- 
dig, die durch die elektrische Maschine erzeugten störenden Effekte ge- 
nauer zu untersuchen und Potentiale für eine Optimierung der vibroa- 
kustischen Eigenschaften aufzuzeigen. Besonderes Interesse gilt dabei 
den kritischen niederfrequenten (<200Hz) Wechselkräften, welche bei 
kleinen Geschwindigkeiten beziehungsweise Raddrehzahlen entstehen. 
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1.2 _Ausgangssituation 


Vibroakustische Auffälligkeiten in der Fahrgastzelle konnten bereits in 
der Entwicklungsphase des elektrischen Radnabenantriebs festgestellt 
werden. Für die hier dargestellte Forschungsarbeit wird die 2. Genera- 
tion des von der Schaeffler Technologies AG & Co. KG entwickelten 
Radnabenantriebs genauer untersucht (E-Wheel-Drive ,,beta“). Abbil- 
dung 1.2 zeigt den Schalldruckpegel bei einer beschleunigten Fahrt mit 
diesem neuartigen Antriebssystem, welches am Ford Fiesta als Ver- 
suchsträger zum Einsatz kam. Diese Abbildung zeigt bereits die vom 
elektrischen Motor ausgehenden hörbaren Effekte. Deutlich zu erken- 
nen sind die niederfrequenten Ordnungen (40., 47., 120., 240.), welche 
von der elektrischen Maschine ausgehen und zu störenden Effekten im 
Fahrzeuginnenraum führen. 
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Abbildung 1.2: Schalldruckpegel im Fahrzeuginnenraum am Akustik- 
rollenprüfstand (A-bewertet) 
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1.3 Ziel und Aufbau der Arbeit 


Aus diesem Grund werden in dieser Arbeit Informationen über das Zu- 
sammenwirken der elektromagnetischen und mechanischen Effekte 
strukturiert aufbereitet, die helfen den gesamten Entstehungsprozess 
der Vibrationen, welche vom Elektromotor ausgehen, zu verstehen und 
Gegenmaßnahmen daraus abzuleiten. 


Ziel der Arbeit ist die Analyse der aus dem Motor stammenden Kräfte 
und die Entwicklung einer Methodik zur Optimierung der vibroakusti- 
schen Eigenschaften von Fahrzeugen mit elektrischen Radnabenan- 
trieb. Besonderer Wert wird hier auf die Möglichkeit eines optimierten 
Steuerungsansatzes gelegt, welcher großes Potential für eine vibroakus- 
tische Optimierung bietet. Dieser soll helfen, den vibroakustischen 
Nachteil gegenüber Fahrzeugen mit zentralem Antrieb zu minimieren. 


Die Wirksamkeit dieses Ansatzes wird in dieser Arbeit simulativ sowie 
experimentell nachgewiesen. 


Um diese Maßnahmen zur Reduzierung der störenden vibroakustischen 
Effekte auszuwählen und diese gezielt anwenden zu können, ist ein 
grundlegendes Verständnis der Quellen, der für die Vibrationen verant- 
wortlichen Wechselkräfte sowie der Übertragungspfade notwendig. 


Der Aufbau dieser Arbeit orientiert sich an der in Abbildung 1.3 ge- 
zeigten Vorgehensweise. 


1 Einleitung 


1. Einleitung 
2. Grundlagen 


3.1 Quellen 3.2 Übertragung 


> Elektromagnetik > Körperschall 


° 4.1 Aufbau 


4. Gesamtsimulation e 4.2 Simulationsergebnisse 
e 4.3 Validierung des Modells 


5. Identifikation von Gegenmaßnahmen 


m - 


o Akustisch optimierter Steuerungsansatz 


6.1 Anwendung der Methodik in der Simulation 
e 6.2 Anwendung der Methodik am Prüfstand 
e 6.3 Betriebsstrategie 
e 6.4 Vergleich verschiedener Methoden 
e 6.5 Messung am Gesamtfahrzeug 


7. Zusammenfassung / Ausblick 
A. Anhang 

- Prüfstandsaufbau 

- zusätzliche Messergebnisse 


Abbildung 1.3: Aufbau der Arbeit 


Zu Beginn der Arbeit (Kapitel 2) wird auf die notwendigen Grundlagen, 
die dem Verständnis der weiteren Ausführungen dienen, eingegangen. 
In Kapitel 3 werden Quellen der störenden Wechselkräfte und deren 
Übertragungswege aufgezeigt und beschrieben. So soll eine solide Ba- 
sis für die in der Arbeit betrachteten Effekte geschaffen werden. Aus- 
gehend vom Stand der Technik werden wesentliche Neuerungen, die 
zur Erklärung der Entstehung der störenden Effekte beitragen, genauer 
erläutert. Aufbauend auf diesen Beschreibungen und Modellen kann ein 
Gesamtmodell (Kapitel 4) erstellt werden, mit dessen Hilfe die wesent- 
lichen, zur vibroakustischen Störung beitragenden Effekte nachgebildet 
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1.3 Ziel und Aufbau der Arbeit 


werden können. Mit Hilfe der zuvor erarbeitenden Informationen wer- 
den im nächsten Kapitel Gegenmaßnahmen (Kapitel 5) identifiziert und 
bewertet. Die Gegenmaßnahme mit dem größten Potential, ein akus- 
tisch optimierter Steuerungsansatz, wird weiterverfolgt, da diese Maß- 
nahme auch auf andere elektrische Antriebe übertragbar ist. Die Wei- 
terentwicklung und Anwendung dieser Optimierung wird in Kapitel 6 
systematisch erklärt. Hierfür wurde eine ganzheitliche Methode entwi- 
ckelt, basierend auf aufgeprägte Schwingungen in den Strömen der 
elektrischen Maschine, die simulativ und experimentell angewandt 
wurde. Beispielhaft wird diese Methodik am System des Radnabenan- 
triebs getestet. Ebenfalls werden die daraus entstehenden Herausforde- 
rungen aufgezeigt und Informationen über das vibroakustische Verhal- 
ten elektrischer Maschinen ableitet. Aus den gewonnenen Daten aus 
Versuch und Experiment wird eine Betriebsstrategie für Fahrzeuge her- 
geleitet, welche die vibroakustischen Eigenschaften des Systems und 
somit des Fahrzeuges verbessert. Abschließend in diesem Kapitel wird 
der vibroakustisch optimierte Steuerungsansatz im Gesamtfahrzeug an- 
gewandt um dessen Wirksamkeit zu zeigen. Kapitel 7 fasst die gesamte 
Arbeit zusammen und zeigt neue offene Fragestellungen sowie weitere 
Anwendungsgebiete auf, welche sich aus den Ergebnissen dieses Pro- 
jektes ergeben haben. Im Anhang finden sich weitergehende Informati- 
onen über die im Projekt entstanden Versuchsaufbauten und Messer- 
gebnisse. 


Diese Arbeit bezieht sich auf das Beispiel eines Fahrzeuges mit elektri- 
schen Radnabenantrieb, woraus aber auch die grundsätzlichen Entste- 
hungsmechanismen anderer elektrischer Maschinen und Optimierungs- 
möglichkeiten abgeleitet werden können. 


1 Einleitung 


Wesentliche Neuerungen oder Erweiterung in der Arbeit zum bis- 
herigen Stand der Technik: 


e Erklärung der Entstehung (Kapitel 3) und Möglichkeiten zur Re- 
duzierung der Wechselkräfte an den Anbindungspunkten im 
System Radnabenantrieb (Kapitel 5) 

e Aufbau einer multiphysikalischen Simulationskette zur Berech- 
nung des vibroakustischen Verhaltens an den Anbindungspunk- 
ten des Systems (Kapitel 4) 

o Neuartige Methode zur schnellen Berechnung der in der 
elektrischen Maschine wirkenden Kräfte in Bezug auf 
Stromoberschwingungen 

e Entwicklung eines akustisch optimierten Ansteuerverfahrens 
(Kapitel 6) 

o Entwicklung und Beschreibung einer Methodik zum Be- 
rechnen der für die Ansteuerung notwendigen Parameter 
basierend auf der in der Arbeit entwickelte Simulations- 
kette 

o Entwicklung und Beschreibung verschiedener Methoden 
zum Berechnen der für die Ansteuerung notwendigen Pa- 
rameter aus experimentellen Daten 

o Analyse des Verhaltens verschiedener Regelungsverfah- 
ren für die akustisch optimierte Ansteuerung 

o Anwendung der Ansteuerungsverfahren an verschiede- 
nen Prüfständen und am Gesamtfahrzeug 

e Methode zur Bewertung und Unterscheidung der verschiedenen 
Ursachen für die vom magnetischen Feld angeregten Vibratio- 
nen im Experiment (Kapitel 6.2.5) 


2 Grundlagen 


2.1 Radnabenantrieb 


Radnabenantriebe sind vollständig in die Felge integrierte elektrische 
Direktantriebe. Leistungselektronik mit Controller, Motor und Bremse 
befinden sich bei dem hier untersuchten Konzept innerhalb des Rades. 
Hierbei ist die Statorplatte starr mit dem Stator, Reibbremse und Rad- 
lager verbunden. Der Rotor ist über das Radlager mit der Statorplatte 
gekoppelt. Über die Radnabe ist die Felge mit dem Reifen verschraubt, 
sodass diese rotieren kann. Abbildung 2.1 zeigt den grundsätzlichen 
Aufbau des in der Arbeit untersuchten elektrischen Radnabenantriebs. 


Leistungs- 
elektronik und 
Controller Statorplatte 


O 


(Reibungsbremse, 


Radlager) Reifen und 


Rotor 
€ Stator Felge 


Abbildung 2.1: Aufbau des Radnabenantriebs 


Montiert ist der Radnabenantrieb beim zu untersuchenden Versuchsträ- 
ger an einer Koppellenkerachse. Diese besteht im Grunde aus einem 
Längslenker, welcher über eine Koppel mit dem gegeniiberliegenden 
Längslenker verbunden ist. Die Längslenker sind drehbar über die 
Längslenkerbuchse gelagert. Die Längslenkerbuchse, Dämpfer und Fe- 
der verbinden das Fahrwerk mit der Karosserie. Abbildung 2.2 verdeut- 
licht den Aufbau. 


2 Grundlagen 


Dämpfer 
Feder 


Längslenker 


x Koppel 
Langslenkerbuchse 


Abbildung 2.2: Aufbau des Gesamtsystems 


Der Motor besteht aus einer dreiphasigen, permanenterregten Syn- 
chronmaschine mit vergrabenen Magneten. Fiir das Erzeugen der 
Strangspannungen und damit der Wicklungsströme kommt ein Umrich- 
ter mit fester Taktfrequenz zum Einsatz. Die Daten des Radnabenan- 
triebes lassen sich aus Tabelle 2.1 ablesen. 


Eigenschaften Wert Einheit 
Axiale Aktivlänge 70,5 mm 
Außendurchmesser Stator 334 mm 
Rotordurchmesser 276 mm 
Luftspaltgröße 0,9 mm 
Statornutzahl 48 - 
Polpaarzahl 20 - 
Leistung (Dauer/Spitze) 33/44 kW 
Drehmoment (Dauer/Spitze) 350/400 Nm 
Gewicht des Motors inkl. Elektronik 53 kg 
Umrichterfrequenz 10 kHz 


Tabelle 2.1: Daten des Radnabenantriebs „beta“ 
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2.2 Wicklung 


2.2 Wicklung 


Im Motor des Radnabenantriebs (Version ,,beta“) kommt eine Zahnspu- 
lenwicklung zum Einsatz. Bei dieser Variante liegen Hin- und Rücklei- 
ter einer Wicklung in benachbarten Nuten. Der Vorteil dieser Wickel- 
technik ist vor allem die Möglichkeit der maschinellen Fertigung, wel- 
che im Vergleich zu den verteilten Wicklungen leichter und kosten- 
günstiger ist. Nachteilig bei diesem Verfahren sind die erhöhten Ober- 
wellenanteile der magnetischen Flussdichte im Luftspalt, was zu erhöh- 
ter Geräuschemission führen kann. Abbildung 2.3 zeigt die verwendete 
Zahnspulenwicklung. 


Wicklung W 


Abbildung 2.3: Wicklungsschema des Radnabenmotors 


11 


2 Grundlagen 


2.3 Magnete 


Als Permanentmagnete kommen Neodym-Eisen-Bor (NdFeB) Mag- 
nete zum Einsatz. Beschrieben werden können diese über die Koerziti- 
vfeldstärke, die relative Permeabilität und die Magnetisierungsrichtung. 
Für die Arbeit werden die relative Permeabilität und die Koerzitivfeld- 
stärke der Magnete für die elektromagnetische FEM Simulation als 
konstant angenommen. 


Physikalische Größe Wert Einheit 
relative Permeabilität u, 1 - 
Koerzitivfeldstärke H, -1034,5 kA/m 


Tabelle 2.2: Daten des Radnabenantriebs 


Die Magnete sind im Rotor in der Halbach-Anordnung platziert. Dies 
sorgt dafiir, dass sich der magnetische Fluss auf der Seite des Luftspal- 
tes zum Stator hin verstarkt, wohingegen dieser auf der Riickseite zum 
Rotorträger stark geschwächt ist. 


Abbildung 2.4: Anordnung der Magnete 
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2.4 Feldorientierte Regelung 


2.4  Feldorientierte Regelung 


Bei der Regelung der elektrischen Maschine des Radnabenantriebes 
kommt die feldorientierte Regelung [Sch15] zum Einsatz, welche für 
den Betrieb von permanent erregten Synchronmaschinen als Traktions- 
antrieb am weitesten verbreitet ist. Das Grundprinzip besteht darin, die 
drei Wechselgrößen der Phasenströme oder Strangspannungen in einem 
rotorfesten System mit Hilfe von zwei senkrecht aufeinander stehenden 
Ersatzgrößen (d-Achse und q-Achse) darzustellen. Abbildung 2.5 zeigt 
eine vereinfachte Darstellung der feldorientierten Regelung. Dabei liegt 
die d-Achse in Richtung der vom Dauermagneten ausgehenden magne- 
tischen Flussverkettung Ypy . Die q-Achse liegt senkrecht dazu. 


& 
NS 


Abbildung 2.5: Darstellung des Bezugsystems einer feldorientierten 
Regelung 


Zum Betrieb der Maschine werden grundsätzlich die Ströme I, und 
Iq eingeregelt, wie Abbildung 2.6 zeigt. Dabei wirkt I} hautsächlich 
drehmomentbildend, wohingegen Jy mit negativem Vorzeichen feld- 
schwächend ist, um die durch das Rotormagnetfeld induzierte Span- 
nung zu reduzieren. Sind die d- und q- Induktivitäten Lg und L, einer 
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2 Grundlagen 


elektrischen Maschine jedoch unterschiedlich, so kann Jg auch zum 
Drehmoment beitragen, wie folgender Zusammenhang zeigt. 


3 3 
M => Wpm lq +5 (La — La) ` la` lq (2.1) 


Abbildung 2.6 zeigt die schematische Darstellung einer feldorientierten 
Regelung mit dessen Regler und Transformation. Vorgegeben werden 
die rotorfesten Ersatzstromgrößen J, und I4 statt die Phasenströme 


Iy, ly und Iw. 


Stromtabellen 


"us: 


Drehzahl g 


Inv. Clarke-Park 1 


dA i Uz Leistungs- 


DA h Pag U, elektronik |4 


Moment- 
vorgabe 


Regler 


= IS 
= & 
ÀA 


Drehzahl ---—-- > 


A0 


At Rotor- 
position 
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Abbildung 2.6: Darstellung einer feldorientierten Regelung 


Häufig werden für die richtige Wahl von q- und d-Strom im Vorfeld 
berechnete Stromtabellen genutzt, welche meist nach der MTPC (,,Ma- 
ximum Torque per Current“)-Strategie [Mey10] berechnet wurden. Die 
Wahl der Verhältnisse zueinander hängt von Drehmoment und Dreh- 
zahl ab. Grundsätzlich wird versucht durch richtige Wahl im jeweiligen 
Betriebsbereich den Motor möglichst energieeffizient zu betreiben. Es 
finden sich in der Literatur viele verschiedene Ausprägungen der feld- 
orientierten Regelung, die allerdings hier nicht ausführlich behandelt 
werden sollen. Eine gute Erklärung und Überblick liefern [Nus10] und 
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2.4 Feldorientierte Regelung 


[Gem15]. Hierzu werden in dieser Arbeit auch verschiedene Rege- 
lungsansätze auf einem Motorenprüfstand für elektrische Maschinen 
und deren Reaktion auf den in der Arbeit entwickelten, akustisch opti- 


mierten Steuerungsansatz untersucht (Kapitel 6.3.4). 


Die Umrechnung der drei Wicklungsströme in das rotierende System 
wird Clarke-Park oder d/q-Transformation genannt und lässt sich fol- 


gendermaßen beschreiben. 


2T An 
2 cos(@) cos(@ — = cos(@ — Saal Ty 
[i |=; 2n i 4 (2.2) 
—sin(@) -sin(9 — =) -sin(0 — 7) 
Die inverse Transformation kann beschrieben werden durch: 
cos(0) a 
I 20 
[v _ |cos(@ =g) —sin(@ = H (2.3) 
Iw An 
os(0 — 7) -sin(0 — =, 
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3 Quellen und 
Übertragungswege 


3.1 Beschreibung der Quellen 


Geräusche und Vibrationen in elektrischen Maschinen können ver- 
schiedene Ursachen haben. Es kann zwischen mechanischen Ursachen, 
fluiddynamischen und elektromagnetischen Ursachen unterschieden 
werden, wie die nachfolgende Darstellung zeigt. 


Geräusch- und Vibrationsquellen in elektrischen Maschinen 


mechanisch 


fluid-dynamisch elektromagnetisch 


Luft-/Flüssigkühlung Magnetostriktion 


Lager 


Lorentzkräfte 


Unwucht Maxwellsche 


Grenzflächenkräfte 


Abbildung 3.1: Übersicht möglicher Geräusch- und Vibrationsquellen 
in elektrischen Maschinen (in Anlehnung an [Vij98]) 


Im Folgenden wird Bezug auf den Radnabenantrieb genommen und ob 
die dargestellten möglichen kritischen Geräuschquellen bei elektri- 
schen Maschinen auch für diese Art von System für die hauptsächlich 
vibroakustischen Störungen im Fahrzeuginnenraum verantwortlich sein 
können. 


Mechanische Quellen, die aus dem Lager oder Unwuchten resultieren, 
können ausgeschlossen werden, da die daraus entstehenden möglichen 
Ordnungen nicht mit den als störend empfundenen Ordnungen überein- 
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stimmen. Der Radnabenantrieb besitzt eine Flüssigkühlung, in der Was- 
ser mit gefrierpunktabsenkenden Zusätzen eingesetzt wird. Aufgrund 
von Vorversuchen ohne Kühlflüssigkeit können auch fluiddynamische 
Quellen ausgeschlossen werden. Mechanische und fluiddynamische 
Ursachen spielen somit beim Radnabenantrieb eine untergeordnete 
Rolle und werden daher in dieser Arbeit nicht weiter betrachtet. 


Die elektromagnetischen Ursachen haben hingegen großen Einfluss 
auf das Schwingungsverhalten des Systems. Dies belegen viele For- 
schungsberichte und Veröffentlichungen zu diesem Thema. Bereits 
[Jor50] beschrieb die an der Geräuschbildung beteiligten physikali- 
schen Vorgänge, wobei als Quelle die elektromagnetischen Wechsel- 
kräfte im Luftspalt ausgemacht wurden. Somit ist sehr wahrscheinlich, 
dass die störenden Effekte auch im Fall des Radnabenantriebs die 
Hauptursache dieser vibroakustischen Störungen sind. Aus diesem 
Grund wird in der Arbeit die Entstehung und Übertragung der elektro- 
magnetischen Kräfte genauer untersucht und erläutert. Diese im elektri- 
schen Motor wirkenden Kräfte können wiederum in Elektro-/Magne- 
tostriktion, Lorenzkräfte und Maxwellsche Grenzflächenkräfte aufge- 
teilt werden (siehe Abbildung 3.1). 


Ursachen, die dem Effekt der Magnetostriktion zugeordnet werden 
können, spielen in der Regel bei vibroakustischen Fragestellungen in 
elektrischen Maschinen eine untergeordnete Rolle [Ari01], auch wenn 
dieser Effekt in speziellen Fällen bei der Anregung von Resonanzen 
[Mel07, Kur18] auftreten kann. In dieser Arbeit wird dieser Anregungs- 
mechanismus nicht weiter betrachtet. 


Ebenso gilt dies für die Lorentzkräfte, die auf die stromdurchflossenen 
Leiter wirken, da die in den Nuten liegenden Leiter fast feldfrei sind 
[Pon10]. 


Hauptquelle der Geräusche und Vibrationen in permanent erregten Syn- 
chronmaschinen sind die Maxwellschen Grenzflächenkräfte. Diese 
und deren Auswirkungen werden in dieser Arbeit genauer beschrieben 
und untersucht. 
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3.1 Beschreibung der Quellen 


Einen Überblick zu den Zusammenhängen von auftretenden störenden 
vibroakustischen Effekten zeiget Abbildung 3.2. 


Regelung 
(Ursprung 2) 
Geometrisches 


misst stellt n A Material- 
ai Wickl - Design 
Ströme Spannung ee) € toleranzen 


Material- 


abhängig Magnetik eigenschaften 


Wicklungsströme K 
voneinander (Ursprung 1) 


Fehl- 
stellungen 
Maxwell‘sche 
Grenzflachenkrafte 


oszillierende Radial- und 
Tangentialkräfte im Luftspalt 


(Kräfte auf Rotor und Stator) 


Mechanischer 
Aufbau des Motors 


Mechanischer — Drehmoment- 
Aufbau des ungleichförmigkeiten 
angeschlossenen — Translatorische 
Systems Riittelkrafte 


störende vibroakustische 
Effekte 


Abbildung 3.2: Uberblick des Zusammenwirkens der verschiedenen 
Ursachen 


Es zeigt vor allem die Abhängigkeit der wirkenden Kräfte (Maxwell- 
sche Grenzflächenkräfte) und der Magnetik von vielen unterschiedli- 
chen Faktoren (Wicklungsschema, geometrischem Design, Materialto- 
leranzen, Materialeigenschaften, Fehlstellungen), die sich wiederum 
wechselseitig mit den Wicklungsströmen beeinflussen (Kapitel 3.1.1), 
welche ebenso von der Regelung abhängig sind. Der untere Teil der 
Darstellung zeigt die Einflüsse der mechanischen Komponenten, die im 
Kapitel 3.2.2 genauer beschrieben werden. 
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3.1.1 Betrachtung der Ursprünge 


Die vom magnetischen Feld hervorgerufenen Wechselkräfte (Maxwell- 
sche Grenzflächenkräfte) können unterschiedliche Ursprünge haben. 


e Ursprung 1: magnetische Design 
e Ursprung 2: Oberschwingungen in den Strömen 


3.1.1.1 Magnetisches Design der elektrischen Maschine als 
Quelle (Ursprung 1) 


Aufgrund des elektromagnetischen Designs einer elektrischen Ma- 
schine entstehen im Luftspalt Wechselkräfte, auch unter der Annahme, 
dass die Phasenströme keine Oberschwingungen auf der Grundwelle 
aufweisen beziehungsweise rein sinusförmig sind. Besonders großen 
Einfluss hat hier das geometrische beziehungsweise magnetische De- 
sign und deren Materialeigenschaften. Die Nuten im Statorblech ver- 
kleinern zum Beispiel den magnetischen Widerstand des nützlichen und 
drehmomenterzeugenden Magnetfeldes und erhöhen damit die Leis- 
tungsdichte der Maschine. Gleichzeitig entstehen dadurch aber im Zu- 
sammenhang mit deren Form und deren nichtlinearen Materialeigen- 
schaften, wie der Sättigung, magnetischen Oberwellen im Luftspalt. 
Diese treten auch bei rein sinusförmigen Strömen auf, was letztendlich 
zusätzliche schwingungserzeugende Kräfte innerhalb des Luftspaltes 
zur Folge hat. Gleiches gilt für die Anordnung der Magnete, der Form 
der Magnete, des Rotorblechs und deren Materialeigenschaften. Dies 
kann zusätzliche Drehmomentschwingungen sowie Anregungen in 
Normalenrichtung (bei Radialflussmaschinen also radiale Richtung) 
zur Folge haben, die vibroakustisch relevant sind. Eine präzisere Auf- 
teilung der Oberschwingungen im Drehmoment findet sich in Kapitel 
3.1.2.1 und der radialen Kraftanregung in Kapitel 3.1.2.2 wieder. 
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Fehlstellungen des Rotors gegenüber dem Stator können das Schwin- 
gungsverhalten der elektrischen Maschine signifikant beeinflussen. 
Durch diese Fehlstellung wird das gesamte magnetische Feld im Luft- 
spalt beeinflusst, sodass dieses nicht mehr symmetrisch über den Um- 
fang ist und sich die dadurch entstehenden Normalenkräfte nicht mehr 
gegenseitig aufheben. Daraus folgt das Entstehen einseitiger Zugkräfte 
auf Rotor und Stator, die wiederum die eigentliche Fehlstellung beein- 
flussen können. Diese einseitige Zugkraft ist in Richtung des kleinsten 
Luftspaltes am Umfang ausgerichtet und rotiert bei exzentrischer Fehl- 
stellung des Rotors gegenüber dem Stator mit. Diese Zugkraft besitzt 
ebenfalls eine oszillierende Komponente, welche zu zusätzlichen Ord- 
nungen im Anregungsspektrum führen kann. Ebenso entstehen durch 
das geänderte Magnetfeld und der damit verbundenen geänderten, nicht 
mehr symmetrischen Kraftverteilung im Luftspalt zusätzliche räumli- 
che Anregungsformen, die das Aufschwingen zusätzlicher mechani- 
scher Resonanzen hervorrufen kann. In bisherigen Forschungsarbeiten 
wurden hauptsächlich dynamische und statische Exzentrizitäten unter- 
sucht [Mic14]. Dass dies starke Auswirkungen auf die Vibrationen am 
Statormantel hat, zeigten [Piel2, Leell, DiT16]. Folglich kann auch 
eine Zunahme des abgestrahlten Luftschalls bei exzentrischer Lagerung 
des Rotors direkt gezeigt werden [Sch07]. Ebenso sind axiale Kippef- 
fekte des Rotors gegenüber dem Stator möglich, welche die vibroakus- 
tischen Eigenschaften weiter negativ beeinflussen können [Kra17]. 


Auch Toleranzen in den Materialeigenschaften, zum Beispiel bei der 
Magnetisierung der Permanentmagnete oder unterschiedliche magneti- 
sche Leitwerte in den Blechpaketen, können Einfluss auf die magneti- 
sche Flussdichte und somit auf das Schwingungsverhalten der Ma- 
schine ähnlich wie eine Fehlstellung haben. Ebenso wirken sich Tole- 
ranzen im geometrischen Design aus. Diese können sich beispielsweise 
auf Zahnlänge oder Zahnbreite beziehen und so das magnetische Feld 
und die damit verbundenen magnetischen Kräfte beeinflussen. 
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3.1.1.2 Schwingungen im Strom als Quelle (Ursprung 2) 


Für Traktionsantriebe mit permanenterregten Synchronmaschinen 
kommt häufig die feldorientierte Regelung zum Einsatz. Obwohl für 
einen konstanten Betriebspunkt auch ein konstanter Sollstrom eingere- 
gelt werden soll, treten ungewollte Oberschwingungen in den Strang- 
strömen auf. Gründe für diese auftretenden Schwingungen sind: 


e _ Oberschwingungsbehaftete induzierte Spannungen durch 
winkelabhängige Flüsse 

e = Nichtlinearitäten im Wechselrichter 

e _ Fehlerbehaftete Strom- und Rotorwinkelsensorik 


Durch die aufgezählten, in der elektrischen Maschine auftretenden Ef- 
fekte, die oft nicht oder nicht vollständig im Regler berücksichtigt wer- 
den, kommt es zu zusätzlichen Oberschwingungen in den Strangströ- 
men I,,, lp, Iy und somit auch in den Ersatzsteuergrößen Iq, Iq [Pel12b, 
Ricl5c, Kar16]. Diese wirken sich letztendlich auf die im Luftspalt der 
elektrischen Maschine wirkenden Wechselkräfte aus und können zu- 
sätzliche vibroakustisch spür- und hörbare Effekte hervorrufen. Ein 
Beispiel dazu findet sich in Abbildung 3.4, in der die Messergebnisse 
bei einem konstanten Betriebspunkt aufgezeigt sind. Die Aufzeichnung 
der Messwerte stammt von einem Versuch am Motorenprüfstand für 
elektrische Maschinen am elektrotechnischen Institut (ETI) des Karls- 
ruher Instituts für Technologie (KIT), welcher mit verschiedenen Be- 
schleunigungssensoren und einer Drehmomentmessstelle ausgestattet 
wurde (Abbildung 3.3). 
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Abbildung 3.3: Versuchsaufbau zur Vermessung der elektrischen Ma- 
schine (1,2: Beschleunigungssensoren; 3: Drehmomentsensor) 


In Abbildung 3.4 sind die Signale über den mechanischen Winkel 
(links) und über die mechanische Ordnung (rechts) aufgetragen. Auf- 
fällig daran sind die auftretenden Schwingungen in den Strömen J, 
(Bild oben) und /, (Bild mitte). Ebenso zeigen sich in dem Signal (Bild 
unten), welches die axiale Beschleunigung an der Statorplatte wieder- 
gibt, ähnliche Überhöhungen in den Amplituden der Ordnungen wie im 
Stromsignal. Es ist bekannt, dass es Zusammenhänge zwischen den 
Stromoberschwingungen und dem vibroakustischen Verhalten gibt. In 
dieser Arbeit wird hierzu eine einfach anwendbare Methodik in Kapitel 
6.3.5 gezeigt, um die beiden angesprochenen Ursachen beziehungs- 
weise Effekte trennen zu können. 
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Abbildung 3.4: Oberschwingungen gleicher Ordnung in den Strömen 
und den Beschleunigungen (Statorplatte z-Richtung) 
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3.1.2 Betrachtung der wirkenden Kräfte 


Die im Luftspalt zwischen Rotor und Stator wirkenden Kräfte können 
in tangentiale und radiale Kraftkomponenten, sowie Drehmomente auf- 
geteilt werden. Bei der Konzentration der Krafte auf einen Punkt pro 
Zahn, wie in Abbildung 3.5 dargestellt, entsteht ebenso ein auf die 
Zahne wirkendes Moment um den Angriffspunkt. 


Abbildung 3.5: Veranschaulichung der radialen und tangentialen 
Kräfte und Momente, die auf die Statorzähne wirken 


3.1.2.1 Wirkung der tangentialen Kräfte 


Durch tangentiale Kraftdichteschwankungen können Schwingungen 
im Drehmoment mit verschiedenen Ordnungen auftreten. Diese können 
folgendermaßen unterschieden werden: 


e Nut-/Rastmoment (Cogging) 
Diese Drehmomentschwingung kommt bei permanenterregten 
Maschinen zustande, wenn der magnetische Widerstand über 
den Rotorwinkel nicht konstant ist. Dies führt auch bei un- 
bestromten Wicklungen zu Schwingungen im Drehmoment 
[Jah96]. Die kleinste Ordnungszahl vr dieser Ungleichförmig- 
keit kann mit Hilfe des kleinsten gemeinsamen Vielfachen der 
Nutzahl N und der Polzahl 2p beschrieben werden [Jur03]. 
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Ebenso können daraus, wie auch in den nachfolgend beschrie- 
benen Ordnungen einer Drehmomentungleichförmigkeit, Viel- 
fache dieser Ordnung entstehen. 


Ur = kgV(N,2p) (3.1) 


Lastpulsation/Drehmomentrippel (Torque Ripple) 

Diese Ungleichmäßigkeit des Drehmoments resultiert aus dem 
Zusammenspiel der Felderregung durch die Bestromung der 
Wicklungen und dem magnetischen Design des Rotors [Jah96]. 
Die daraus entstehende Ordnungszahl ergibt sich aus dem 
kleinsten gemeinsamen Vielfachen der Polzahl 2p und der An- 
zahl der Wiederholungen der symmetrischen Wicklungen w um 
den Statorumfang (siehe Abbildung 4.11). 


uv, = kgV (w, 2p) (3.2) 


Weitere Drehmomentungleichförmigkeiten 

Grundsätzlich können durch Asymmetrien im magnetischen De- 
sign weitere Drehmomentungleichförmigkeiten auftreten, da 
diese Änderungen im magnetischen Feld hervorrufen, was sich 
auf die entstehenden Kräfte im Luftspalt auswirkt. Deren Ord- 
nungszahlen lassen sich wiederum über die kleinsten gemeinsa- 
men Vielfachen berechnen. In dieser Arbeit werden derartige 
Asymmetrien allerdings nicht weiter betrachtet. Gründe für der- 
artige Asymmetrien können Exzentrizitäten sowie geometrische 
Toleranzen oder Toleranzen in den Materialeigenschaften sein. 


Unter anderem durch die beschriebenen Oberschwingungen im Dreh- 
moment werden die am Motor angeschlossenen Strukturen zum 
Schwingen angeregt. 
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3.1.2.2 Wirkung der radialen Kräfte 


Ebenso kann das Vibrieren der Strukturen durch radiale Kraftdichte- 
schwankungen beziehungsweise radiale Wechselkräfte hervorgerufen 
werden. In bisherigen Forschungsarbeiten über die akustischen Eigen- 
schaften elektrischer Maschinen werden vor allem die Auswirkungen 
der Radialkräfte auf die Statorstruktur [Krol2] und anschließend auf 
das Abstrahlen von Luftschall direkt vom Statorgehäuse untersucht und 
beschrieben [Bral3a, Sch07]. Um die Anregung kritischer Eigenmoden 
in einem System zu identifizieren wird häufig die räumliche und zeitli- 
che Kraftdichteverteilung im Luftspalt über den Rotorwinkel (Abbil- 
dung 3.6 - links). mit Finite-Elemente Simulation berechnet und an- 
schließend mit Hilfe einer 2D-Fouriertransformation in eine Darstel- 
lungsform (Abbildung 3.6 - rechts) gebracht, um kritische 
Resonanzstellen und Drehzahlen abschätzen zu können. 


Kraftdichteverlaufüber den Darstellung der Kraftanregung als 
Rotorwinkel 2D-Fourierreihe 


s 
sm w 
bo 
5 200 
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2 
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Abbildung 3.6: Darstellung der Kraftdichteanregung bei 400 Nm 
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Bewegung des Rotor [°] 


Diese Darstellungen trennen räumliche und zeitliche Anteile der Anre- 
gung. So ist auf der linken Seite der Abbildung die Kraftanregung über 
6 Zähne am Umfang für eine elektrische Periode (18° mechanisch) zu 
sehen. Die Ordnungsdarstellung (Abbildung rechts) hingegen zeigt die 
Amplituden der zeitlichen Ordnungen v;, und die räumlichen Ordnun- 
gen v,, auf. Passt die Form der Kraftanregung (Raumordnung v,.,) zu 
der Raumordnung der Stator-Eigenmode v,.; (Abbildung 3.7), kann die 
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Drehzahl n, bei der diese Eigenmode in der Resonanzfrequenz angeregt 
wird, über die dazugehörige zeitliche Anregungsordnung v;, berechnet 
werden. 


_ fe 
n= 


3.3 
Ui (3.3) 


Aus dem Spektrum der Anregung in Abbildung 3.6 (rechts) ist ersicht- 
lich, dass relevante Anregungsformen beziehungsweise Anregungsord- 
nungen bei v,a = 0; v;4 = 8; va = 16;... und weiteren Vielfachen 
von acht auftreten. 


fe =524 Hz fe =525 Hz fs =1177 Hz fx =1182 Hz 
Oval (vu, = 2) Oval (v,g = 2) Dreieck (v,z = 3) Dreieck (u,~ = 3) 


Abbildung 3.7: Relevante Strukturmoden fiir Radialkraftdichte 


Vergleicht man diese mit dem hier relevanten Frequenzbereich erkennt 
man, dass hier keine entsprechenden Strukturmoden existieren, die 
durch eine derartige Kraftform angeregt werden. Dieser Anregungsme- 
chanismus spielt also in diesem System nur eine untergeordnete Rolle 
und wird aus diesem Grund in dieser Arbeit nicht naher untersucht. 


Ein weiterer wichtiger Effekt ist das Entstehen von Wechselkräften auf- 
grund Asymmetrien in den mechanischen Übertragungswegen. 
Dadurch können auch bei symmetrischen radialen Kräften, die auf Sta- 
tor oder Rotor wirken, resultierende Kräfte am Anbindungspunkt der 
elektrischen Maschine entstehen [Str17b]. Da dieser Effekt auch beim 
System Radnabenantrieb zu erwarten ist, wird dieser genauer in Kapitel 
3.2.2 beschrieben. 
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3.1.3 Modellierung der Quellen 


3.1.3.1 Berechnung der Kräfte aus der magnetischen Fluss- 
dichte 


Die magnetischen Kräfte (Maxwellsche Grenzflächenkräfte) wirken im 
Luftspalt auf Rotor und Stator. Sie können aus der magnetischen Fluss- 
dichte B berechnet werden, welche wiederum von drei Faktoren abhän- 
gig ist: 


e Winkel im statorfesten System 6, 
e Lage des Rotors beziehungsweise Rotorwinkel 6 
e Bestromung der Statorwicklungen [,,, Iy, hy bzw. Ig, lą 


Am Ubergang vom Blechpaket zum Luftspalt zwischen Rotor und Sta- 
tor kann diese in einer zweidimensionalen Darstellung in radiale und 
tangentiale Flussdichte aufgeteilt werden. Die auf Stator und Rotor wir- 
kenden Kraftdichten f, und fẹ können aufgrund der großen Unter- 
schiede der Permeabilitätszahlen zwischen Eisen und Luft (pure « 
Ukisen) vereinfacht über den Maxwellschen Spannungstensor beschrie- 
ben und folgendermaßen berechnet werden [Zbo1 1]: 


il 
fä 1 (5 (B? - Be) 
"T B,By ae 


Aus dieser Beschreibung können lokale Kraftdichten im Luftspalt in 
Abhängigkeit des Statorwinkels berechnet werden. 


Abbildung 3.8 zeigt beispielhaft die radiale Kraftdichte bei einem Be- 
triebspunkt für einen bestimmten Rotorwinkel. Umgerechnet wird die 
Kraftdichte auf eine tangentiale und radiale Kraft, die jeweils auf einen 
Zahn am Stator oder einen Pol am Rotor wirken. 
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N 
Fri =f" frdsy + Sax 


SLi+1 


Abbildung 3.8: 2D-Darstellung der Kraftdichte im Luftspalt 


Die tangential und radial angreifenden Kräfte F,; und F,; an den 
Statorzähnen i können mit Hilfe der Strecke im Luftspalt s, über einen 
Zahn (s};+1 — Szi) und der axial wirkenden Strecke Say berechnet wer- 
den. 


Es gilt somit: 


SLi+1 

Fri = | feds, ' Sax (3.5) 
SLi 
SLi+1 

Fri = | frdSL * Sax (3.6) 
SLi 


Durch die Verschiebung der radialen Krafte auf den Mittelpunkt des 
Zahnes entsteht aufgrund asymmetrischer Verteilung der Kraftdichte 
ein Drehmoment, welches am Zahn i angreift. Durch die sehr kleine 
Krümmung des Luftspaltes über einen Zahn kann diese vernachlässigt 
beziehungsweise die Strecke sz; bis s,;,,als Gerade angenommen wer- 
den. 
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Sodass gilt: 
SLi+1 
Mix [fe (i= 82) dsr Sax 3.7) 
SLi 


Untersuchungen [Kes16] haben gezeigt, dass dieses Moment sehr klein 
ist und die mechanischen Ubertragungswege des Momentes zu resultie- 
renden Effekten sehr schwach ausgeprägt sind. Somit ist dieser Einfluss 
auf das vibroakustische Verhalten sehr gering. 


Die gleiche Methodik wird für die auf den Rotor wirkenden Kräfte an- 
gewendet. 


3.1.3.2 Modellierungsmöglichkeiten der magnetischen 
Flussdichte 


Wie bereits beschrieben, ist die magnetische Flussdichte im Luftspalt 
die Grundlage für die Berechnung der wirkenden Kräfte. Durch ver- 
schiedene Methoden kann die magnetische Flussdichte im Luftspalt be- 
rechnet werden. Im Folgenden wird kurz auf die verschiedenen Model- 
lierungsmöglichkeiten eingegangen. 


Analytische Modelle 


Durch analytische Berechnungsverfahren kann die magnetische Fluss- 
dichte B berechnet werden. Für erste analytische Modelle können ver- 
einfachte Bedingungen angenommen werden. Hierbei geht man von ei- 
ner Statoroberfläche ohne Nuten und mit punktförmigen Leitern sowie 
von einem gleichbleibenden Luftspalt aus. Auch Sättigungseffekte wer- 
den hierbei nicht berücksichtigt. Daraus lässt sich auf einfache Weise 
die magnetische Flussdichte berechnen und damit auch Kräfte in radia- 
ler Richtung abschätzen. Mit Hilfe dieser Methode können bereits erste 
Aussagen über die entstehenden Ordnungen getroffen werden. 


Eine Weiterentwicklung dieser analytischen Methode ist die konforme 
Abbildung [Rab96, He15], womit auch Luftspaltunterschiede zwischen 
Rotor und Stator und somit die Geometrie der Zähne und des Rotors 
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abgebildet werden können. Ebenso gibt es Ansätze, welche Sättigungs- 
effekte berücksichtigen [Chel6]. Für erste Abschätzungen können 
diese Modelle gut aufzeigen, welche Ordnungszahlen sich für Oberwel- 
len sich im Luftspalt ergeben, woraus sich auch die Anregungsordnun- 
gen gut ableiten lassen. Diese Modelle zu erweitern, sodass daraus auch 
die Amplituden berechnet werden können, ist sehr anspruchsvoll und 
aufwändig, da geometrische und materialabhängige Eigenschaften in 
mathematisch beschreibbare Funktionen dargestellt und aufgelöst wer- 
den müssen. So muss dieses Modell häufig auch schon bei kleinen Än- 
derungen in der Geometrie neu hergeleitet werden. 


FEM-Modell 


Genaueren Aufschluss über die im Luftspalt wirkenden Kräfte kann vor 
allem eine elektromagnetische FEM-Simulation liefern. Ein direkter 
Vergleich mit dem Experiment beziehungsweise eine direkte Validie- 
rung ist mit dem im Forschungsprojekt zur Verfügung stehenden Mit- 
teln nicht möglich gewesen und stellt sich grundsätzlich als sehr 
schwierig und aufwändig dar. Diese Kräfte können lediglich indirekt 
über die magnetische Flussdichte oder deren Auswirkungen (z.B. 
Strukturschwingung) gemessen werden. Aufgrund der Möglichkeit mit 
dieser Methode grundsätzlich alle Effekte modellieren zu können, wird 
diese als Referenz zum Vergleich der verschieden Modellierungsmög- 
lichkeiten und für die spätere Berechnung der wirkenden Kräfte ge- 
nutzt, da mit dieser Methode eine gute Annäherung zum Experiment zu 
erwarten ist. Nachteilig bei dieser Methode ist die dafür notwendige 
Rechenzeit, auch schon bei 2D-Modellen. Hinzu kommt, dass bei einer 
Erweiterung des Modells auf 3D eine circa 50-fach höhere Rechenzeit 
zu erwarten ist. 


Ein Vergleich beim Elektromotor des Radnabenantriebs zwischen der 
analytischen Herleitung über der konformen Abbildung der magneti- 
schen Flussdichte nach [Zar06] und der elektromagnetischen FEM-Be- 
rechnung zeigt große Unterschiede [He15] und damit auch große Un- 
terschiede in der daraus resultierenden Kraftdichteverteilung im Luft- 
spalt. 
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In der elektromagnetischen FEM-Berechnung werden im Vergleich zur 
angewendeten konformen Abbildung mehr Einflussgrößen berücksich- 
tigt, sodass für nachfolgende Berechnungen das FEM-Modell als 
Grundlage verwendet wird. Weitere Vergleiche zu Berechnungsmetho- 
den der Anregenden Kräfte finden sich in [Pil18]. 


Feldrekonstruktionsmethode 


Eine auf die elektromagnetische FEM aufbauende Methodik ist die Fel- 
drekonstruktionsmethode [Kho08, Zhu06]. Bei dieser Vorgehensweise 
wird die Berechnung des magnetischen Feldes mit Hilfe von FEM-Be- 
rechnungen für verschiedene Betriebspunkte durchgeführt, extrahiert 
und abgespeichert. Im einfachsten Fall wird hier beispielsweise das 
Feld der elektrischen Maschine mit und ohne Bestromung simuliert. 
Ohne Berücksichtigung von Nichtlinearitäten kann für das Feld somit 
das Superpositionsprinzip angenommen werden, wodurch sich das Feld 
in seine Abhängigkeiten zerlegen lässt. Dieses kann dann im Anschluss 
für beliebige Betriebspunkte oder auch Stromformen aus den bereits 
gespeicherten Daten ohne zusätzliche FEM-Rechnung rekonstruiert 
werden. Vorteile bietet diese Methode in Bezug auf Rechenzeit, vor al- 
lem wenn, wie in dieser Arbeit, viele verschiedene Stromformen ge- 
rechnet werden sollen. Diese reduziert sich mit dieser Methode auf ei- 
nen Bruchteil im Gegensatz zur FEM. Besonders gut geeignet ist diese 
Vorgehensweise für Motoren mit Oberflächenmagneten. Eine Erweite- 
rung dieser Methode auch für Motoren mit vergrabenen Magneten, wie 
sie im Radnabenantrieb vorkommen, findet sich in [Öze16]. Nachteilig 
sind jedoch die erforderlichen Vorberechnungen, aufgrund vieler nicht- 
linearer Effekte, welche bei derartigen Maschinen mitberücksichtigt 
werden müssen, um eine ausreichende Genauigkeit zu erreichen. 


Kraftrekonstruktionsmethode 


Für diese Arbeit wurden die Vorteile der FEM-Methode und der an- 
schließenden Reduzierung der Informationen auf die wirkenden Kräfte 
genutzt. Durch Vorausberechnung verschiedener Betriebspunkte ist es 
möglich mit nur wenigen Daten einen großen Betriebsbereich der 
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elektrischen Maschine abzudecken. Da die vibroakustischen Eigen- 
schaften im Mittelpunkt der Untersuchungen der Arbeit stehen, wurden 
die wirkenden Kräfte in den Frequenz- beziehungsweise Ordnungsbe- 
reich umgerechnet, sodass auch Interpolation zwischen den Betriebs- 
punkten möglich ist. Um auch die durch Stromoberschwingungen ver- 
ursachten Effekte nachbilden zu können wurde das Übertragungsver- 
halten zwischen den Strömen J, und J, analysiert und modelliert. Dies 
wird in Kapitel 3.1.3.3 genauer beschrieben. 


Die verschiedenen Methoden zur Berechnung können mit Tabelle 3.1 
zusammengefasst werden. Dabei bezieht sich die Rechenzeit auf die 
Zeit, die zur Berechnung eines neuen Betriebspunktes oder einer geän- 
derte Stromform notwendig ist. 


Berechnungsmethode Genauigkeit Rechenzeit 
Analytisch ungenau sehr kurz 
Finite-Elemente-Methode Sehr gut sehr lange 
Feldrekonstruktionsmethode mittel kurz 
Kraftrekonstruktionsmethode gut sehr kurz 


Tabelle 3.1: Zusammenfassung der Methoden 
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3.1.3.3 Kraftrekonstruktion 


Abbildung 3.9 zeigt den Aufbau des elektromagnetischen Modells, wel- 
ches auf der Kraftrekonstruktion basiert. 


PRot > Radialkräfte 


—>|  Tangentialkräfte 


tangentiale Kraft 
Iq > 
Gizi © am i-ten Zahn 
Id 
Gizi 
Tq radiale Kraft 
Gz am i-ten Zahn 
Id 
Gri 


Abbildung 3.9: Aufbau des elektromagnetischen Modells zur Berech- 
nung der wirkenden Kräfte auf Statorzähne 


Als Eingangsparameter für das Modell dienen der Rotorwinkel, Dreh- 
moment und Drehzahl sowie die Oberschwingungen der Ströme. Die 
konstante Bestromung Iq und J, für einen statischen Betriebspunkt ruft 
in Zusammenhang mit dem Permanentmagneten harmonische Kräfte an 
der Zahnoberfläche hervor. Die vorab mit FEM gerechneten Betriebs- 
punkte können extrahiert und mit Hilfe von Tabellen, wie Abbildung 
3.11 zeigt, dargestellt und abgespeichert werden. Betriebspunkte zwi- 
schen den Stützstellen werden interpoliert. Ebenso rufen die harmoni- 
schen Anteile der Ersatzströme I, und J, zusätzliche harmonische An- 
teile im Luftspalt hervor, welche mit Hilfe der tangentialen Übertra- 


: I : : I : Pr 
gungsfunktionen G,},, Gig, sowie der radialen G.7,, G/; auf die Zähne 
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i beschrieben werden können. Die harmonischen Anteile und deren da- 
zugehörigen Übertragungsfunktionen können somit je nach Stromform 
beliebige Ordnungen in die Kraftanregung einleiten. 


Einfluss der Lage des Luftspaltes für Kraftdichteberechnung 


Da die Kräfte aus der magnetischen Flussdichte im Luftspalt berechnet 
werden, ergeben sich mehrere mögliche Radien der abzugreifenden 
Größe. Abbildung 3.10 zeigt diese Unterschiede auf. Für diese Arbeit 
wird der Radius (139,89 mm) nahe des Stators gewählt, da die auf den 
Stator wirkenden Kräfte den wesentlich größeren Einfluss auf das 
vibroakustische Verhalten haben. 


—— 139,15mm 
—— 139,30mm 
—— 139,45mm 
—— 139,60mm 
—— 139,75mm 
——139,89mm 


mag. Flussdichte [T] 


Variation der 
Lufstpaltradien 


0 10 20 30 40 
Statorwinkel [°] 


Abbildung 3.10: Einfluss des Radius auf die magnetische Flussdichte 
in radialer Richtung im Betriebspunkt 100 Nm/100rpm 


Ergebnisse FEM bei statischen Betriebspunkten mit konstanter 
Bestromung 


Dargestellt sind nachfolgend in Abbildung 3.11 die auf den Stator und 
Rotor wirkenden Kräfte und Momente und deren Abhängigkeiten von 
Drehzahl und Drehmoment für den Radnabenantrieb. 
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Abbildung 3.11: Übersicht der am stärksten ausgeprägten Ordnungen 
am Stator 
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Ergebnisse FEM bei statischen Betriebspunkten mit harmonischen 
Anteilen in den Stromgrößen 


Um auch die Auswirkung der harmonischen Stromschwankungen auf 
die Kräfte mit Hilfe der in Abbildung 3.9 dargestellten Methode rekon- 
struieren zu können, wurden die dabei entstehenden Kräfte in verschie- 
denen Betriebspunkten mit Hilfe gezielt aufgeprägter Harmonischen in 
den Strömen J, und Iq analysiert und Übertragungsfunktionen daraus 
abgeleitet. Hierbei wurden als harmonische Testsignale I„arm mit den 
Amplituden â , Phasenlage y und Ordnung v für die Ströme I, und Ig 
erzeugt. Als /konst werden die für den jeweilig untersuchten Betriebs- 
punkt passenden Parameter für J, und I4 bezeichnet. 


Inarm = Ikonst + @ sin(v "0 + y) (3.8) 


Das daraus dargestellte Amplitudenverhältnis a berechnet sich aus der 
Differenz der darzustellenden Größen mit (x,) und ohne (xo) harmo- 
nischem Testsignal im Verhältnis zur Amplitude A der aufprägten 


Schwingung. 
Ku 
a = = | (3.9) 
Ebenso kann die Phasenverschiebung 6 ermittelt werden: 
ô = y —arg(x, — Xo) (3.10) 
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Die Übertragungseigenschaften für die harmonischen Ordnungen kön- 
nen grundsätzlich abhängig von folgenden Faktoren sein: 


w N w 
Ss a Ss 
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Abbildung 3.12: Abhängigkeit von Stromharmonischen in I, 


Beispielhaft wird in Abbildung 3.12 die Abhängigkeit des Drehmo- 
ments von harmonischen Oberschwingungen im Strom dargestellt. Die 
Ordnungsauflösung der Darstellung beträgt 20, sodass erkennbar ist, 
dass sich hauptsächlich die aufgeprägte Harmonische im Spektrum des 
Drehmomentes wiederfindet. In den oberen Diagrammen zeigt sich die 
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Abhängigkeit des Drehmoments von der Stromschwingungsamplitude 
am Beispiel der 100. Ordnung. Diese wurde gewählt, da hier keine Be- 
einflussung durch überlagerte Drehmomentschwingungen zu erwarten 
sind. Es kann die Änderung des Drehmoments im Ordnungsbereich 
(links) und das sich daraus ergebende Amplitudenverhältnis a (Glei- 
chung 3.9) für die aufgeprägte Schwingung (rechts) betrachtet werden. 
Darunter sind äquivalent zu den oberen Darstellungsformen die Aus- 
wirkungen auf die Tangentialkräfte dargestellt. Die Nebenordnungen 
60. und 140. Ordnung wirken auf die Zähne, liefern aber keinen Beitrag 
zum Drehmoment, wie in Abbildung 3.14 zu erkennen ist. Dies liegt 
daran, dass diese Ordnungen in der 2D-FFT Betrachtung der Kraft- 
dichteverteilung (Abbildung 3.13) nicht auf der 0. Raumordnung lie- 
gen. 


Darstellung der Kraftanregung als 
2D-Fourierreihe 


10* 

10? 

10° 
-96 -48 0 48 96 


anregende räumliche Ordnung v,4 


] 


N 
m? 


anregende zeitliche Ordnung v;, 
Kraftdichte [- 


Abbildung 3.13: Darstellung der tangentialen Kraftanregung als 2D- 
Fourierreihe 


Eine strukturdynamische Anregung von Resonanzen durch diese Kräfte 
sind ebenfalls nicht zu erwarten, da Eigenformen, in denen gegenläufi- 
ges Schwingen der Statorzähne in Tangentialrichtung auftritt, erst bei 
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sehr hohen Frequenzen möglich sind. Somit können diese Nebenord- 
nungen für spätere Berechnungen vernachlässigt werden. 


Diese Auswertung wurde für verschiedene Variationen durchgeführt, 
um mögliche Einflüsse erkennen zu können. Beispielhaft zeigt hier Ab- 
bildung 3.14 die Einflüsse harmonischer Stromoberschwingungen auf 
I, auf die Amplitudenverstärkung a ausgehend von folgenden Betriebs- 
punkteinstellungen: 


e Betriebspunkt: M = 100 Nm /n = 100 rpm 
e Stromoberschwingung: â = 10 A/y = 0° /v = 220 


Dabei zeigen sich folgende Zusammenhänge: 


e Sehr schwache Abhängigkeit von der Amplitude 
> Annahme eines konstanten Verstärkungsfaktors und Phasen- 
verschiebung möglich 
e Sehr schwache Abhängigkeit von der Phasenlage 
> Annahme eines konstanten Verstärkungsfaktors und einer 
konstanten Phasenverschiebung möglich 
e Schwache Abhängigkeit von der Ordnung 
> Annahme eines konstanten Verstärkungsfaktors und Phasen- 
verschiebung möglich 
e Abhängigkeit vom Konstantanteil der Ströme 
> Da vor allem die tangentialen Kräfte starke Auswirkungen 
auf das vibroakustische Verhalten zeigen, können die in Ta- 
belle 3.2 dargestellten Parameter für die Simulation genutzt 
werden 
e Abhängigkeit vom betrachteten Zahn 
> Es werden gerade und ungerade Zähne unterschieden 
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Abbildung 3.14: Übersicht der ermittelten Parameter für die Übertra- 
gungsfunktion bei harmonischer Anregung in I, 
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Gleiche Untersuchungen wurden mit harmonisch überlagerten Strömen 
in Jq durchgeführt und auch die auf den Flächen der Pole angreifenden 
Kräfte wurden ausgewertet. Die gemittelten Parameter für die Übertra- 
gungsfunktionen G/4,, G/4,, G4, und G/4; zur Kraftrekonstruktionsme- 
thode mit ihren Parametern (Amplitudenverstärkung a und Phasenver- 
schiebung 6), welche in Abbildung 3.9 dargestellt sind, finden sich in 
Tabelle 3.2 wieder. Dabei wird zwischen geraden und ungeraden Zäh- 
nen unterschieden. Grundsätzlich ist es möglich die Methodik im Hin- 
blick auf weitere auch schwächere Abhängigkeiten bei Bedarf zu er- 
weitern. Für die folgende Arbeit wird sich jedoch auf die dargestellten 


Parameter beschränkt. 


Amplituden- Phasen- 
Verstärkung a | verschiebung ô 


tzi ungerade 0,38 0 


cid gerade 0,05 180 
tal ungerade 0,01 0 
ga gerade 0,05 180 
rZi ungerade 0,20 0 
Gia. gerade 0,55 0 

zZ ungerade 0,55 0 


Tabelle 3.2: Parameter für die Rekonstruktion der Kräfte am Stator 


Eine analoge Vorgehensweise kann für die wirkenden Kräfte auf die 
Rotorpole durchgeführt werden. 
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Eine Validierung der Kraftrekonstruktionsmethode wird in Abbildung 
3.15 gezeigt. Hierzu wurden im Strom J, unterschiedliche Amplituden 
â bei unterschiedlichen Ordnungen v eingeprägt. 


fab 100N / 100rpm; 6A / 60. Ordnung joo 100N / 100rpm; 6A / 240. Ordnung 
— FEM — FEM 
—-Kraftrekonstruktion —-Kraftrekonstruktion 
A 
50 N á 50 \ \ 
= \ \ = \ \ 
Z \| & \ 
& 0 \ E 0 
£ \| & J / 
En} / x 
N / / 
-50 Ñ N -50 
\ J J j 
\S =, Pa 
-100 -100 
0 5 10 15 0 5 10 15 
mechanischer Winkel [°] mechanischer Winkel [°] 
fon 100N / 100rpm; 10A / 220. Ordnung iò 200N / 100rpm; 10A / 100. Ordnung 
— FEM — FEM 
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Abbildung 3.15: Darstellung der tangentialen Krafte auf den 1. Zahn 
zum Vergleich der Kraftrekonstruktionsmethode mit der FEM bei ver- 
schiedenen Parametern 


Die Ergebnisse zeigen in verschiedenen Betriebspunkten gute Uber- 
einstimmungen mit der FEM, sodass diese Methode fiir die weiteren 
Berechnungen verwendet werden kann. Die Radialkrafte zeigen eben- 
falls eine sehr gute Ubereinstimmung. 
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3.2 Beschreibung der Übertragungswege 


3.2.1 Überblick der Übertragungswege 


Aufgrund der komplexen Übertragungswege im Verbund Motormecha- 
nik und Fahrwerk werden in den nachfolgenden Absätzen diese Wege 
auf vereinfachte aber anschauliche Weise beleuchtet und erklärt, wie 
sie bei einer Antriebstopologie mit elektrischem Radnabenantrieb im 
Automobil vorkommen. Dabei sind nicht die Antriebskräfte im Fokus, 
sondern die oszillierenden Kräfte, welche im Luftspalt des Motors wir- 
ken und zu vibroakustisch störenden Effekten führen können. 


Da in diesem Kapitel der Motor des Radnabenantriebs im Fokussteht, 
wird nicht das übliche Fahrzeug-Koordinatensystem (KOS) verwendet, 
sondern ein auf die elektrische Maschine ausgerichtetes Koordinaten- 
system. Die Zusammenhänge zwischen den Koordinatensystemen sind: 


e x-Achse Motor-KOS = x-Achse Fahrzeug-KOS (Längsachse) 
e y-Achse Motor-KOS = z-Achse Fahrzeug-KOS (Vertikalachse) 
e z-Achse-Motor-KOS = -y-Achse Fahrzeug-KOS (Querachse) 


3.2.1.1 Wirkende Radialkräfte auf die Rotoreinheit 


Abbildung 3.16 zeigt dabei den Weg der Kräfte beziehungsweise des 
Körperschalls von wirkenden Radialkräften auf die Rotoreinheit zu den 
Anbindungspunkten. Aufgrund der vor allem auf die Rotorpole wirken- 
den radialen Kräfte (1.) entstehen am Anbindungspunkt des Rotors, 
also am Radlagerzapfen (2.) resultierende Kräfte. Grund dafür können 
Exzentrizitäten oder mechanische und magnetische Asymmetrien im 
Rotor sein. Die sich mit dem Rotor drehenden resultierenden harmoni- 
schen Kräfte Frot beliebiger Ordnung v beziehungsweise Winkelge- 
schwindigkeiten vw übertragen sich dabei mit folgender Beziehung zu 
harmonischen Kräften auf den Radlagerzapfen Fp,, (x-Richtung) und 
Frry (y-Richtung). 
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(2) = Fao ‘eee (3.11) 


Frry cos(wt) 


Mit Frot = Aroesin(vwt) 


| = AporSin(vwt) aa (3.12) 


Fry cos(wt) 


ae 1 on — 1)wt) — cos((v + so 


a, sin((v — 1)wt) + sin((v + Dat 


(3.13) 


F RLy 


Ebenso entstehen aufgrund der Radialkräfte am Rotor im Zusammen- 
spiel mit dessen Geometrie resultierende Kräfte in z-Richtung (siehe 
Kapitel 3.2.2.2). Diese werden über den Längslenker (3.) an die Karos- 
serie (4.) über Feder, Dämpfer und Längslenkerbuchse weitergeleitet. 


Karosserie 


4. Kräfte auf 
Karosserie — 

| — 2. Kräfte auf Radlager 

und über Statorplatte 


__— 1. Radialkräfte auf Rotor 


Karosserie | 


Gis 
7 7 
3. Kräfte über 
Längslenker 


Radlagerzapfen 
Rotor 

Felge 

Reifen 

Straße 


Abbildung 3.16: Übertragungsweg Rotor radial 


3.2.1.2 Wirkende tangentiale Kraftdichte auf Rotor 


Die in tangentialer Richtung auf den Rotor wirkenden Kräfte durchlau- 
fen zwei Übertragungspfade (Abbildung 3.17). Das schließlich im Ro- 
tor aufgebrachte Drehmoment (2.) wird über die Felge (3.) und den Rei- 
fen (4.) danach auf die Fahrbahn übertragen. Dort wird das Drehmo- 
ment in eine translatorisch wirkende Kraft (5.) umgewandelt. Die 
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dadurch entstehende translatorische Gegenkraft (6.), wird wieder zu- 
rück über den Reifen (7.), Felge (8.), Rotor (9.), Radlager/Statorplatte 
(10.), an den Längslenker (11.) des Fahrwerks weitergeleitet. Dort wer- 
den die Wechselkräfte über die Feder- und Dämpfereinheit, hauptsäch- 
lich aber über die Längslenkerbuchse an die Karosserie übertragen. 
Grundsätzlich werden diese oszillierenden Kräfte durch den sehr „wei- 
chen“ Übertragungspfad schon bei niedrigen Geschwindigkeiten und 
Ordnungen entkoppelt. 


LT tialkräfte im Rotor — 
OERE E a - 3. Drehmoment über Felge 


2. Drehmoment über Rotor _ 


> 4. Drehmoment über Reifen 
Rotor 
Felge 
Reifen 


5. Kraft in Reifenaufstandsfläche 
Straße 


Karosserie 


11. Kräfte über Längslenker 


12.2 Beschleunigungskraft 
in y-Richtung Radlagerzapfen 
| Längs- 


. | lenker 
Karosserie 


10. Kräfte über Radlager und 
Statorplatte 


9. Kräfte über Rotor 
. 8. Kräfte über Felge 
12.1 BescMeunigungskraft in 


x-Richtung 7. Kräfte über Reifen 


6. Kräfte in 
Reifenaufstandsfläche 


Abbildung 3.17: Übertragungsweg Rotor tangential 


3.2.1.3 Wirkende radiale Kraftdichte auf Stator 


Abbildung 3.18 zeigt den Übertragungsweg der auf den Stator wirken- 
den oszillierenden Radialkräfte. Diese verursachen zum einen resultie- 
rende Kräfte bei der Anbindung der Statorplatte zum Fahrwerk (2.1). 
Zum anderen führen die Wechselkräfte zu einer Bewegung der 
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Statoraußenwand, wodurch Schall (2.2) abgestrahlt wird. Die Wechsel- 
kräfte an der Befestigung der Statorplatte zum Fahrwerk werden 
schließlich über den Längslenker (3.) und dessen Anbindungspunkte 
über Feder, Dämpfer und Längslenkerbuchse an die Karosserie weiter- 
gegeben. 


Karosserie 


2.2. Schallabstrahlung durch 


4. Kräfte auf Bewegung der Statorstruktur 


Karosserie ~ i 
| = 1. Radialkräfte auf Rotor 


=\——— 2.1 Kräfte auf Anbindung 
der Statorplatte zum 
Längslenker 


Statorplatte 


za 
Karosserie | 
| 


= Kräfte über Stator 
Längslenker 


Abbildung 3.18: Übertragungsweg Stator radial 


3.2.14 Wirkende tangentiale Kraftdichte auf Stator 


Tangential oszillierende Kräfte (1.), die auf den Stator wirken, verursa- 
chen ein Drehmoment (2.), welches über den Längslenker Kraft- 
schwankungen an der Längslenkerbuchse hervorrufen (3.), wie Abbil- 
dung 3.19 zeigt. Diese werden an die Karosserie weitergeleitet. 


Karosserie 


__— 1. Tangentialkräfte auf Rotor 


4. Kräfte auf — 


Karosserie __— 2. Drehmoment über 


Stator/Statorplatte auf 


karoseri 3 
arosserie Längslenker 


Statorplatte 
3. Kräfte über 


A Stator 
Längslenker 


Abbildung 3.19: Übertragungsweg: Rotor tangential 
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3.2.2 Modellierung der Ubertragungswege 


Die für das mechanische Modell notwendigen Steifigkeits- und Mas- 
senmatrizen wurden mit Hilfe der Finiten Elemente Methode (FEM) 
erstellt. Dazu wurden für jede Einheit (Stator, Rotor, Fahrwerk) verein- 
fachte Annahmen zur Modellerstellung getroffen. 


3.2.2.1 Statoreinheit 


Das Gehäuse bildet die äußere Hülle des Stators. Im Inneren ist das 
Statorblech mit 48 Zähnen zu finden. Im Statorblech befinden sich die 
Wicklungen, welche um dessen Zähne gewickelt sind. Der Kühlkörper 
für den Radnabenantrieb besitzt gefräste Kanäle für die Kühlmittelfüh- 
rung. Der Radzapfen bildet mit den Wälzkörpern und der dazugehöri- 
gen Radnabe das Radlager, welches den Rotor drehbar lagert. Verbun- 
den ist der Radzapfen mit der Statorplatte, welche dann an der Flansch- 
platte des Fahrwerks befestigt wird. 


Abbildung 3.20: Aufbau Statoreinheit 


Die Statorplatte dient einerseits zur Verbindung der restlichen Stato- 
reinheit mit dem Rotor, muss aber auch die Möglichkeit der Kabel- 
durchführung für Sensoren, Seildurchführung für die Bremse sowie die 
Befestigung der Bremsanlage bieten. Aus diesem Grund ist die Geo- 
metrie der Statorplatte asymmetrisch konstruiert. 
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Folgende Annahmen wurden am FE-Modell getroffen, welche durch 
den Vergleich der Statoreinheit anhand einer experimentellen Schwing- 
formanalyse verifiziert wurden [Mal14]. 


== Modell V1 
2500 
= Modell V2 
'N2000 
E 
N = Modell V3 
gio 
31000 = Experiment 
> 
500 = Modell V4 
(angepasst) 
0 


(2,0) l (2,1) i (3,0) . (3,1) i (4,0) l (4,1) 
Modenbezeichnung 


Abbildung 3.21: Vergleich der berechneten mit den experimentell be- 
stimmten Frequenzen [Mal14] 


Für das Blechpaket (BP) wurde ein Material mit transversalisotropen 
Eigenschaften verwendet [van12]. Die Vergussmasse wurde durch Er- 
satzmassepunkte modelliert und die Steigung der Kühlnuten im Kühl- 
körper wurde ebenfalls vernachlässigt. 


Modalanalyse 


Mit Hilfe der FEM-Simulation können die Eigenformen und die Eigen- 
frequenzen auch mit geänderten Randbedingungen des Stators be- 
stimmt werden, welche für das Schwingungsverhalten und den Wech- 
selkräften zwischen Elektromotor und Fahrwerk entscheidend sind. Als 
Randbedingung für die Schwingungsanalyse in der FE-Simulation 
wurde eine sehr hohe Federsteifigkeit c gewählt, an welcher die Stator- 
platte in der Realität mit dem Fahrwerk verbunden ist, um so vom Fahr- 
werk unabhängige Ergebnisse zu erzeugen. Das andere Ende dieser Fe- 
der wurde fest mit der Umgebung verbunden, sodass sich dieser Punkt 
nicht bewegen kann. Mit diesen Randbedingungen können einerseits in 
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der Modalanalyse die Eigenfrequenzen und die Eigenformen relativ re- 
alitätsnah bestimmt werden. Andererseits kann anhand der Verschie- 
bung x des Anbindungspunktes des Stators zur Feder die Kräfte Fp ein- 
fach ermittelt werden, die an das Fahrwerk weitergeleitet werden. 


Mode 
Frequenz [Hz] 
Bezeichnung 


M1 M2 M3 M4 M5 
70 84 135 175 425 
Kippen um x Kippen um y Verschiebung z Torsion um z Verschiebung x 


M6 
437 
Verschiebung y 


M9 M10 
525 980 1098 
Oval Membran 1 Membran 2 


M13 M15 MI6 
1182 1247 1324 
Membran 2 Dreieck Dreieck Membran 3+ Membran 3+ 


Abbildung 3.22: Modalanalyse Statoreinheit 


In Abbildung 3.22 sind die wichtigsten ermittelten Eigenfrequenzen 
und Eigenformen des Stators dargestellt. Neben den Eigenformen, bei 
denen der Stator gegenüber des Anbindungspunktes kippt (M1, M2), 
sich verschiebt (M3, M5, M6) oder sich dreht (M4) und den Eigenfor- 
men des Stators selbst (M7, M8, M12, M13), gibt es eine Vielzahl von 
Eigenmoden, in denen sich vor allem die Statorplatte bewegt, welche 
hier als Membranmoden bezeichnet werden. Diese unterscheiden sich 
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in der Anzahl der sich ausbildenden Maxima innerhalb der Statorplatte. 
In den ersten Moden bilden sich hier noch klar erkennbare Muster aus, 
die mit einem Maximum (M9) beziehungsweise zwei Maxima (M10, 
M11) beschrieben werden können. Ab drei Maxima (M16) ist aller- 
dings kaum noch ein Muster erkennbar, da hier die asymmetrische For- 
mung der Platte und dessen Aussparungen eine immer größere Rolle 
spielen. Wie sich im späteren Kapitel zeigt, hat die Form der Stator- 
platte einen erheblichen Einfluss auf die Übertragungsfunktionen der 
Zahnkräfte zur Fahrwerksanbindung. Zu beachten sind die zusammen- 
gehörigen Strukturmodenpaare, die sich teilweise stark in ihrer Reso- 
nanzfrequenz unterscheiden. 


Wirkung der radialen Kräfte auf die Statoreinheit 


Aufgrund der asymmetrischen Geometrie der Statoreinheit (vor allem 
der Statorplatte) entstehen durch Radialkräfte, auch wenn diese sym- 
metrisch über den Statorumfang auf dessen Zähne wirken, resultierende 
Kräfte am Anbindungspunkt des Radnabenantriebs. Dieses Phänomen 
wird nachfolgend genauer erläutert (siehe auch [Str17b]). 


Um das Übertragungsverhalten des Statorsystems etwas gründlicher zu 
betrachten, wird mithilfe der aus der FEM berechneten Massenmatrix 
und der Steifigkeitsmatrix K des Stators ein Gleichungssystem aufge- 
stellt. Damit lässt sich das System folgendermaßen beschreiben. 


Më + Dž + Kx = F. (3.14) 
Mit 


cx = Fp (3.15) 


Die Kraft Fra (abhängig von der betrachteten Kraftrichtung æ am Aus- 
gang) wird dabei an das Fahrwerk über den Anbindungspunkt weiter- 
geleitet, wohingegen F,.; die auf den Stator wirkenden Radialkräfte auf 
den Zahn i beschreibt, wie Abbildung 3.23 verdeutlicht: 
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Zahn i=1 


\ In der Arbeit beschriebene 


A Möglichkeiten für a: 
r: radial zur anregenden Kraft 
N \ x: x-Richtung 
A allgemein: F;, y: y-Richtung 


z: z-Richtung 


Abbildung 3.23: Übersicht der Kraftbeschreibung 


Die Dämpfungsmatrix D wird so gewählt, dass die Systemdämpfung 
dem Standarddämpfungsfaktor für verschraubte Stahlkomponenten von 
0,04 entspricht [Mue00]. Durch Umformung und Einsetzten der Feder- 
steifigkeit erhält man folgende allgemeine Ubertragungsfunktionen Gra 
in Abhängigkeit der betrachteten Kraftrichtungen æ und der aufge- 
brachten Kraft F,;, in radialer Richtung auf den Zahn i. 


Fra (5) a c 


= 3.16 
Fals) Ms?+Ds+K ote) 


Gra(s) = 


Die Anregung beziehungsweise die Krafteinleitung F,; in das System 
erfolgt radial auf die Zahnflächen über einen Punkt der mittig auf der 
Zahnfläche liegt und mit dieser verbunden ist. Zusätzlich ist die Reak- 
tionskraft Fra am Anbindungspunkt in radialer Richtung, also Fp, dar- 
gestellt. 
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3 Quellen und Übertragungswege 


Die radiale Reaktionskraft Fry am Befestigungspunkt des Radnabenan- 
triebes wird in der Übertragungsfunktion G,,.(s) sowie in der Übertra- 
gungsfunktion G,.,(s) in z-Richtung mit den wirkenden Kräften Fpz be- 
schrieben. Betrachtet man das System statisch (0 Hz), erkennt man, 
dass sich die Kräfte auf gegenüberliegenden Zähnen (Zahnpaar) aufhe- 
ben, da die Kraft am Anbindungspunkt gleich der Kraft am Einleitungs- 
punkt ist (Übertragungsfunktion hat den Wert 1 bei 0 Hz - siehe Abbil- 
dung 3.24). Es entstehen somit keine resultierenden Kräfte in x- oder y- 
Richtung. Betrachtet man jedoch höhere Frequenzen, ist dies nicht der 
Fall. Zur Veranschaulichung wurde für jeden Zahn die Übertragungs- 
funktion ermittelt. Abbildung 3.24 zeigt den Verlauf der Übertragungs- 
funktion G,,(s) im Frequenzbereich bis 2000 Hz der radialen Kraft am 
Anbindungspunkt des Stators in Bezug auf radialer Krafteinwirkung am 
Zahn. 


Frr(s) 


Grr(s) = = (s) (3.17) 
_ Fry (s) 
Gra (s) == F. (s) (3 : 1 8) 


Es ist erkennbar, dass die Peaks der Übertragungsfunktionen nicht 
exakt übereinander liegen. Besonders groß ist dieser Unterschied im 
Bereich um 75 Hz und um 430 Hz. Zu erklären ist dies durch die asym- 
metrische Struktur der Statorplatte, was sich bereits in leicht unter- 
schiedlichen Resonanzfrequenzen (M1/M2 und M5/M6) in Abbildung 
3.22 geäußert hat. 
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3.2 Beschreibung der Übertragungswege 
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Abbildung 3.24: Radiale Übertragungswege (rotatorisch) Stator 


Abbildung 3.25 stellt den Verlauf der Übertragungsfunktion in x-, y- 
und z-Richtung dar. Dabei ist zu erkennen, dass auch in z-Richtung re- 
sultierende Wechselkräfte entstehen, welche an das Fahrwerk weiterge- 
geben werden. 


FR 
© 


Verstärkungsfaktor [N/N] 
Verstärkungsfaktor [N/N] 


0 1000 2000 


0 1000 2000 0 1000 2000 


Frequenz [Hz] Frequenz [Hz] Frequenz [Hz] 


Abbildung 3.25: Radiale Ubertragungswege (translatorisch) Stator 


Grund fiir die resultierenden Krafte in z-Richtung ist die Einwirkung 
der radialen Kraft auf den Statorzahn, der einen Hebelarm beziiglich 
des Anbindungspunktes besitzt und somit die angeschlossenen Massen 
zum Schwingen in z-Richtung anregt. Das hierdurch entstehende Bie- 
gemoment regt vor allem die Membranbewegung der Statoreinheit an 
und ruft schließlich im Anbindungspunkt eine Reaktionskraft in z-Rich- 
tung hervor, wie Abbildung 3.26 verdeutlicht. 
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3 Quellen und Übertragungswege 


Zahn i=1 
\ = 


Pi 
7 


Abbildung 3.26: Skizze radial wirkende Kräfte auf Stator 


Um das Übertragungsverhalten des Gesamtsystems darzustellen, be- 
dient man sich nun der berechneten Übertragungsfunktionen der einzel- 


nen Zähne und summiert diese auf, um eine Aussage der resultierenden 
Kräfte am Anbindungspunkt des Radnabenantriebes zu bekommen. 
Aufgrund des Superpositionsprinzips bei linearen Systemen ist dies zu- 
lässig. Die Übertragungsfunktion des Gesamtsystems wird mit Hilfe 


der komplexen Addition berechnet. 


Daraus ergibt sich für die einzelnen Richtungen folgender Zusammen- 


hang: 
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(3.19) 


(3.20) 


(3.21) 


3.2 Beschreibung der Übertragungswege 


= G u G, = G 
Z 0.15 = 2.6 
Z Z Z 
= 0.2 
5 0.1 5 8 0.4 
yg” % 0.15 2» 
g 5 5, 
B0 
= 0.05 gs O4 Boz 
= E 0.05 = 
? 2 20 
200 1000 20002 0 1000 2002 0 1000 2000 
Frequenz [Hz] Frequenz [Hz] Frequenz [Hz] 


Abbildung 3.27: Resultierende Übertragungsfunktionen Stator 


In den resultierenden Übertragungsfunktionen G, .und Gy ist zu erken- 
nen, dass vor allem bei ca. 75Hz, 430Hz und um 1000 Hz Peaks auftre- 
ten, welche zu resultierenden Wechselkräften in x- und y- Richtung füh- 
ren. Verantwortlich dafür ist die Asymmetrie der Statorplatte und des 
Radzapfens. Diese Asymmetrie führt zu Moden, die für die Überhöhun- 
gen für Gy und Gy verantwortlich sind. 


Die Funktion G, zeigt das Entstehen der Wechselkräfte in z-Richtung 
aufgrund der beschriebenen Hebelwirkung. Ausschlaggebend für die 
besonders dominanten Peaks sind die Membranmoden, die in der Ab- 
bildung 3.22 dargestellt sind. 


Ebenso zeigt sich, dass durch eine symmetrische Auslegung der Stator- 
platte sich Reaktionskräfte in x-, und y- Richtung stark verringern las- 
sen [Str17b]. 


Die tangentialen Kräfte hingegen, die auf die Statoreinheit wirken, re- 
gen vor allem die Torsionsmode sehr stark an, was zu einer überhöhten 
Krafteinleitung von Drehmomentschwingungen in das Fahrwerk führt. 
In Kapitel 4 wird im Gesamtsimulationsmodell mit Rotor und Stator 
dieser Übertragungsweg näher erläutert. 
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3 Quellen und Übertragungswege 


3.2.2.2 Rotoreinheit 


Die Rotoreinheit ist über die Radnabe auf dem Radlagerzapfen der Sta- 
toreinheit drehend gelagert. Darauf befestigt ist der Rotorträger und der 
geblechte Rotor, welcher die Permanentmagneten beinhaltet. Frontal 
wird ein Ring zur weiteren Befestigung des Rotorblechs aufgeschraubt. 


Abbildung 3.28: Aufbau Rotoreinheit 


Auch hier wurden vereinfachte Annahmen für die FEM-Simulation ge- 
troffen, welche anschließend stichpunktartig aufgelistet sind. Genauere 
Beschreibungen der Vereinfachungen finden sich in [Paul4]: 


e Einheitsdicke am Ring 

e Kreisrunde Außenfläche am Rotorblech 

e Geometrische Vereinfachungen am Rotorhalter 
o Vereinfachung der Strebengeometrie 
o  Entfall von kleinen Befestigungsbohrungen 
o Entfall von Anfasungen 

e Zusammenfassung von Magnetpaketen 

e Vereinfachung der Schraubenverbindungen 
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3.2 Beschreibung der Übertragungswege 


Modalanalyse 


Ebenso wie bei der Modalanalyse der Statoreinheit wird die Rotorein- 
heit mit einer hohen Federsteifigkeit an die Umgebung gekoppelt. Al- 
lerdings wird hier der Drehfreiheitsgrad um die Rotationsachse nicht 
gekoppelt, da diese Randbedingungen den späteren Randbedingungen 
im Betrieb am nächsten kommen. 


Mode M1 M2 M3 
Frequenz [Hz] 365 365 650 
Beschreibung Kippen x Kippen y Translation z 


M4 M5 M6 M7 
1524 1528 1956 1956 
Quetschen Quentschen Translation y Translation x 


Abbildung 3.29: Modalanalyse Rotoreinheit 


Wirkung der radialen Kräfte auf den Rotor 


Analog zu den radial wirkenden Kräften auf den Stator wird nun das 
Übertragungsverhalten der Radialkräfte auf den Rotor untersucht. Als 
Eingang für die Kräfte dient ein mit den Polflächen verbundener Refe- 
renzpunkt (Abbildung 3.33). 
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3 Quellen und Übertragungswege 
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Abbildung 3.30: Radiale Ubertragungswege (rotatorisch) Rotor 


Es fallt in Abbildung 3.30 auf, dass bei den 40 Ubertragungspfaden des 
Rotors viel kleinere Abweichungen zueinander zu erkennen sind als bei 
der Statoreinheit. Dies liegt an der achssymmetrischen Konstruktion 
des Rotors. Die erkennbar geringen Abweichungen können aufgrund 
kleiner asymmetrischer Abweichungen in der Konstruktion und asym- 
metrischer Vernetzung in der FEM erklärt werden. 


G 


rx Ty 


10° 


Verstarkungsfaktor [N/N] 
E 
Verstarkungsfaktor [N/N] 
Verstärkungsfaktor [N/N] 


0 1000 2000 0 1000 2000 0 1000 2000 
Frequenz [Hz] Frequenz [Hz] Frequenz [Hz] 


Abbildung 3.31: Radiale Übertragungswege (translatorisch) Rotor 


In Abbildung 3.31 ist vor allem der Übertragungsweg in z-Richtung 
auffällig. Zum einen ist erkennbar, dass durch das Aufbringen radial 
wirkender Kräfte auf das angebundene Radlager in z-Richtung entste- 
hen. Ebenso wie beim Stator entstehen diese Kräfte aufgrund der Bie- 
gung des äußeren Rings am Rotor, wie in Abbildung 3.33 angedeutet. 
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3.2 Beschreibung der Übertragungswege 
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Abbildung 3.32: Gesamt resultierende Übertragungsfunktionen des 
Rotors 


Grundsätzlich sind die resultierenden Übertragungsfunktionen im Ver- 
gleich zu den resultierenden Übertragungsfunktionen der Statoreinheit 
weniger stark ausgeprägt. Dass trotzdem Überhöhungen in x- und y- 
Richtung erkennbar sind, kann aufgrund kleiner asymmetrischer Ab- 
weichungen in der Konstruktion und der Vernetzung erklärt werden. 
Die entstehenden Kräfte in z-Richtung sind ähnlich wie beim Stator auf 
die Hebelwirkung der Kräfte und der dabei angeregten Resonanzen zu- 
rückzuführen. 


Abbildung 3.33: Skizze radial wirkender Kräfte auf Rotor 
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3 Quellen und Übertragungswege 


3.2.2.3 Fahrwerk 


Das Fahrwerk des Versuchsträgers besteht aus mehreren Blechen, wel- 
che durch eine Schweißverbindung miteinander verbunden sind. Im 
FEM Modell wurden diese Verbindungen durch kinematische Verbin- 
dungen auf den angrenzenden Flächen realisiert. Abbildung 3.34 ver- 
deutlicht den Aufbau der halben Achskonstruktion. 


Halterung 
für Dämpfer 


Radflansch 


Verstärkungs 
blech 


Längslenkerbuchse Längslenker 


Abbildung 3.34: Aufbau Fahrwerk 


Die bei der Modellierung des Fahrwerkes verwendeten Vereinfachun- 
gen sind nachfolgend stichpunktartig aufgelistet. Es sei aber auf 
[Knö15] verwiesen, wo diese Vereinfachungen und die Betrachtung der 
Übertragungswege ausführlicher dargestellt sind. 


Ersatzsteifigkeit für Längslenkerbuchse 

e Modale Dämpfung für Starrkörperbewegung um die Längslen- 
kerbuchse (LLB) 

e Ersatzsteifigkeit für Aufbaufeder 

Geschwindigkeitsabhängige Dämpferkonstante als Däm-pferer- 

satz 


Nachbildung der Schweißverbindungen durch Klebeverbindun- 
gen 


Mit der Nutzung des Modells in der Gesamtsimulation kann jedoch ge- 
zeigt werden, wie die grundsätzlichen Übertragungswege des Gesamt- 
systems aussehen und welche Optimierungsmöglichkeiten sich daraus 
ergeben. 
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4 


4.1 


Gesamtsimulation 


Allgemeines 


Frühere Veröffentlichungen und Forschungsberichte zeigen viele Mög- 
lichkeiten, ein multiphysikalisches Simulationsmodell in Hinblick auf 
das vibroakustische Verhalten einer elektrischen Maschine aufzubauen, 
in denen die anregenden Kräfte durch die Magnetik und die mechani- 
schen Übertragungspfade abgebildet sind. Dadurch können verschie- 
dene Analyse- oder Optimierungsziele der elektrischen Maschine be- 
reits in der Simulation genauer untersucht werden. Aktuelle Literatur 
und Veröffentlichungen sind nachfolgend aufgezählt: 


[Dup18]: Akustisches Modell zur PWM -Optimierung 

[Put16, Wan18]: Optimierung der Motortopologie auf Grund- 
lage eines Multiphysikalischen Simulationsmodells 

[LeB16]: Schnelle Berechnung des akustischen Verhaltens 
[Sail6a, Sail6b]: Berechnung der Vibrationen am Gehäuse über 
die experimentell gemessenen Transferfunktionen mit Hilfe ei- 
ner Vibrationssynthese 

[Ric15b, Ric15a, Dup12, Tor11]: Schallabstrahlung eines Trak- 
tionsantriebes für elektrisch angetriebene Fahrzeuge 

[Adal4]: Betrachtungen zum Einfluss des Regelkreises auf die 
Schallabstrahlung 

[Dup14]: Einfluss von Fehlern auf die Schallabstrahlung 
[Fan14]: Einfluss der magnetischen Kräfte auf das vibroakusti- 
sche Verhalten eines Antriebstranges mit Getriebe 

[Bral3a]: Analytische Berechnung der Schallabstrahlung basie- 
rend auf elektromagnetischer FEM 

[Bral2]: Multiphysikalisches Modell basierend auf der 
„Lumped Element“-Methode 

{[Hum12, San09]: Optimierung des magnetischen Design mit 
Blick auf das vibroakustische Verhalten des Systems 
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e [Jeal2]: Schwerpunkt auf Homogenisierung von heterogenen 
Materialien in der Struktursimulation 

e  [Pell2a]: Einfluss von zusätzlichen Harmonischen in den Sta- 
torströmen auf das Gesamtsystem 

e [Torl2]: Einfluss von Exzentrizität und Drehmomentwelligkeit 
auf das vibroakustische Verhalten mit Hilfe der Feldrekonstruk- 
tionsmethode 

e [Jea08]: Mehr-Kriterien-Optimierung einer elektrischen Ma- 
schine in Bezug auf das vibroakustische Verhalten 

e [Sch07]: Einfluss des Gehäusemateriales auf das vibroakusti- 
sche Verhalten 


In dieser Arbeit wird das in Kapitel 3.1.3 beschriebene elektromagneti- 
sche Modell (Kraftrekonstruktion) mit den aus der mechanischen FEM 
gewonnen Übertragungsfunktionen eines Viertelfahrzeugmodells ge- 
koppelt. Damit können Effekte von der Entstehung störender Schwin- 
gungen bis zur Übertragung der Kräfte zu den Anbindungspunkten gut 
erklärt und schnell simuliert werden. Somit bietet sich so die Möglich- 
keit schnell und einfach neue Optimierungsansätze in der Ansteuerung 
zu testen, vor allem im Hinblick auf einen akustisch optimierten Steue- 
rungsansatz. 


Vorteile des in dieser Forschungsarbeit entwickelten multiphysikali- 
schen Simulationsmodells sind: 


e Übersichtlichkeit des Modells 

o Effekte können sehr gut und einfach dargestellt werden. 
Dies dient vor allem der Nachvollziehbarkeit der Entste- 
hungsmechanismen. 

o Einfache Änderungsmöglichkeiten bei geänderten Rand- 
bedingungen (z.B. Bestromung, Anbindungssteifigkeiten, 
etc.) 

e Sehr schnelle Berechnung (Variationsmöglichkeiten) 
e Riickrechenbarkeit für optimale Bestromung (Alleinstel- 
lungsmerkmal) 
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4.1 Allgemeines 


Das verwendete Modell gliedert sich auf in ein elektromagnetisches 
Modell mit den Strömen I, und I4 als Eingänge und ein mechanisches 
Modell bei dem die wirkenden Kräfte im Luftspalt als Eingangsgröße 
dienen. Abbildung 4.1 zeigt den Aufbau des Simulationsmodells. 
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Abbildung 4.1: Aufbau Gesamtsimulation 
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4.2 Aufbau der Gesamtsimulation 


4.2.1 Elektromagnetische Simulation 


Für das Gesamtmodell wurden die Ursachen für die Kraftschwankun- 
gen im elektromagnetischen Modell getrennt (siehe Kapitel 3.1.1.1 und 
3.1.1.2). Die konstanten Anteile der Ströme werden durch den zu tes- 
tenden Betriebspunkt (Drehmoment und Drehzahl) bestimmt. Die da- 
zugehörige Ansteuertabelle für die /,- und /4-Ströme findet sich in Ab- 


bildung 4.2 wieder. 
400 
= = 
z 7 2 -100 
2 Bd 
= = 200 
-400 
1000 1000 1000 1000 
Drehzahl [rpm] Oe Drehzahl un 
renza] rpm renza! rpm 
p 0 -1000 Moment [Nm] [pm] 0 -1000 Moment [Nm] 


Abbildung 4.2: Ansteuertabellen für J, und Ia 


Dabei entstehende Oberschwingungen in den Strömen können mit Hilfe 
der in Kapitel 3.1.3 entstandenen Werte (Abbildung 3.11) nachgebildet 
werden. Nicht berechnete Betriebspunkte zwischen den berechneten 
Stützstellen werden linear interpoliert. Die bisherigen Untersuchungen 
(Kapitel 3.1.3.3) haben gezeigt, dass für die harmonischen Anteile in 
den Ersatzströmen ein lineares Übertragungsverhalten angenommen 
werden kann. 
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4.2 Aufbau der Gesamtsimulation 


4.2.2 Mechanische Simulation 


Aufgrund der Übertragung der störenden Phänomene durch die Struktur 
im mechanischen System und der dadurch verursachten kleinen Bewe- 
gung werden hier linearisierte Zusammenhänge angenommen. Grund- 
sätzlich werden alle genutzten Übertragungsfunktionen aus einem 
FEM-Gesamtmodell mit gekoppelten Komponenten ermittelt. Die Wir- 
kung der Kräfte auf die Statorzähne wird mit Hilfe der aus der FEM 
gewonnen Übertragungsfunktionen berechnet. Ebenso gilt dies für die 
auf das Radlager wirkenden Kräfte, deren Ursache vor allem auf den 
Rotor angreifende Kräfte sind. Wie in Kapitel 3.2.1.1 beschrieben, wer- 
den die Frequenzen der mitrotierenden Kräfte in x- und y-Richtung ent- 
sprechend aufgeteilt. Somit ergeben sich daraus die Nebenordnungen 
mit v + 1 und v — 1. Gleiches gilt für die angreifenden Drehmomente 
am Radlager um x und y. Wie auch später im Experiment wurden nur 
die Komponenten Statoreinheit, Rotoreinheit und Fahrwerk untersucht. 


Anbindung 
Dämpfer 


Anbindung 
Aufbaufeder 


Anbindung 
Längslenkerbuchse (LLB) 


Abbildung 4.3: Mechanisches Modell 
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4.3 Simulationsergebnisse 


Im folgenden Kapitel wird eine Übersicht der im System des Radna- 
benantriebs möglichen Eigenschwingungen und Übertragungsfunktio- 
nen dargestellt und diese nach deren Einfluss auf das vibroakustische 
Verhalten bewertet. 


4.3.1 Modalanalyse 


Die sich aus dem Gesamtsystem ergebenden Eigenfrequenzen sind teil- 
weise anhängig von den Randbedingungen. In diesem Kapitel werden 
ausgehend von dem gezeigten Gesamtmodell ohne Ersatzmodelle für 
Reifen und Felge Schritt für Schritt zusätzliche Randbedingungen hin- 
zugefügt und deren Schwingformen sowie Übertragungswege darge- 
stellt. Hierfür wurden die in Abbildung 4.4 beschriebenen Randbedin- 
gungen gewählt. 


Feder Creaer = 20 N/mm 


Punktmasse: 
m = 0,02t 


600 0 0 
I=| 0 60 0 en? 
0 


0 0 95 


LLBY Reifensteifigkeit: 


Cr11 = 300 N/mm 
Cr22 = 300 N/mm 
Cr33 = 300 N/mm 
Crag = 19- 106 Nr 


LLBx 


Ciipx = 1064 N nm 


Cıupy = 331 N 
LLBy /mm Crss = 19- 106 Be 


Abbildung 4.4: Gesamtmodell mit Darstellung der Randbedingungen 
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4.3 Simulationsergebnisse 


Die folgende Tabelle zeigt die Unterschiede der variierten Randbedin- 


gungen: 


Nr 


Randbedingung und 
Erweiterung der RB 


Ohne Komplettrad 

nur: 

Steifigkeit Längslenkerbuchse 
[Kn615] 

Steifigkeit Aufbaufeder ohne Dämp- 
fer 


Mit Komplettrad abgehoben 
zusätzlich: 

Punktmasse als Ersatzmasse für 
Komplettrad 


Mit Komplettrad auf Fahrbahn 
zusätzlich: 

Steifigkeit zwischen Punktmasse 
mit Umgebung zur Nachbildung der 
Reifensteifigkeit 

(wie Abbildung 4.4) 


Darstellung 


Tabelle 4.1: Überblick der geänderten Randbedingungen 


Tabelle 4.2 zeigt die Abhängig- oder Unabhängigkeiten der Eigenfre- 
quenzen von den verschiedenen Randbedingungen. Es ist deutlich zu 
erkennen, dass vor allem die Eigenfrequenzen bei niedrigen Frequen- 
zen sehr stark abhängig sind. 
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Nr 


Frequenz 


42: 47 
35; 53 
39; 55 


47 
45 
45 


65 
81 
82 


110 
68 
69 


214 
152 
152 


443 
377 
377 


Bezeichnung 


Rotation Rotor 


Kippen gesamt 


Kippen gesamt 


Rotation Stator mit 
Fahrwerk 


Kippen Stator 


Kippen Stator 


Membran Statorplatte 
mit Translation Stator 


Translation Rotor 1 


55666558, 


4.3 Simulationsergebnisse 


2 412 Translation Rotor 2 


520 
520 Oval 1 
520 


521 

521 Oval 2 
521 
>1000 
665 
665 
>1000 
680 Kippen Rotor 2 
680 


WN -e 


Kippen Rotor 1 


WU N-e U Ne UN 


1 903-1119 
2 862-1107 
3 862-1108 


Membran Moden der 
Statorplatte 


psbb 


Tabelle 4.2: Übersicht der Eigenfrequenzen und deren Abhängigkeiten 
im Gesamtsystem mit den in Tabelle 4.1 dargestellten unterschiedli- 
chen Randbedingungen 


Analog zu der in Kapitel 4.2.2 gezeigten Möglichkeit zur Berechnung 
und Darstellung von resultierenden Kräften zeigt Abbildung 4.5 diese 
Übertragungsfunktionen. Es sind die von den Rotorpolen radial wirken- 
den Kräfte zu den radial wirkenden Kräften am Radlager dargestellt, 
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4 Gesamtsimulation 


welche an den Radzapfen weitergegeben werden. Dieser Übertragungs- 
weg wurde bereits in Kapitel 3.2.1.1 genauer beschreiben. 
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Abbildung 4.5: Darstellung der radialen Ubertragungswege des Rotors 
auf das Radlager im Gesamtsystem und deren Einfluss auf der in Ta- 
belle 4.1 dargestellten Randbedingungen 


Es zeigen sich Abhängigkeiten von den Randbedingungen. Besonders 
deutlich wird dies, wenn dazu die Übertragungsfunktion zu den jewei- 
ligen Anbindungspunkten an der Karosserie im Gesamtsystem gebildet 
wird (Abbildung 4.7). 
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4.3.2  Übertragungswege des Gesamtsystems 


Ebenso können die Übertragungswege, von den wirkenden Kräften auf 
Rotor und Stator hin zu den Anbindungspunkten des Fahrwerksgebildet 
werden, wie Abbildung 4.6 zeigen soll. 


FDämpfer 


Freder 


Fitpy 


Fitpx 


Abbildung 4.6: Veranschaulichung der zu den Ubertragungsfunktio- 
nen gehörenden Anregungs- und Anbindungskräfte 


So ergeben sich fiir die radiale und tangentiale Anregungsrichtungen 
folgende Ubertragungsfunktionen für jeden Zahn beziehungsweise Ro- 
torpol i: 


Gi tps (s) 
Fi papper (S itDampfer 
Girdampfer(s) = —— _ Fipämpfer (8) (4.1) 
F(s) 

_ Fireaer(S) _ Fi Feder (s) 
Girreaer (S) = Erw Gitreaer(S) = El (4.2) 

= Fi pampfer (S) _ Fittpx(S) 
Girtipx(s) = Eee Gittiex(S) = Fs) (4.3) 
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Fi tay (S) FiLLBy (s) 
Gir =z Gi = — 4.4 
ruby (S) FF) itLLBey(S) F.,(5) (4.4) 


Die resultierenden Kräfte (am Beispiel der Statorkräfte zu den Anbin- 
dungspunkten mit 48 Zähnen) können somit folgendermaßen berechnet 


werden: 
48 Gepamp fer 
48 
Gypamp fer = a _ G (4.5) 
i=1 = i,tDämpfer 
i=1 
48 48 
GrFeder = Ni Gtreder = > Gi treaer (4.6) 
et i=1 
48 48 
Gruusx = >, e Gtirex = > Gi tLLBx (4.7) 
tt i=1 
48 48 
Grusy = >, Ga GtLLBy = > Gi tLLBy (4.8) 
i=1 £ 
i=1 


Betrachtet man das Gesamtübertragungsverhalten der an den Anbin- 
dungspunkten des Fahrwerks auf die Karosserie entstehenden resultie- 
renden Kräfte aufgrund radialer Anregungskräfte am Rotor in Abbil- 
dung 4.7, ist zu erkennen, dass das Einbringen der Reifensteifigkeit nur 
noch wenig Einfluss auf die vom Radlager übertragenen Kräfte hat. 
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Resultierende Resultierende 
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Abbildung 4.7: Resultierende Übertragungsfunktionen bei Variation 
der Randbedingungen 
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Die nachfolgend dargestellten Übertragungsfunktionen (Abbildung 
4.8) zeigen die Übertragungswege der tangential anregenden Kräfte für 
Stator und Rotor zu den Anbindungspunkten des Fahrwerks an die Ka- 
rosserie, wie sie bereits in Kapitel 3.2.1.2 und 3.2.1.4 beschrieben wur- 
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Abbildung 4.8: Vergleich der Ubertragungswege mit tangentialen An- 
regungskräften 
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Außerdem können die radialen Übertragungswege (Resultierende 
Übertragungsfunktion) aus Kapitel 3.2.1.1 (blau) und 3.2.1.3 (rot) mit- 
einander unter der Randbedingung (3) verglichen werden, wie die nach- 
folgende Abbildung 4.9 zeigt. 
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Abbildung 4.9: Vergleich der Ubertragungswege mit radialen Anre- 
gungskraften 
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Im Vergleich zu den radialen Übertragungswegen sind die tangentialen 
Übertragungswege sehr viel stärker ausgeprägt. Aufgrund der hohen ra- 
dialen Kräfte, die auf Rotor und Stator wirken, besteht jedoch trotzdem 
die Möglichkeit, dass diese im Innenraum spür- und hörbare Effekte 
hervorrufen. Außerdem ist erkennbar, dass die tangentialen Übertra- 
gungswege beim Stator wesentlich stärker ausgeprägt sind als beim Ro- 
tor. Folgende Schlussfolgerungen lassen sich aus den gezeigten Über- 
tragungswegen ableiten: 


1. Die tangentialen Kräfte auf den Stator sind ausschlaggebend 
für die vibroakustischen Eigenschaften, die an den Anbin- 
dungspunkten und schließlich in der Fahrzeugkabine gemes- 
sen werden können. 

2. Die tangentialen Kräfte auf den Rotor spielen eine untergeord- 
nete Rolle. 

3. Radiale Kräfte, die auf Rotor und Stator wirken, können Ein- 
fluss auf das vibroakustische Verhalten haben. 


Abhängigkeiten der Übertragungsfunktionen von den Randbedin- 
gungen 


Ebenso werden hier die entstehenden Unterschiede der Übertragungs- 
funktionen für die verschiedenen Randbedingungen aufgezeigt. Der am 
stärksten ausgeprägte Übertragungspfad (Tangentialkräfte zum Aus- 
gangspunkt) wird aus diesem Grund in Abbildung 4.10 dargestellt. 
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Abbildung 4.10: Abhangigkeit des tangentialen Ubertragungsweges 
des Stators von den Randbedingungen 


Es ist zu erkennen, dass das Ubertragungsverhalten vor allem bei nied- 
rigen Frequenzen abhängig von den Werten der Randbedingungen ist, 
welche durchaus bei verschiedenen Betriebspunkten variieren können. 
Gründe hierfür sind: 


e Reifensteifigkeit/Reifenmassen 

o  Reifentemperatur 

o Wechsel der Reifen 

o Alterung und Abnutzung der Reifen 
e Dämpferkonstante 

o  Frequenzabhiangigkeiten 

o Alterung 

o Temperatur 
e Steifigkeit Längslenkerbuchse 

o Alterung /Temperatur 
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4.3.3 Motorordnungen 


Mit Hilfe der vorher dargestellten Effekte und den Modellen kann für 
das System des Radnabenantriebs Tabelle 4.3 erstellt werden, welche 
die Hauptursachen der zeitlichen Ordnungen der wirkenden Kräfte am 
Anbindungspunkt näher beschreibt. Zusätzlich treten zu den beschrie- 
benen Hauptordnungen auch harmonische Ordnungen auf, indem diese 
mit į multipliziert werden. Dabei ist i positiv und ganzzahlig. Die Höhe 
der wirkenden Amplituden hingegen kann nur durch elektromagneti- 
sche Berechnungen genau ermittelt werden, wie sie beispielsweise die 


FEM liefert. 
Allgemeine Ordnung 
Berechnun E-Wheel Drive 
me beta 
en 40-i 
N'i 48-i 
NFRA 48-it1 


kgV(w, 2p) -i kgV (24,40) «i 


kgV(N,2p)-i  kgV (48,40) -i 


Hauptsächliche 
Ursprungskraft 


radiale Kräfte auf Stator 

> asym. Übertragungswege 
des Stators 

> resultierende Axialkräfte 


radiale Kräfte auf Rotor 
> resultierende Axialkräfte 


radiale Kräfte auf Rotor 
> asym. ÜW des Rotors 
+ Modulierung durch Rotor 


tangentiale Kräfte 
(Lastpulsation) 

> Übertragungspfad Stator 
stark ausgeprägt 


tangentiale Kräfte 
(Nutrastmoment) 

> Übertragungspfad Stator 
stark ausgeprägt 


Tabelle 4.3: Zuordnung der entstehenden Ordnungen 
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Tabelle 4.3 zeigt die entstehenden Ordnungen am Anbindungspunkt der 
elektrischen Maschine bei symmetrischer elektromagnetischer Anre- 
gung für die verwendeten Wicklungsanordnungen. Dabei unterschei- 
den sich Lastpulsation und Nutrastmoment nur in den zugrunde liegen- 
den kleinsten gemeinsamen Vielfachen zwischen Polpaarzahl 2p, Nut- 
zahl N und der Wiederholung der Wicklungsanordnung w. Das 
Nutrastmoment kann dabei direkt aus der Geometrie der elektrischen 
Maschine abgelesen werden, wohingegen bei der Lastpulsation auch 
das Wicklungsschema eine entscheidende Rolle spielt. Entscheidend ist 
hier die Wiederholung der Wicklungsanordnung um den Umfang (siehe 
Abbildung 4.11). 


Beispiel: Wicklungen E-Wheel Drive “beta” 


4 Wicklungs- 
8tel Stator anordnung 


J 
aa“ Eae ere 
pe as? 
D . Ta 


eo 
y or z 
Se A N 


S N N ® A 
Z D Wicklung U 
z Wicklung V 
Wicklung W 
>w = 24 


Beispiel: Wicklungen E-Wheel Drive “gamma” 
20tel Stator 


icklungs na a 
an” znoranund = 
AC! 
an nnd 
K 


ane 


Wicklung U 
Wicklung V 
Wicklung W 


=>w = 60 


Abbildung 4.11: Wiederholung der Wicklungsanordnung 
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In realen Systemen kommt es zusätzlich zu Abweichungen dieser sym- 
metrischen Anregung zum Beispiel aufgrund von Exzentrizitäten. 
Dadurch sind die gegenüberliegenden Kräfte, welche auf Rotorpole und 
Statorzähne wirken, nicht mehr gleich [Kes16, Kral7]. Es ändern sich 


die Amplituden der Hauptanregungen. Dies zeigt sich in Abbildung 
4.12. 
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Abbildung 4.12: Darstellung der Kräfte über den Winkel auf einem 
Zahn bei geändertem Luftspalt, statische Exzentrizität, Betriebspunkt 
100Nm/100rpm 


Vor allem die Amplitude der 40. Ordnung in den Radialkräften ist stark 
abhängig von der Luftspaltgröße. Ähnliches gilt auch für die 48. Ord- 
nung der Radialkräfte, die auf die Rotorpole wirken. Zusätzlich entste- 
hen bei dynamischen Exzentrizitäten Nebenordnungen beziehungs- 
weise Seitenbänder in den anregenden Kräften, deren Amplituden sich 
je nach betrachteten Zahn oder Pol stark unterscheiden können [Mar13, 
Piel2]. Im Folgenden wird allerdings zur Vereinfachung von einem 
symmetrischen System ausgegangen. 
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4.3.4 Wirkende Kräfte 


Mit den bereits beschriebenen Mechanismen werden nun die wirkenden 
Krafte an den Anbindungspunkten des Fahrwerks zur Karosserie nach- 
gebildet. Dabei wird, wie zuvor beschrieben, von einem elektromagne- 
tisch symmetrischen Modell ausgegangen. 


Für die Berechnung beziehungsweise Nachbildung der verschiedenen 
Ordnungen müssen folgende Fälle unterschieden werden, wie bereits 
Tabelle 4.3 zeigte: 


1. Fall: Vielfaches der Polpaarzahl 
=> erregt durch radial wirkende Kräfte auf Stator mit dessen 
relevanten Übertragungswegen 
Beispiel am Radnabenantrieb (beta): 
40. Ordnung, 80. Ordnung, 120. Ordnung, ... 


2. Fall: Vielfaches der Statornutzahl 
=> erregt durch radial wirkende Kräfte auf Rotor mit dessen 
relevanten Übertragungswegen in axiale Richtung 
Beispiel am Radnabenantrieb (beta): 
48. Ordnung, 96. Ordnung, ... 


3. Fall: Vielfaches der Statornutzahl +/- 1 
=> erregt durch radial wirkende Kräfte auf Rotor und den da- 
raus entstehende rotierenden resultierenden in nicht axialer 
Richtung 
Beispiel am Radnabenantrieb (beta): 
47. Ordnung, 49. Ordnung, 95. Ordnung, ... 


4. Fall: Kleinstes gemeinsames Vielfaches der Polpaarzahl und der 
Wiederholung der Wicklungsanordnung 
=> erregt durch tangentiale Kräfte auf Rotor und Stator 
Beispiel am Radnabenantrieb: 


120. Ordnung, 240. Ordnung, ... 
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Bei der Beschreibung der Übertragungspfade soll Abbildung 4.13 un- 
terstützten. Dabei sind die Übertragungsfunktionen nun als Verschie- 
bung zu Kraft definiert. Es sind nur die translatorischen Raumrichtun- 
gen und Kräfte für die Anbindung des Radlagers aufgezeigt, um die 
angewandte Methodik vereinfacht erklären zu können. Für die späteren 
gezeigten Berechnungen wurden allerdings ebenso die dazugehörigen 
rotatorischen Effekte beziehungsweise die auftretenden Drehmomente 
mitberücksichtigt. 


m> Ubertragungsart (m = mechanisch) 


mst? Bauteil und Richtung (St = Stator tangential, Rr = Rotor radial, L = Radlager ) 


Z LF — Ausgangsort (F = Anbindungspunkt Fahrwerk; L = Anbindungspunkt Radlager) 


> Eingang der aufgebrachten Kraft (Zi = i-ter Zahn; Pj = j-ter Rotorpol; Lx = Radlager x-Richtung) 


Eingang Ausgang Beschreibung 


i=Zahn 
F = Ausgang am Fahrwerk 


msr mst 
Gz i_F Gzi F 


j=Pol 
L = Raumrichtung Radlager 


mRr 
Gp J-E 


Lz = axiale Raumrichtung Radlager 
L = alle nicht axialen Raumrichtungen 
F = Ausgang am Fahrwerk 


mLz mL 
Giz F GEF 


j = Pol 
F = Ausgang am Fahrwerk 


mRt 
Gpj F 


Abbildung 4.13: Nomenklatur und Veranschaulichung der Übertra- 
gungspfade 
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4.3 Simulationsergebnisse 


1. Fall: Schwingungen stammen von Radialkräften am Stator. 


Die am Anbindungspunkt des Fahrwerks wirkenden Kräfte zur Karos- 
serie Fr für die beschriebenen Ordnungen können aus der Summe der 
anregenden Kräfte in radialer Richtung F7; auf den Zahn i mit den da- 
zugehörigen mechanischen Übertragungswegen G ar beschrieben 

werden. Zur Berechnung der wirkenden Kräfte wird die gerade be- 
schriebene Auslenkung mit der dort vorliegenden Steifigkeit cp multi- 
pliziert. Die anregenden Kräfte sind zudem von Motormoment M und 


Drehzahl n abhängig. 
48 
Fr Fan = y F7,(M,n) > GH Cr (4.9) 
i=1 


2. Fall: Schwingungen stammen von den resultierenden Axialkräften 
des Rotors (Kapitel 3.2.2.2), welche über das Radlager übertragen wer- 
den und von Radialkräften am Rotor stammen. 


Die vom Radlager übertragene Kraft in axiale Richtung F,, bildet sich 
aus der Summe der auf die Pole P radial wirkenden Kräfte Fp; mit des- 
sen axial wirkenden Übertragungswegen G Pilz und deren Multiplika- 


tion mit der Radlagersteifigkeit c,. 


Kraft auf Radlager durch Radialkräfte Rotor (nur für Axialkräfte) 


„=> Fane Gis ta” ti (4.10) 


Somit kann die wirkende Kraft Fp Fauz am Anbindungspunkt des Fahr- 
werkes durch die axial wirkende Radlagerkraft F,,, multipliziert mit der 
Übertragungsfunktion Gre. und der am Anbindungspunkt vorliegen- 
den Steifigkeit cr berechnet werden. 
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Axialkraft aus Radlager zu Kraft auf die Anbindungspunkte des Fahr- 
werks zur Karosserie 


Fr Fauz = Fiz’ Girt ‘Cr (4.11) 


3. Fall: Schwingungen stammen aus rotierenden Kraften am Radlager 


Die vom Radlager übertragene rotierende Kraft F; rot bildet sich aus 

den resultierenden Kraften (Kapitel 3.2.2.2), welche durch die auf die 

Pole radial wirkenden Kräfte Fp; mit deren nicht axial wirkenden Uber- 
mRr 


tragungswegen Gp; , entstehen. Durch Multiplikation mit der Steifig- 
keit des Radlagers c,, kann die Kraft berechnet werden 


Kraft auf Radlager durch Radialkräfte Rotor (für alle Richtungen außer 
axial) 


40 
Fi rot = = Fp,(M,n) i GPL “CL (4.12) 
j=1 


Da in diesem Schritt rotierende Kräfte beschrieben sind und diese sich 
auf ein stehendendes Koordinatensystem übertragen, wird die übertra- 
gene Schwingung in ihrer Ordnung moduliert (siehe Kapitel 3.2.1.1). 
So kann der Zusammenhang zwischen der rotierenden Kraft F; rot und 
der translatorischen Kraft F; trans folgendermaßen beschrieben werden 


Brana 
L.trang "77 Rot sin((v — 1)wt) + sin((v + Dot 


Kraft an den Anbindungspunkten aus Radlager (für alle Richtungen au- 
er axial) 


Fr raus = Fi trans ' GIF "CL (4.14) 
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4. Fall: Schwingungen aus Tangentialkräften 


Hinzu kommen die eingeleiteten Schwingungen aus den tangential wir- 
kenden Kräften von Rotor Ff j und Stator Fé, mit deren Ubertragungs- 


wegen Gp} $ und Gg , die folgendermaßen berechnet werden können: 


40 48 
Fr rans = ) FE CHE + ) FE -GES ce (4.15) 
j=l i=1 


Mit Hilfe der eben dargestellten Modellierungsvorschriften können die 
wirkenden Krafte an den Anbindungspunkten des Fahrwerks zur Ka- 
rosserie schnell und effektiv berechnet werden. Abbildung 4.13 zeigt 
beispielhaft die Berechnung der wirkenden Kräfte auf die Langslenker- 
buchsen in x- und z-Richtung. 
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Abbildung 4.14: Wirkende Krafte auf Anbindungspunkte 


Mit Hilfe dieses Modells können verschiedene Erkenntnisse gewonnen 


werden: 


Anteile der wirkenden Kräfte aus Tangential- und Radialkräften 


e Im symmetrischen Modell sind die Anteile der durch rotierende 
Radialkräfte auf den Rotor hervorgerufenen Effekte gering (47., 


49.) 


e Auch die daraus entstehenden Axialkräfte, die durch die Geo- 
metrie des Rotors entstehen, sind in den Auswirkungen auf die 
Kräfte der Längslenkerbuchse sehr klein (48.) 
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4.3 Simulationsergebnisse 


Verhalten der wirkenden Kräfte bei geänderter Bestromung und 
geänderten Betriebspunkten 


Abbildung 4.15 zeigt das Verhalten der 120. Ordnung bei Änderung des 
Betriebspunktes. Es ist zu erkennen, dass besonders in den Resonanz- 
stellen die auftretenden Effekte besonders stark hervortreten. Aller- 
dings hat auch die Bestromung bzw. das geforderte Drehmoment Aus- 
wirkungen auf die am Anbindungspunkt wirkenden Kräfte. So ist bei 
300 Nm ein Minimum zu erkennen, was sich, wie Abbildung 3.11 zeigt, 
auf eine hier schwächere Anregung der Tangential- sowie Radialkräfte 
zurückführen lässt. 


Wirkung auf LLBx der 120. Ordnung 


600 


400 
200 200 


Drehzahl [rpm] 300 0 


100 


Drehmoment [Nm] 


Abbildung 4.15: Wirkende Kräfte bei geänderten Betriebspunkten 
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4 Gesamtsimulation 


4.4 Validierung des Modells 


Grundsätzlich werden für die Validierung der Simulation verschiedene 
Vergleiche herangezogen. Zum einen werden die gemessenen Eigen- 
frequenzen des Systems mit den berechneten verglichen, um vor allem 
Informationen über die Statoreinheit des Systems genauer zu erhalten. 
Für eine Aussage der Gültigkeit des Gesamtmodells wird zudem das 
Übertragungsverhalten aus Simulation und Messung gegenübergestellt. 


4.4.1 Eigenfrequenzen des mechanischen 
Systems 


Das Simulationsmodell des mechanischen Systems wurde anhand einer 
Betriebsschwingungsanalyse durch Gegenüberstellung der Lage der Ei- 
genfrequenzen mit dem Modell verglichen. Dabei war der Motor starr 
am Prüfstandsadapter angeschraubt und mit triaxialen Beschleuni- 
gungssensoren versehen (Abbildung 4.16). 


Aus Betriebsschwingungsanalyse (Versuch): 


Kippen 1 Kippen2 Membran Torsion Oval Verschiebung 
65 Hz 79 Hz 131 Hz 314 Hz 388 Hz 458 Hz 
“eB I~ AN hy AN ER 
RY t K y} K ye fe VY K ER 

WSA Ara LA N AL "a 

Aus FEM (nur Stator): 

70 Hz 84 Hz 135 Hz 175 Hz 524 Hz 437 Hz 


Abbildung 4.16: Vergleich der Eigenfrequenzen 


Grundsätzlich ist beim Vergleich der Eigenfrequenzen zwischen Simu- 
lation und experimenteller Betriebschwingungsanalyse eine ausrei- 
chend gute Übereinstimmung zu erkennen, sodass daraus geschlossen 
werden kann, dass die wesentlichen Effekte auch mit dem in der Arbeit 
entwickelten mechanischen Modell erklärt und nachgebildet werden 
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4.4 Validierung des Modells 


können. Eine amplituden- und frequenzgetreue Nachbildung aller im 
Antrieb vorhanden Übertragungswege ist mit diesem Modell allerdings 
nicht möglich. Dies liegt vor allem an der Modellierung der Verguss- 
masse im Stator mit deren komplexen Dämpfungseigenschaften 
[Mil15]. Außerdem unterscheiden sich zwischen Simulation und Expe- 
riment die Torsionsmoden. Dies liegt vor allem an der modellierten An- 
bindungsart, in der lediglich ein kleiner Radius um die Befestigungs- 
schrauben des Radnabenantriebs an die feste Umgebung in der Simula- 
tion gekoppelt wurde. In der Realität ist die Verbindung zwischen 
Radnabenantrieb und Befestigungswinkel deutlich steifer. 


4.4.2 Eingebrachte Störungen 


Ebenso können die aufgebrachten Kräfte mit Hilfe eines vereinfachten 
Aufbaus gemessen und mit dem Modell verglichen werden. 


Kraft in x-Richtung [N] Kraft in y-Richtung[N] Kraftin z-Richtung [N] 


Drehzahl [rpm] Drehzahl [rpm] Drehzahl [rpm] 


Frequenz [Hz] 


Abbildung 4.17: Versuchsaufbau und Festlegung des Koordinatensys- 
tems 
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4 Gesamtsimulation 


Eine Darstellung und Vergleich der translatorischen Kräfte auf die 
Kraftmessdose, welche unter der Statorplatte sitzt, zwischen Simulation 
und Experiment zeigt nachfolgende Darstellung. 


Kraft in x-Richtung (Messung) 
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Abbildung 4.18: Vergleich der wirkenden Kräfte zwischen Simulation 
und Messung 
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4.4 Validierung des Modells 


In Abbildung 4.18 sind Gemeinsamkeiten erkennbar, wie beispiels- 
weise der Übertragungspfad der 240. Ordnung. Teilweise treten aber 
auch stärkere Abweichungen auf. Dies liegt unter anderem an folgen- 
den Gründen: 


e Strom mit Oberschwingungen in Messung und ohne Ober- 
schwingungen in Simulation 

e Annahme der Symmetrie in elektromagnetischem Design (Er- 
klärung der stärkeren Abweichung der 47. Ordnung) im Zusam- 
menspiel einer einseitigen Lagerung 

e Vereinfachung der Geometrie, Verbindungselemente und Mate- 
rialeigenschaften des mechanischen Modells 


Trotzdem können mit der Simulation die grundsätzlichen Mechanismen 
der störenden Kräfte gezeigt werden. 


4.4.3 Ubertragungsfunktionen 


Die Simulationskette für harmonische Schwingungen in den Strömen 
kann effektiv mit dem Experiment verglichen werden, indem man die 
gemessenen und simulierten Übertragungswege gegenüberstellt. Um 
das Übertragungsverhalten im Experiment zu bestimmen, werden 
Schwingungen auf die Ersatzströme aufgeprägt und mit Hilfe der Mes- 
sung an den Ausgängen deren Übertragungsverhalten berechnet. Um 
eine höhere Genauigkeit zu erzielen, wird diese Messung mehrmals mit 
unterschiedlichen Amplituden und Phasenlagen durchgeführt und die 
Punkte für die Übertragungsfunktion mit Hilfe der kleinsten Quadrate 
Methode ermittelt. 


Als Versuchsträger für den Radnabenantrieb wurde ein am Institut für 
Fahrzeugsystemtechnik vibroakustisch optimierter Prüfstandsadapter 
für Fahrwerke entwickelt und aufgebaut, welcher am Akustik-Allrad- 
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4 Gesamtsimulation 


Rollenprüfstand verwendet wurde. An diesem wurde, wie im Versuchs- 
träger (Ford Fiesta), eine Koppellenkerachse mit Radnabenantrieb ver- 
baut. Um auf die Ströme die gewünschten Oberschwingungen aufzu- 
prägen, wurde eine flexible Ansteuereinheit genutzt (siehe Anhang). 


Messpunkt 
Dämpferkraft 


Messpunkt 
Beschleunigung 
LLB i 


Ansteuereinheit 5 
zur Aufpragung 
= von 
Schwingungen 
in die Ströme 


DD N 


Abbildung 4.19: Aufbau für Messung der Übertragungswege 


Mit dem Messaufbau in Abbildung 4.19 werden die ordnungsabhängi- 
gen Parameter eines linearen Modells bestimmt. Dazu wird bei unter- 
schiedlichen Drehzahlen der in Abbildung 4.20 dargestellte Messablauf 
mit und ohne aufgeprägte Schwingungen durchgeführt, indem auf die 
Sollstromvorgabe die jeweilig zu untersuchende Ordnung aufgeprägt 
wird. Die Messungen benötigen für jeden Berechnungspunkt vier me- 
chanische Umdrehungen 
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4.4 Validierung des Modells 


Messablauf für eine aufgeprägte Ordnung bei einer Drehzahl 
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Abbildung 4.20: Messablauf für Parametersuche 


Aus den Messungen wird der Verstärkungsfaktor a für die betrachtete 
Ordnung bestimmt. Hierzu werden aus den Fenstern (jeweils eine me- 
chanische Umdrehung) der Messignale die Fouriertransformierten für 
Y (Yre ohne aufgeprägte Schwingungen in den Strömen; Yauf = mit 
aufgeprägten Schwingungen) und analog dazu die Fouriertransformier- 
ten der Ströme Irer und luf für die zu untersuchende Ordnung be- 
stimmt. Durch die Differenzbildung der Messung mit aufgeprägten 
Schwingungen mit der Referenzmessung werden sowohl Schwingun- 
gen, die ohne aufgeprägte Ströme vorhanden sind, als auch unbeabsich- 
tigt auftretende Abweichungen im Stromregelverfahren mit berücksich- 
tigt. 


= lYaut — YRef | 


= (4.16) 
Haut = Irel 


Ebenso lässt sich die Phasenverschiebung y aus der Differenz bestim- 
men: 


y= arg(Yauf z Yret) = arg(/aur aa Tret) (4.17) 
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4 Gesamtsimulation 


Abbildung 4.21 verdeutlicht diese Vorgehensweise und den dazugehö- 
rigen Signalfluss. 


Sollströme 
(mit/ohne) 
aufgeprägten 


Schwingungen 


Regelung 


Istströme 
(mit/ohne) 
aufgeprägten 
Schwingungen 


lauf oder IRer 


Messgrößen 
(mit/ohne) 
aufgeprägten 


Gesamtsystem 


Schwingungen 
n 
Yauf Oder Yrer 


Abbildung 4.21: Bestimmung der Übertragungswege 


Mit dieser Vorgehensweise kann folgendes Ergebnis dargestellt wer- 
den, welches den Verstärkungsfaktor a der aufgeprägten Schwingung 
in /, zur Kraftmessstelle der Dämpferbefestigung zeigt. 


Verstarkungsfaktor a [N/A] 


200 400. 
Frequenz [Hz] 


—e— 20. Ordnung 
—©— 40. Ordnung 
—6— 47. Ordnung 
-© 80. Ordnung 
—©- 120. Ordnung 
—8— 240. Ordnung 


600 800 


Abbildung 4.22: Berechnete Verstärkungsfaktoren zum Anbindungs- 
punkt des Dämpfers 
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4.4 Validierung des Modells 


Es zeigen sich die berechneten Werte, welche aus den einzelnen Auf- 
prägungen bestimmter Ordnungen gewonnen wurden und somit das 
Übertragungsverhalten beschreiben. Außerdem ist erkennbar, dass es 
nicht relevant ist, welche Ordnung für die Messung des Übertragungs- 
verhaltens aufgeprägt wurde. 


In Abbildung 4.23 werden die aus dem Experiment berechneten den aus 
der Simulation stammenden Übertragungsfunktionen gegenüberge- 
stellt. 
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Abbildung 4.23: Vergleich der Ubertragungswege 
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4 Gesamtsimulation 


Die Ergebnisse der Validierung für die Gesamtsimulation zeigen Ab- 
weichungen zwischen Experiment und Simulation. Vor allem bei der 
Dämpferkraft ergeben sich sehr große Unterschiede. Die höheren Werte 
aus der Simulation können durch den vereinfachten Modellierungsan- 
satz des Dämpfers erklärt werden, welcher im Modell nur eine ge- 
schwindigkeitsabhängige Gegenkraft erzeugt. So ergeben sich zum 
Beispiel mit einer 6-fachen Verkleinerung der geschwindigkeitsabhän- 
gigen Dämpfung relativ gut übereinstimmende Ergebnisse (Abbildung 
4.24 rechts). Im Experiment ist der Dämpfer aufgrund seines Aufbaus 
am Anschlusspunkt zur Karosserie durch ein Elastomer entkoppelt, 
welches allerdings auch, sowie der Dämpfer selbst, frequenzabhängige 
Eigenschaften besitzt. Dadurch werden die angreifenden Kräfte, vor al- 
lem bei höheren Frequenzen, entkoppelt beziehungsweise weniger stark 
übertragen. 


ohne Korrektur mit Korrektur (empirisch) 6 
102 10? 
— Simulation 


— Simulation 


o0 
OD Messungen do Messungen 


Verstarkungsfaktor [N/A] 


£ 
Verstärkungsfaktor [N/A] 


Q! 

PPE 

o 

O 100 200 300 400 500 600 O 100 200 300 400 500 600 
Frequenz [Hz] Frequenz [Hz] 


Abbildung 4.24: Vergleich der Übertragungsfunktionen zum Dämpfer 
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4.4 Validierung des Modells 


Insgesamt können die teilweise unterschiedlichen Werte durch die fol- 
genden getroffenen Annahmen erklärt werden: 


Vereinfachte Verbindungssteifigkeit des Reifens zur Straße 
Fehlende Strukturdynamik des Reifens durch vereinfachten Mo- 
dellierungsansatz aus nur einer Masse und einer Steifigkeit 
Fehlende Exzentrizität in der elektrischen Maschine 

Annahme des linearen Dämpferansatzes für die Übertragung der 
Kraft ist zu stark vereinfacht; somit fehlt das strukturdynamische 
Verhalten des Dämpfers in den Übertragungsfunktionen 
Annahme des linearen Federansatzes für die Übertragung der 
Kraft ist zu stark vereinfacht; es fehlt das strukturdynamische 
Verhalten der Feder selbst und des Federbeinkopflagers 
Fehlende Frequenzabhängigkeit für Anbindung der Längslen- 
kerbuchse 

Anbindungspunkte zur Karosserie starr (im Experiment Karos- 
serie oder Prüfstandsadapter) 


Trotzdem können mit dem Modell die wesentlichen störenden Effekte 
nachgebildet werden, die in elektrischen Maschinen vorkommen. Dabei 
ist zu beachten, dass Betrag und Phasenlage der Übertragungswege sich 
durchaus unterscheiden, nicht aber die grundsätzlichen Entstehungsme- 


chanismen. 
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5 Identifikation von 
Gegenmaßnahmen 


Aus den bisher dargestellten Informationen und Zusammenhängen kön- 
nen folgende Gegenmaßnahmen abgeleitet werden. 


5.1  Elektromagnetisches Design 


Das elektromagnetische Design ist für die Charakteristik einer elektri- 
schen Maschine ausschlaggebend. Oft werden elektrische Maschinen 
besonders für hohe Leistungen beziehungsweise einer hohen Leistungs- 
dichte ausgelegt, was Material und Kosten spart, wohingegen die vibro- 
akustischen Eigenschaften in dieser Entwicklungsphase oft nicht be- 
rücksichtigt werden. Dies erfolgt meist iterativ mithilfe von elektro- 
magnetischen FEM-Simulationen. Diese Methodik kann aber auch 
angewandt werden, um die auf die Komponenten einer elektrischen 
Maschine wirkenden Kräfte in Hinblick auf ihre vibroakustischen Ei- 
genschaften zu verbessern. Folgende grundsätzliche Möglichkeiten 
sind in der aktuellen Literatur bekannt: 


e Wahl und Auswirkung geeigneten Pol-Nutenzahl-Kombina- 
tion auch in Hinblick auf Vibroakustik [Ver15, Wan13, Isl10, 
Sunl 1, LiZ09] 

e Staffelung im Rotor [Wan16, Junl1] 

e  Schrägung von Rotor oder Statornuten [Blu14, Jur03] 

e Asymmetrische Anordnung der Rotormagnete [Bia02] 

e Optimierung der Rotormagnetausrichtung und Design der 
Permanentmagnete im Rotor [Sun05, Ishl1, Lee16] 

e Optimierung der Statorgeometrie [Pet15, Gurl2, Kro12] 

e Optimierung der äußeren Rotorgeometrie [Put15] 


Eine Ubersicht, wie stark allgemein die vibroakustischen Eigenschaften 
mit der elektromagnetischen Auslegung zusammenhängen und welche 
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5 Identifikation von Gegenmaßnahmen 


Optimierungspotenziale sich beim Anwenden mehrerer Optimierungs- 
möglichkeiten ergeben, zeigt [Kro1l0] oder einschlägige Fachliteratur 
[Gie06]. 


5.2 Änderungen in der mechanischen 
Struktur 


Durch Änderungen in der Geometrie des Radnabenantriebes ist es zu- 
dem möglich das Übertragungsverhalten des mechanischen Aufbaus 
des Gesamtsystems stark zu beeinflussen. Das Ziel, dieses Verhalten im 
Hinblick auf das vibroakustische Verhalten zu optimieren, ist vielver- 
sprechend. Eine Umsetzung dieser Verbesserungen ist allerdings auf- 
grund der ständigen Weiterentwicklung und Änderung der Geometrie, 
z. B. wegen geänderter Montage oder vereinfachter bzw. kostengünsti- 
gerer Produktion der Komponenten, sehr schwierig und auf andere 
Strukturen meist nicht oder nur sehr schwer übertragbar. Sie kann sinn- 
vollerweise nur direkt in der Entwicklungsabteilung durchgeführt wer- 
den. Wie sich aber bereits in den vorherigen Kapiteln gezeigt hat, ist 
zum Beispiel eine achssymmetrische Auslegung des Systems [Str17b] 
für das vibroakustische Verhalten des Systems von Vorteil. 


5.3 Einbringen von Dämpfung in 
Statoplatte und Stator 


Vor allem bei der Anregung der Resonanzstellen beziehungsweise Ei- 
genformen durch die elektromagnetischen Kräfte spielt die Statorplatte 
und der Stator selbst mit dessen Gehäuse eine sehr große Rolle, da diese 
auf die Entstehung kritischer Resonanzen einen erheblichen Einfluss 
hat. Es liegt nahe, dass durch Einbringen von Dämpfung in diesen Kom- 
ponenten kritische Resonanzen gedämpft werden können und so eine 
Verbesserung der vibroakustischen Eigenschaften des Systems erzielt 


102 
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wird. Aufgrund der kompakten Bauweise des Radnabenantriebs ist das 
Einbringen von Dämpfung immer mit konstruktiven Änderungen ver- 
bunden und somit für spätere Versionen meist schwierig übertragbar 
und mit hohen Kosten verbunden. Eine weitere Möglichkeit über die 
Vergussmasse oder dem verwendeten Elektroblech Dämpfung im Sys- 
tem hinzuzufügen, würde nur geringe Vorteile beim speziellen Aufbau 
des Radnabenantriebs aufgrund seiner hohen Polpaarzahl bieten. Eine 
Statorverformung wird hier nur mit hohen Raum-ordnungen angeregt. 
Resonanzen, die dieser Anregungsform entsprechen werden erst bei 
sehr viel höheren Frequenzen angeregt, welche somit unkritisch sind 
und in dieser Arbeit nicht weiter betrachtet wurden. Vielmehr bewegt 
sich der Stator in den kritischen Schwingformen eher als Starrkörper. 


5.4 Änderungen am Fahrwerk 


Das bei den Versuchsträgern verwendete Fahrwerkkonzept ist nicht op- 
timal für das Übertragungsverhalten von vibroakustischen Störsigna- 
len, die aus dem Radnabenantrieb stammen. So werden diese sehr stark 
in den Fahrzeuginnenraum eingeleitet. Ein Grund hierfür ist die not- 
wendige Erhöhung der Feder- und Dämpfersteifigkeiten, welche aus 
fahrdynamischer Sichtweise wegen der erhöhten ungefederten Masse 
durch den integrierten Antrieb im Rad notwendig ist. Dies verstärkt das 
Übertragungsverhalten der Wechselkräfte und leitet diese besser an die 
Karosserie weiter. Eine geänderte Fahrwerkgeometrie kann sich durch- 
aus positiv auf die vibroakustischen Eigenschaften auswirken. Aller- 
dings sind für die Änderung oder Entwicklung eines neuen Fahrwerk- 
kozeptes viele weitere Randbedingungen ausschlaggebend und im Rah- 
men dieses Forschungsprojektes nicht vorhersehbar. 
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5.5 Änderungen der Anbindung des 
Fahrwerks 


Ebenfalls besteht die Möglichkeit direkt an der Anbindung des Fahr- 
werks zur Karosserie eine gezielte vibroakustische Entkopplung für stö- 
rende Frequenzen vorzunehmen. Ebenso wie bei Änderungen im Fahr- 
werk, sind hierbei allerdings viele weitere Randbedingungen zu beach- 
ten, sodass entwickelte Gegenmaßnahmen für andere Fahrwerk- 
konzepte oder Fahrzeugtypen unter Umständen nicht anwendbar sind. 


5.6 Aktive Kompensation durch 
zusätzlichen Aktor 


Mit dem Aufbringen eines zusätzlichen Aktors (z.B. Unwuchtmotor, 
Piezoaktor) können gezielt zusätzliche Schwingungen ins System ein- 
gebracht werden, welche den durch die elektrische Maschine erzeugten 
Schwingungen entgegenwirken und somit das vibroakustische Verhal- 
ten der Maschine verbessern. Notwendig hierfür sind allerdings zusätz- 
liche Komponenten, welche für den späteren Einsatz zusätzliche Kos- 
ten verursachen würden. Hinzu kommt, dass bei konstruktiven Ände- 
rungen im Fahrwerk oder am Radnabenantrieb selbst die 
Randbedingungen stark geändert würden und somit die Wirksamkeit 
neu geprüft werden müsste. 


5.7 Kompensation durch 
Ansteuerungsmethodik 


Durch geänderte Ansteuerungsmethoden lassen sich vor allem die dy- 
namischen und vibroakustischen Eigenschaften einer elektrischen Ma- 


schine positiv beeinflussen, welche in dieser Arbeit besonders von Be- 
deutung sind. Naheliegend ist eine Verbesserung des Geräusch- und 
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5.8 Bewertung der Gegenmaßnahmen 


Vibrationsverhaltens durch Reduzierung der Oberschwingungen im 
Strom. Durch den Wegfall dieser Oberschwingungen werden die in Ka- 
pitel 3.1.1 beschriebenen Phänomene und besonders das von den Dreh- 
momentoberschwingungen verursachte Drehmomentrippel verringert. 
Weiterhin wird häufig versucht, die Drehmomentungleichförmigkeit zu 
reduzieren, um somit das Schwingen am Motor angeschlossener Struk- 
turen zu reduzieren. Durch die Kopplung von mechanischen und elekt- 
romagnetischen Simulationen kann der Zusammenhang zwischen der 
Ansteuerungsmethode und mechanischen Schwingungen aufgezeigt 
und optimiert werden. Eine ausführliche Beschreibung zum Stand der 
Technik dieser Optimierungsmöglichkeit findet sich in Kapitel 6. Vor- 
teil dieser Gegenmaßnahme ist die Anwendbarkeit auch bei geänderten 
Konstruktionen der Mechanik am Radnabenantrieb oder am Fahrwerk 
mit nur leichten Anpassungen. 


5.8 Bewertung der Gegenmaßnahmen 


Als Kriterien bei der Auswahl der untersuchten Gegenmaßnahme wur- 
den folgende Punkte berücksichtigt und bewertet: 


e Übertragbarkeit auf andere Designs 

o Kann die Methode auf eine weiterentwickelte Version des 
Radnabenantriebs übertragen werden? 

o Kann die Methode auf andere Traktionsantriebe im automo- 
bilen Bereich oder auf elektrische Maschinen in anderen Be- 
reichen angewandt werden? 

e Zusätzliche Kosten für Endprodukt 

o Ist für eine Anwendung der Optimierung zusätzliche Hard- 
ware erforderlich (Sensoren, Aktoren, zusätzliche Rechen- 
leistung für die Datenverarbeitung) 

e Forschungs- und Entwicklungsaufwand 

o Wie aufwändig ist die Entwicklung und Implementierung? 

o Wie aufwändig ist eine Anpassung an geänderte Randbedin- 
gungen oder andere elektrische Maschinen? 


105 


5 Identifikation von Gegenmaßnahmen 


Maßnahme Übertrag- Kosten? 
barkeit! 

5.1 Anpassung des elektromagnetischen - + 

Designs 

5.2 Änderungen in der mechanischen F + 

Struktur 

5.3 Einbringen von Dämpfung in Stator- 5 + 

platte 

5.4 Änderungen am Fahrwerk o + 

5.5 Änderungen der Anbindung des o + 


Fahrwerks zur Karosserie 


5.6 Aktive Kompensation durch zusätzli- + E 
chen Aktor 

5.7 Kompensation durch Ansteuerungs- $ n 
methodik 

1) zur Übertragbarkeit: 2) zu Kosten: 

- = schlechte Übertragbarkeit - = erhöhte Kosten für Produkt zu erwarten 
o = mittlere Übertragbarkeit + = keine erhöhte Kosten für das Produkt 


+ = gute Übertragbarkeit 


Tabelle 5.1: Bewertung der Gegenmaßnahmen 


Für die Arbeit wurde sich aufgrund der guten Bewertung (Tabelle 5.1) 
für die Entwicklung einer Ansteuerungsmethodik entschieden, welche 
die durch die elektrische Maschine des Radnabenantriebs eingebrach- 
ten Schwingungen im Gesamtsystem reduziert (Kapitel 6). Gründe 
hierfür sind vor allem die Übertragbarkeit auch auf andere Designs und 
das Potential ohne zusätzlichen Kostenaufwand für das Endprodukt die 
vibroakustischen Eigenschaften zu verbessern. Um ein besseres Ver- 
ständnis für die Effekte zu entwickeln und die Ansteuermethodik 
schnell und effizient testen zu können, wird das in Kapitel 3 und 4 ent- 
wickelte Simulationsmodell und die für die Validierung aufgebauten 
Experimente genutzt. 
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6 Akustisch optimierter 
Steuerungsansatz 


6.1 Allgemeines 


Reduzierung der Stromoberschwingungen 


Es ist bekannt, dass durch die Reduzierung der Oberschwingungen zur 
Grundwelle im Strom die vibroakustischen Störeinflüsse oder die Dreh- 
momentungleichförmigkeit in den meisten Betriebspunkten und für die 
meisten Ordnungen verkleinert werden können. In der Literatur gibt es 
dazu bereits verschiedenen Ansätze: 


e  [Ric15c]: Kompensation deren Stromoberschwingungen 
durch rotorwinkelabhängige Integratoren, dessen Werte im 
Betrieb zu der eigentlich errechneten Spannung addiert wer- 
den und somit zu einer starken Reduzierung führen. Kein Be- 
zug zur Vibroakustik oder Drehmomentrippel. 

e [Lul5, Qil5] : Kompensation der Stromharmonischen durch 
Aufprägen zusätzlicher Harmonischer in der Spannung. Bezug 
auf das Drehmomentrippel, aber keine Messung zur tatsächli- 
chen Reduzierung. 

e [Pell2b]: Aufzeigen des vibroakustischen Potenzials durch 
Verringerung der Stromoberschwingungen. Reduzierung der 
Ströme durch paralleler Regelung, welche die Harmonischen, 
ähnlich der feldorientierten Regelung, minimiert. 

e [Geo13]: Reduzierung der Vibrationen durch Kompensation 
der Stromharmonischen über Beobachter (ähnlich [Kar16]) 
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6 Akustisch optimierter Steuerungsansatz 


Verbesserung der vibroakustischen Eigenschaften durch Reduzierung 
des Drehmomentrippels mit Hilfe der Stromformung 


In der Literatur finden sich bereits Arbeiten, die durch gezielte Beein- 
flussung der Motorströme das vibroakustische Verhalten verbessern, 
indem sie diese gezielt abweichend zur Grundwelle formen und somit 
zusätzlich die durch die aus vibroakustischer Sicht ungünstige elektro- 
magnetische Auslegung eingebrachten Schwingungen versuchen zu 
kompensieren. Sie beschränken sich aber vor allem auf die Kompensa- 
tion der Drehmomentungleichförmigkeiten. Die dafür angewandten 
Methoden sind jedoch durchaus unterschiedlich, auch wenn sie letzt- 
endlich alle auf dem Prinzip der Stromformung basieren. Vor allem 
wird versucht, die Drehmomentschwingungen im Zeitbereich durch 
Rückkopplung einer Drehmomentmessstelle oder mit Hilfe eines geeig- 
neten Modells zur Berechnung des Moments zu reduzieren. 


e [Fen18]: Drehmomentrippelkompensation durch die Nutzung 
von Drehschwingungsharmonischen 

e [Str17a]: Untersuchung verschiedener Methoden zur Realisie- 
rung einer aktiven Vibrations- und Geräuschauslöschung an 
einem elektrischen Radnabenantrieb 

e [Bral3b]: Reduzierung des Drehmomentrippels vor allem im 
Sättigungsbereich durch einen iterativen Lernprozess 

e [Nakl1]: Drehmomentkompensation mit Online-Berechnung 
der Ungleichförmigkeit aus einem Modell durch entsprechen- 
des Entgegensteuern 

e [Cha99]: Vergleich von Stromformen im Hinblick auf ver- 
schiedene Kriterien bei Motoren mit Oberflachenmagneten 
(zum Beispiel: „Minimum Torque Ripple“) 

e [Ho196]: Identifizierung der Form der Drehmomentungleich- 
förmigkeiten durch ein Modell, welches diese nachbilden kann 
und Anpassung der Ströme 
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6.1 Allgemeines 


Verbesserung der vibroakustischen Eigenschaften durch Stromfor- 
mung, welche sich nicht auf die Reduzierung des Drehmomentrippels 
beschränken 


Ebenfalls ist es möglich durch gezieltes Einbringen von Schwingungen 
in den Motorströmen (engl. Harmonic Current Injection; Abkürzung: 
HCI) die vibroakustischen Eigenschaften positiv zu beeinflussen. 


e [Harl8]: Radialkraftreduzierung und dadurch Reduzierung der 
Luftschallabstrahlung durch HCI 

e [Wan18]: Auswirkung geänderter Stromform auf die Radial- 
kräfte 

e [Oks16]: Verschiedene Methoden zum Aufprägen von 
Schwingungen werden untersucht. Der Schwerpunkt der 
Kompensation liegt auf der Reduzierung des Geschwindig- 
keitrippels, da dieses über die vorhandene Sensorik erfasst 
werden kann 

e [Ben14, Ben15]: Kompensation von Vibrationen einer ange- 
schlossenen Getriebestufe mit rückwirkenden Sensor und Ac- 
tive Vibration Control Algorithmus 

e [Kroll]: Multikriterielle Optimierung; unter anderem Dreh- 
moment und Radialkräfte durch Stromformung 


Die bisherigen Forschungsarbeiten können in Systemen, die mit und 
ohne vibroakustischer Rückwirkung arbeiten, aufgeteilt werden. Allge- 
mein wird die Möglichkeit mit rückwirkender Größe zur Kompensation 
vibroakustischer Störungen Active-Vibration-Control genannt (AVC). 
Diese Verfahren bieten sowohl Vor- und Nachteile. 


Vorteile: 


e Möglichkeit zur Kompensation von Schwingungen anderen 
Ursprungs (z. B. Verzahnungsanregung im Getriebe) 

e Möglichkeit auf geänderte Betriebszustände zu reagieren 

e Detektion anderer Informationen über das System 
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6 Akustisch optimierter Steuerungsansatz 


Nachteil: 


e Zusätzliche Sensorik für vibroakustische Rückwirkung not- 
wendig (Aufwand und Kosten) 


Aufgrund der fehlenden Sensorik oder Anbindungsmöglichkeiten der 
Sensorik im Gesamtsystemen der meisten Antriebssysteme, kon- 
zentriert sich diese Arbeit auf Verfahren ohne rückwirkendende vibro- 
akustische Informationen. 


Abgrenzung 


In dieser Arbeit werden die harmonischen Schwingungen direkt auf die 
Sollgrößen für die Ersatzströme I, und I, im rotierenden Koordinaten- 
system aufgeprägt. Im Gegensatz zu den häufig in der Literatur genutz- 
ten Methoden der Aufprägung auf die Spannung, was auch oft im Drei- 
phasensystem geschicht, bietet dies folgende Vorteile: 


e Keine Umrechnung der Ordnungen zwischen Wirkung und Auf- 
prägung notwendig, im Gegensatz zum Aufprägen im Dreipha- 
sensystem 

e Bezug auf den Strom liefert gleichzeitig Informationen über das 
magnetische Feld und den magnetischen Fluss, da diese vonei- 
nander abhängig sind 

e Durch starke lineare Zusammenhänge im rotierenden Koordina- 
tensystem können die aus der Magnetik und den Stromober- 
schwingungen stammenden Störungen einfacher getrennt wer- 
den 

e Ableitung einer Betriebsstrategie und Interpolation zwischen 
den getesteten Betriebspunkten ist möglich 


Zu dieser Art der optimierten Ansteuerung in Bezug auf deren Auswir- 
kungen in einem komplexen mechanischen Gesamtsystem in Versuch 
und Simulation sind keine aktuellen Forschungen bekannt. Auch wei- 
tergehende Fragestellungen in Bezug auf die Reaktion der Regelung auf 
die Sollvorgaben mit Oberschwingungen wurden noch nicht untersucht. 
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6.1 Allgemeines 


Somit wird in dieser Arbeit besonderes auf die harmonischen Aufprä- 
gungen im Soll-Stromsignal im rotierenden System eingegangen und 
wie sich diese in einem komplexen Gesamtsystem auf die vibroakusti- 
schen Eigenschaften auswirken. Ebenso wird dargestellt, wie sensitiv 


diese Methodik in Bezug auf Änderungen im mechanischen System ist 
und welche sinnvollen Betriebsstrategien dabei angesetzt werden kön- 


nen. 


Ansatz 


Die zusätzlich aufgeprägten Schwingungen werden den konstanten 
Strömen J, und Iq mit den Amplituden Aj, Ag und den Phasenwinkeln 


Pq. Pa für die zu kompensierende Ordnung v hinzugefügt. 


Ig = Iq + Aq’ sin(v: 6 + Gq) 


Ty" = la +Aa' sin(v -0 + Pq) 


(6.1) 


(6.2) 


Dargestellt werden kann dieser Zusammenhang auch im Signalfluss- 


bild: 
g 
x 
vo 
> 
Ae 
Yq + sin | 
>| 
Ag a E ] 
» Sollvorgabe 
u ———— 
Iq » * I” 
d 
> 
ie 
Pa + sın 
>| 
Ag >| x 
> Sollvorgabe 
I + -—> 
q > I * 


Feldorientierte 
Regelung 


Stellgröße 
Tr — 
Ud 


Stellgröße 


q 


Ug 


Abbildung 6.1: Signalflussbild der Aufprägemethode 
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6 Akustisch optimierter Steuerungsansatz 


Der Zusammenhang der Strangströme Jy mit den aufgeprägten Schwin- 
gungen kann mit Hilfe der Gleichung (2.3) folgendermaßen dargestellt 


werden: 
Ty = Iyeruna + luober 
Mit 
Iyeruna = 6080 I, — sin 0 I, 
und 


Iyober = Aa ' cos 0  sin(vO + pa) — Ag: sind: sin(vé + P4) 
ergibt sich 


Iyoper = 5Aalsin(@(v — 1) + gq) + sin + 1) + pa4)] - 


= Aqlcos(0@ -1)+ P4) = cos(0 (v +1)+ Pa)l (6.3) 


Es ist zu erkennen, dass sich die aufgeprägte Schwingung im rotorfes- 
ten Koordinatensystem von der Ordnung v in zwei Ordnungen v + 
1 und v — 1 im statorfesten Koordinatensystem projiziert. 


Ebenso können hier zwei Spezialfälle gebildet werden. 


Aufprägung der Schwingungen ausschließlich in I4 => Ag = 0 

1 
luoper == 5 Aqlcos(8(v -1)+ P4) — cos(O(v + 1) + P4)l (6.4) 
Durch diese Methode werden zwei gleich große Schwingungen in den 
Ordnungen v + 1 und v — 1 aufgeprägt, welche sich lediglich in der 


Phase um 9, um 7 unterscheiden. 


Aufprägung der Schwingungen ausschließlich in I, => Ag = 0 


luover = “Aa [sin(6@(v — 1) + pa) + sin(6(v + 1) + pa)] (6.5) 
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6.1 Allgemeines 


Durch diese Aufprägung werden zwei gleich große Schwingungen mit 
gleicher Phase in den Ordnungen v + 1 und v — 1 moduliert. 


dq, 
Koordinaten 
O o 2% 30, 4w, 5w 6% 7a 8o 9w, 10w 11w, 12w, 13w 


we wf bt Pe, 


Koordinaten O Wy 20% 30, 40, 50 6% 7% Bos Wy 10% 11w, 120, 130, 


Abbildung 6.2: Veranschaulichung der Transformation 


Um die Parameter für die aufgeprägten Schwingungen Ag, Aq, Pq, Pa 
zu berechnen, welche das vibroakustische Verhalten positiv beeinflus- 
sen, können für Versuch und Simulation unterschiedliche Vorgehens- 
weisen herangezogen werden, welche in den nachfolgenden Kapiteln 
beschrieben werden. 
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6 Akustisch optimierter Steuerungsansatz 


6.2 Anwendung der Methodik in der 
Simulation 


Für die Anwendung von aufgeprägten Schwingungen in der Simulation 
wurde das in Kapitel 4 entwickelte Modell verwendet. Grundsätzlich 
bietet sich so die Möglichkeit, die Amplitude und Phasenlage der auf- 
geprägten Schwingungen zur Kompensation rechnerisch zu bestimmen. 


6.2.1 Darstellung ohne am Rotor wirkende 
Kräfte 


In dem Kapitel soll gezeigt werden, wie sich aufgeprägte Schwingun- 
gen grundsätzlich verhalten und welche Kompensationsmöglichkeiten 
sich daraus ergeben. Dazu werden die auf den Rotor wirkenden Kräfte 
vorerst nicht berücksichtigt, da diese einen geringen Einfluss auf das 
vibroakustische Verhalten haben. Erst in einem nachfolgenden Kapitel 
werden auch die auf den Rotor wirkenden Kräfte mit betrachtet. 


Die Zusammenhänge zur Berechnung der Kompensation können bei 
Vernachlässigung der Kräfte auf den Rotor je nach betrachteter Ord- 
nung vereinfacht werden. Um die zur Kompensation notwendigen Pa- 
rameter zu ermitteln, werden die an den Verbindungspunkten des Fahr- 
werks zur Karosserie auftretenden Störkräfte Fp durch die elektromag- 


netischen Übertragungswege Gy geteilt. 


p Fr _ Sisil—Fa Gai Fa Cz] 6.6) 
Iq ~~ Aq I 1 : 
Ge Dilz Gzir + G2, CZF 


Somit kann A, und p, aus P;, in Abhängigkeit von Drehmoment und 
Drehzahl berechnet werden: 


Aq =|Piq| und P4 = Pig) 
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6.2 Anwendung der Methodik in der Simulation 


Wendet man die berechneten Parameter der aufgeprägten Schwingun- 
gen für /, an, ergibt sich daraus die vollständige Kompensation der stö- 
renden Wechselkräfte am betrachteten Ort, wie Abbildung 6.3 zeigt. 
Dabei kompensieren sich die Kräfte aufgrund von aufgeprägten 
Schwingungen und den aus der magnetischen Topologie stammenden 
Schwingungen vollständig. Da diese Kompensationsparameter dreh- 
zahlabhängig sind, sind diese im Diagramm für ein ausgewähltes Dreh- 
moment dargestellt. 


200 Nm, 40. Ordnung, Ausgang: Dämpfer 
20 F T T T 1 — 


J m .d i R —— ohne HCI 
„mit HCI“ vollständig kompensiert et 


a 
= 
= 10} | 
m 
x 
0 L | 


100 200 300 400 500 
Frequenz [Hz] 


a 
œ 
© 
Phasenlage [°] 


1 L 1 1 1 L 0 
200 300 400 500 600 700 800 


Drehzahl [rpm] 


Abbildung 6.3: Kompensation der am Dämpfer wirkenden Kräfte zur 
Karosserie mit dessen Kompensationsparameter für 74 


Ebenfalls ergibt sich daraus die Möglichkeit, die Kräfte mit Hilfe des 
la-Stroms zu kompensieren. Allerdings zeigt sich in Abbildung 6.4, 
dass die Übertragungswege sehr viel weniger stark ausgebildet sind. 
Die hohen Amplituden sind in der Realität sehr schwierig umzusetzen, 
vor allem bei hohen Frequenzen. Für die Berechnung der für die Kom- 
pensation notwendigen Parameter gilt analog zu Formel (6.6): 


48 t pmst _ pr pmr 
= Fr 2 li Grip Fzi Cnr | (6.7) 
Id — "Id 548 [pid mSt Id mSt ; 
Gr le * Gap + GH, GRE 
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6 Akustisch optimierter Steuerungsansatz 


Aa = [Pial und 94 = (Pra) 


200 Nm, 40. Ordnung, Ausgang: Dämpfer 


N 
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= 
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oo 
© 
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Abbildung 6.4: Kompensation der am Dämpfer wirkenden Kräfte zur 
Karosserie mit dessen Kompensationsparameter für J4 


Aufgrund der sehr schwach ausgeprägten Übertragungswege wird für 
weitere Betrachtungen die Aufprägung in I4 nicht weiterverfolgt. 


In einer weiteren Darstellung (Abbildung 6.5) sind mehrere Ordnungen 
in einem Betriebspunkt gleichzeitig kompensiert und anschließend die 
Ausgabegrößen bei Drehzahl 100 rpm (Abbildung 6.6) genauer be- 
trachtet. Dabei ist zu erkennen, dass den störend wirkenden Kräften auf- 
grund des betrachteten Drehmoments die Hauptursachen zugeordnet 
werden können. 
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6.2 Anwen 


dung der Methodik in der Simulation 


200 Nm, 40. Ordnung, Ausgang: Dämpfer 
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7 aus Radialkräften 
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Abbildung 6.5: Kompensation am Dämpfer wirkenden Kräfte mehre- 
rer Ordnungen 
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6 Akustisch optimierter Steuerungsansatz 


Im Drehmomentverlauf (unten) zeigt sich, dass vor allem die 120. und 
240. Ordnung kompensiert werden, da diese Drehmomentschwingun- 
gen hauptsächlich die wirkenden Kräfte am Anbindungspunkt (hier: am 
Dämpferanbindungspunkt) verursachen. Um auch die 40. und 80. Ord- 
nung der wirkenden Kräfte am Anbindungspunkt kompensieren zu kön- 
nen, werden zusätzliche Drehmomentungleichförmigkeiten aufgeprägt, 
welche dann dort kompensierend wirken. Der für die Kompensation 
notwendige /,-Strom kann in Abbildung 6.6 betrachtet werden. Er zeigt 
die für die Drehmomentbildung notwendigen Oberschwingungen. 
Diese vibroakustisch kompensierend wirkenden Oberschwingungen 
können zu zusätzlichen Verlusten in der elektrischen Maschine führen. 
Diese werden in Kapitel 6.4.2 genauer beschrieben. 


0 5 10 15 20 25 30 35 
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190 L 1 L 
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N 
© 
oO 


Abbildung 6.6: Ausgangsgrößen zu der Kompensation mehrerer Ord- 
nungen 
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6.2 Anwendung der Methodik in der Simulation 


6.2.2 Darstellung mit allen wirkenden Kräften 


Im folgenden Abschnitt werden die durch den Rotor verursachten stö- 
renden Kräfte mitbetrachtet. 


Allgemein können die Zusammenhänge zur Berechnung der Kompen- 
sationsparameter vereinfacht beschrieben werden mit: 


Störende Kräfte 


Piq = Uebertragungsfunktion auf gepraegte Schwingung 


(6.8) 


Dazu werden die störenden Kräfte und die Übertragungsfunktion für 
die aufgeprägte Schwingung berechnet. Für die Übertragungsfunktion 
werden nun alle vibroakustisch relevant wirkende Kräfte auf Stator und 
Rotor betrachtet. 
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Abbildung 6.7: Kompensation 120. Ordnung (Dampfer und LLBz) 
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6 Akustisch optimierter Steuerungsansatz 


Mit dieser Methode können Ordnungen, die nicht Vielfache der 40. 
Ordnung sind, also von Kräften, die auf den Rotor ausgeübt werden, 
dargestellt und kompensiert werden (Abbildung 6.8). 
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200 Nm, 48. Ordnung, Ausgang: LLBz 


> —— ohne HCI 
Z zl —mitHCI |] 
ka 
w 
m 
eS 
0 
100 200 300 400 500 600 
Frequenz [Hz] 
— 360 — 
270 © 
Gy 
4180 = 
o 
Joo 2 
I „= 0 £ 
600 700 


Drehzahl [rpm] 


Abbildung 6.8: Kompensation der 47. und 48. Ordnung an LLBz 


Aus den Ergebnissen der Kompensation in der Simulation können fol- 
gende Aussagen getroffen werden: 
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6.2 Anwendung der Methodik in der Simulation 


e Hauptverursacher der vibroakustischen Störungen ist das 
Drehmomentrippel. Dieses kann durch das Aufprägen harmo- 
nischer Schwingungen gut kompensiert werden. 

e Auch die vibroakustischen Störungen, die nicht aus dem Dreh- 
moment stammen, können durch ein zusätzliches Aufprägen 
von Harmonischen in den I,-Strömen kompensiert werden. 


e Schwingungen in Jg-Strémen aufzuprägen, um eine Kompen- 
sation der Schwingungen zu realisieren, ist im System Radna- 
benantrieb nicht zielführend, da die resultierenden Übertra- 
gungsfunktionen zu den Anbindungspunkten zu schwach sind. 
Somit müssten sehr hohe Amplituden gestellt werden. Jedoch 
kann dieses Kompensationspotential eine Rolle spielen, wenn 
es um direkte Schallabstrahlung vom Statorgehäuse geht. 

e Die berechnete Amplitude in den /,-Strömen liegt in allen Be- 
triebspunkten unter 5 A. Dies ist ein Wert, welcher mit dem 
vorhandenen Regler realistisch anwendbar ist und zudem 
keine zusätzlichen hohen Verluste beziehungsweise starke Re- 
duzierung des Wirkungsgrades zur Folge haben wird. 

e Die für die Kompensation notwendigen Parameterwerte unter- 
scheiden sich in den jeweiligen Betriebspunkten (Drehzahl 
und Drehmoment). 

e Aufgrund der in der Simulation vereinfacht angenommenen 
Randbedingungen und Modellierungsparameter können die 
berechneten Werte für Amplitude und Phasenlage nicht ein- 
fach für das Experiment übernommen werden. 


Aus diesen Gründen wurde für die Forschungsarbeit eine Methodik ent- 
wickelt, um die für die Kompensation notwendigen Parameter im Ex- 
periment schnell und einfach ermitteln zu können (siehe Kapitel 6.3) 
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6.3 Anwendung der Methodik am 
Prüfstand 


Die Ergebnisse der Simulation zeigen eine grundsätzliche Übereinstim- 
mung mit dem Experiment, um die Entstehungsmechanismen verstehen 
zu können, wie in Kapitel 4.4 beschrieben. Jedoch führt die notwendige 
Vereinfachung der Simulation zu Abweichungen in Amplitude und 
Phasenlage, welche große Auswirkungen auf die zu kompensierenden 
Größen bei aufgeprägten Schwingungen haben können. Deswegen wird 
in dieser Arbeit eine Messmethodik zur Parametersuche vorgestellt, 
welche es ermöglicht, auch ohne Information aus der Simulation die 
Parameter für eine Kompensation der störenden Schwingungen zu be- 
stimmen. 


Dieses Messverfahren wurde in drei verschiedenen Stufen angewandt: 
e Motorenprüfstand (Beschreibung Aufbau im Anhang): 


An einem Motorenprüfstand wurden die grundsätzlichen Ef- 
fekte, welche bei einem geänderten Steuerungsansatz mit HCI 
auftreten, untersucht. Dabei wurden Schwingungen im Drehmo- 
ment und Beschleunigungen am Statorgehäuse gemessen und 
kompensiert. 


e Prüfstand Fahrwerk (Beschreibung Aufbau im Anhang): 


Die zuvor ausführlich getestete optimierte Ansteuermethode 
wurde an einem im Projekt entwickelten Prüfaufbau mit Fahr- 
werk getestet. Dabei wurden die Kräfte am Befestigungspunkt 
des Dämpfers und der Feder gemessen. Ebenso lassen sich 
Rückschlüsse auf die an der Längslenkerbuchse angreifenden 
Kräfte über die Analyse der Signale der dort angebrachter Be- 
schleunigungssensoren ziehen. Aufgrund der dort zur Verfü- 
gung stehenden Mittel konnten nur Betriebspunkte mit sehr ge- 
ringem Drehmoment getestet werden. 


122 


6.3 Anwendung der Methodik am Prüfstand 


e Gesamtfahrzeug (Beschreibung Aufbau im Anhang) 


Um die grundsätzliche Wirksamkeit der in dieser Arbeit entwi- 
ckelten Methodik zu verifizieren, wurde dieser Ansatz auch im 
Gesamtfahrzeug getestet. Aufgrund der zur Verfügung stehen- 
den Mittel kommt der Nachfolger (Version „gamma‘“‘) des zuvor 
getesteten Radnabenantriebs zum Einsatz. 


6.3.1 Messablauf 


Der grundsätzliche Messablauf für Experimente mit aufgeprägten 
Schwingungen kann folgendermaßen beschrieben werden. Die Parame- 
ter werden immer nur nach einer vollen mechanischen Umdrehung des 
Radnabenantriebes geändert. Ebenso wird das betrachtete Fenster, wel- 
ches später für die Auswertung in den Frequenzbereich übertragen und 
ausgewertet wird, genau auf eine mechanische Umdrehung angepasst. 
Daraus ergeben sich folgende Vorteile: 


e Ordnungsauflösung von einer mechanischen Ordnung 
e Keine „Leakage“ -Effekte durch Fensterung oder Auflösung 


Anschließend wird eine Umdrehung abgewartet, sodass sich der neu 
eingestellte Istwert einschwingen kann. Zwischen Änderungen der Pa- 
rameter für die aufgeprägten Schwingungen wird eine Referenzmes- 
sung gemacht, sodass auch eine Änderung des vibroakustischen Ver- 
haltens des Systems während eines Testlaufes bemerkt werden kann. 
Somit werden immer abwechselnd zwei Umdrehungen mit und zwei 
Umdrehungen ohne HCI abgefahren. Abbildung 6.9 veranschaulicht 
diese Vorgehensweise. 
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Messablauf für eine aufgeprägte Ordnung 
L 


Einschwingphase 
Referenzmessung 
Referenzmessung 

Einschwingphase 
aufgeprägte 
Schwingung 

Messung 
aufgeprägte 
Schwingung 
Einschwingphase 
Referenzmessung 


Strom [A] (/, oder I,) 


>t 
i i+1. i+2. i+3. i+4. i+5. 
Umdrehung Umdrehung Umdrehung Umdrehung Umdrehung Umdrehung 


Abbildung 6.9: Messablauf für Parametersuche 


6.3.2 Methode zur Parametersuche 


Vorversuche haben gezeigt, dass die Bestimmung der Parameter für 
HCI ohne geeignete Methode sehr aufwändig ist. Vollvektorielle Ver- 
suchspläne benötigen sehr lange Prüfstandszeiten. Somit ist an größe- 
ren Prüfständen die mehrdimensionale Parametersuche mit sehr hohem 
Aufwand und Kosten verbunden oder aufgrund der Parameteranzahl 
nicht mehr durchführbar. Auch beim Problem der Parametersuche kann 
die hohe Anzahl an Parameter und Dimensionen nur mit Hilfe von ef- 
fektiveren Methoden oder der Simulation gelöst werden. So ist es mög- 
lich, Hinweise auf eine Vereinfachung oder Verkleinerung der Dimen- 
sionen zu finden. 


Die zu suchenden Parameter sind: 


e Amplituden der aufgeprägten Schwingung Ayjycı 
e Phasenlagen der aufgeprägten Schwingung ryc 
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Diese können eine Abhängigkeit von folgenden Faktoren besitzen: 


e Ordnung 

e Drehzahl 

e Drehmoment 

e Betrachtete zu kompensierende physikalische Größe 


Hinzu können Änderung der mechanischen Übertragungswege und 
Randbedingungen durch folgende Faktoren kommen: 


e Änderungen der Dämmungseigenschaften der Materialien im 
Fahrwerk aufgrund Temperatur (z. B. Längslenkerbuchse) 

e Radlast 

e Fehlstellungen des Motors 


Aus diesem Grund wurde ein Verfahren zur Parametersuche von auf- 
geprägten Schwingungen entwickelt, welches im Nachfolgenden be- 
schrieben wird. Hiermit ist es möglich Phasenlage und Amplitude in 
nur wenigen Schritten zu berechnen. Grundsätzlich beruht diese Me- 
thode auf dem mathematischen Zusammenhang der Überlagerung 
zweier Schwingungen in linearen Systemen. Es wird also von einer li- 
nearen Modellvorstellung ausgegangen. 


Abbildung 6.10 zeigt die störende Schwingung, welche sich aus den 
Effekten der Elektromagnetik inklusive mechanischer Übertragungs- 
wege ohne Stromoberschwingungen rz.s zusammensetzt und den Stro- 
moberschwingungen lyc; mit deren Übertragungswege Gjycı. Die re- 
sultierende Schwingung an der Messstelle wird mit y bezeichnet. 


Ina | Gmcı(S) s 
y 


Tges 


Abbildung 6.10: Darstellung des zu betrachtenden Systems 
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Das System bei Betrachtung nur einer Ordnung kann somit auch fol- 
gendermaßen als Zeiger beschrieben werden: 


Aracı ` lGrucr| ; et Pracrt#6ıncn + Arges" el Prges 6.9) 
= Ay - el Gy i 
mit 
i=vV-1 
und 
Incr = Arner: e'Plncı 
thes = Arges . el’Prges 
y=Ay:e'® 


Daraus folgt: 


2 
(ArncılGincıl cos(Prncr + $Gıncı) + Arges COS Prges) F 


F 
y ; . 2 
(ArncılGincıl sin(@yycr + 4GrHcr) + Arges Sin Prges) (6.10) 
A G sin + 4G +A sin 
er ( ract|Giycr| Sin(Pracı uct) + Arges ae (6.11) 
ArncılGracıl Cos(Pracı + &Gıncı) + Arges COS Prges 


Die Unbekannten der Gleichungen |G;acıl, Prmcr Ar, 9, können somit 
nach zwei Messabläufen bestimmt werden. Der Messablauf besteht aus 
einer Messung ohne aufgeprägte Schwingungen, also Ajycı = 0, wel- 
che den Referenzwert für y bildet und einer Messung mit aufgeprägten 
Schwingungen (Aycı # 0). Anschließend kann das Gleichungssystem 
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gelöst werden und alle Unbekannten |Gyacıl, mcr Ar, Yr werden be- 
rechnet. Danach werden die gesuchten Parameter Ajycr und Pjucı“ für 
die Kompensation und somit die Reduzierung der Vibrationen berech- 
net. 


Ay(Arucı » Pmr” |Gincrl, Prucy Ar, Pr) = 0 (6.12) 


In der nächsten Messung werden die berechneten Parameter A,ycı“ und 
PrHcı genutzt, um das erwartete Ergebnis zu verifizieren und dessen 
Genauigkeit zu verbessern. Anschließend stehen mehr Messwerte zur 
Verfügung als für die Berechnung notwendig. Hierfür wird das Least- 
Square-Verfahren genutzt. Dabei wird der quadratische Fehler zwi- 
schen dem berechneten Wert der aufgeprägten Amplitude A, und Pha- 
senlage yy und den dazugehörigen Erwartungswerten E(A,) und 
E(p,,) durch Variation der geschätzten Übertragungsfunktion mit deren 
Vergrößerungsfunktion |G;ycı| und deren Phasenverschiebung @rucı 
und der geschätzten Störung, beschrieben durch A, und @,, minimiert. 


YA, - EA) 


ie ii a (6.13) 
THCILPIHCL4r, Pr Xo = E(p D? 
y V. 


Abbildung 6.11 zeigt den gerade beschriebenen Prozess der Parameter- 
bestimmung in einer Ablaufskizze. 
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frei gewählt: 
Amici 
> — Messung 1 
Pracı 
yı Übergabe der gemessenen Größe 


v 


Auflösen der Berechnung der 
Gleichung Kompensationsparameter 


A i 
> THCI — Messung 1+i 


Yı+i Übergabe der gemessenen Größe 
Y 


Berechnung der 


a ee > Berechnung LSQ Parameter durch 
Minimierung der Fehler 


i=i+1 nein 


Nächster 


Kompensation : 
Betriebspunkt 


optimal? 


Abbildung 6.11: Ablauf fiir Parametersuche 


Abbildung 6.12 zeigt die fiir das Forschungsprojekt entwickelte 
LabView-Messroutine mit Teilen des dazugehörigen Matlab-Codes (in 
der Mitte der Abbildung), der es ermöglicht direkt nach jeder Messung 
die geschätzte Fläche zu sehen, wie in Abbildung 6.13 (rechts) gezeigt 
wird. 
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Abbildung 6.12: Messroutine in LabView mit Matlab-Code 


Als Beispiel für die Gültigkeit dieser Vorgehensweise ist in Abbildung 
6.13 ein Vergleich der geschätzten Fläche (rechts oben), welche aus der 
vorgestellten Methodik gewonnen wurde, mit einer gemessenen Fläche 
aus 3600 Messabläufen (links oben) dargestellt. Zudem wird unten 
links der Ablauf der Messungen dargestellt, woraus sich beim Untersu- 
chen eines größeren Drehzahlbereiches das Diagramm unten rechts 
ergibt. Daraus lässt sich das Potenzial zur Kompensation ableiten. 
Übertragen wird die Messung mit der besten Kompensation des be- 
trachteten Messpunktes. 
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Geplottete Fläche Geschätzte Fläche 


Über 3600 Messungen mit 3 Messungen 
(nach dargestellter Methodik) 


n 
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Abbildung 6.13: Darstellung der Vorgehensweise beim Suchen der 

Parameter und Vergleich mit den gemessenen und geschätzten Flä- 

chen (oben) an der Längslenkerbuchse in x-Richtung im Betriebs- 
punkt 68 km/h der 40. Ordnung 


Abbildung 6.13 zeigt somit den Prozess der Bestimmung der Parame- 
ter. Die Anfangsamplitude A;ycı = 1,1 A und Phase #;ycı = 50 ° sind 
frei gewählt, um Ajycı und ryc“ berechnen zu können. Nach der 
zweiten Berechnung ist zwar schon eine Reduzierung der Vibrationen 
mit aufgeprägten Schwingungen messbar, allerdings ist in diesem 
Schritt noch nicht die erwartete fast vollständige Kompensation einge- 
treten. Der Gründe hierfür können Messungenauigkeiten des Sensors 
und der Messwerterfassung sowie Ungenauigkeiten im Timing der 
Messungen sein. In diesem Beispiel sind ab der 3. Messung die Ergeb- 
nisse der berechneten Parameter ausreichend gut, um die Vibrationen 
auch spürbar zu reduzieren. 
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6.3.3 Kompensationsmessungen 


Die Parametersuche und die daraus erfolgten Kompensationsmessun- 
gen wurden am System des Radnabenantriebs mit Fahrwerk am Akus- 
tik-Allrad-Rollenprüfstand (AARPS) des Instituts für Fahrzeugsystem- 
technik (FAST) getestet. Dabei wurde die vorgestellte Methode zur Pa- 
rametersuche für verschiedene Ordnungen o mit verschiedenen 
Drehzahlen n angewendet. Unter anderem wurden als Ausgänge die 
Kraft am Dämpfer, Feder und die Beschleunigungen in x- und z- Rich- 
tung an der Längslenkerbuchse gemessen, wie Abbildung 6.14 zeigt. 


A Kraftmessdose an Dämpfer 


Kraftmessdose an Feder 


Beschleunigungssensor 
in x- und z- Richtung 
(LLBx und LLBz) 


Abbildung 6.14: Lage der ausgewerteten Sensoren 


Ausgewählt zur Kompensation wurden die am stärksten ausgeprägten 
Ordnungen (20., 40., 47., 80., 120., 240.). Getestet wurde die Kompen- 
sation in einem akustisch interessanten Betriebsbereich von 6-80 km/h. 
Beispielhaft dargestellte Messungen finden sich nachfolgend für die 
Beschleunigung an der Längslenkerbuchse in z-Richtung (Abbildung 
6.15). 


Grundsätzlich ist in den Messungen zu erkennen, dass eine sehr große 
Amplitude für aufgeprägte Schwingungen in I, gestellt werden muss, 
um Auswirkungen auf das vibroakustische Verhalten in Form von Be- 
schleunigungen oder wirkenden Kräften messen zu können. Dies zeigte 
sich bereits in der Simulation (Kapitel 6.2.1). Hinzu kommt, dass es 
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Koppeleffekte zwischen den Steuergrößen gibt, sodass bei aufgepräg- 
ten Sollströmen in J, mit der gleichen Frequenz Stromschwingungen in 
Iq messbar waren und umgekehrt. Eine genauere Analyse der Übertra- 
gungswege ist in Kapitel 6.3.5 zu finden. Deswegen sind nachfolgend 
ausschließlich Kompensationsmessungen mit aufgeprägten Schwingen 
in I, dargestellt. 


In Abbildung 6.15 sind im oberen Diagramm die einzustellenden Para- 
meter für HCI zu finden. Aufgetragen sind die Ist- und Sollwerte der 
gestellten Oberschwingungen der 47. Ordnung für die im unteren 
Schaubild dargestellte Kompensation (grün). Die orangene Kurve in 
diesem Diagramm zeigt die Amplitude der Beschleunigung ohne HCI 
mit ihren Einzelmessungen. Die Menge der Einzelmessungen entstehen 
aufgrund der Parametersuche, die vor jeder Messung mit HCI eine Re- 
ferenzmessung durchführt. 
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Abbildung 6.15: Parameter der HCI für Längslenkerbuchse (z) bei ei- 
ner Kompensation mit I, Strömen 
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Ähnliches Kompensationspotential findet sich in den weiteren Ordnun- 
gen ebenfalls wieder (Abbildung 6.16). Nur das Kompensieren der 20. 
Ordnung funktioniert nicht zufriedenstellend. Grund hierfür ist aller- 
dings die nicht periodische Anregung der 20. Ordnung, da hier 
Amplitude und Phasenlage variieren. Die Ursache liegt im verwendeten 
Winkelmessverfahren. 
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Abbildung 6.16: Übersicht des Potentials zur Kompensation an der 
Längslenkerbuchse (z) 


Weitere Messungen für andere Messpunkte am Prüfaufbau finden sich 
im Anhang wieder (Abbildung 7.5, Abbildung 7.6, Abbildung 7.7). 


Ebenso wurde am Motorenprüfstand mit dieser Methodik die Kompen- 
sation der Beschleunigungen an der Statorplatte in axialer und radialer 
Richtung und des Gehäuses (Abbildung 7.8) sowie des Drehmomentes 
über die Messwelle untersucht. Es zeigte sich, dass mit der Aufprägung 
von Harmonischen in den /,-Strömen alle Ordnungen, die im Drehmo- 
ment messbar waren, kompensieren lassen. Die Kompensation der Be- 
schleunigungen hingegen gelang nicht in allen Betriebspunkten und 
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Richtungen. Grund hierfür sind teilweise die notwendigen sehr hohen 
Amplituden, die vom Regler nicht mehr gestellt werden konnten oder 
auch andere negative vibroakustische Effekte hervorrufen (z.B. Neben- 
ordnungen, siehe Kapitel 6.3.4) 


6.3.4 Verhalten von Stromregelverfahren 


Um das Verhalten verschiedener Stromregelverfahren auf das Aufprä- 
gen der Schwingungen auf die Sollgrößen für die Erstatzströme genauer 
zu untersuchen, wurden drei unterschiedliche Arten der Stromregelver- 
fahren analysiert, die alle auf dem Prinzip der feldorientierten Rege- 
lung beruhen. Bei dieser Untersuchung wird das Übertragungsverhalten 
zwischen Soll- und Istwerten in den Strömen genauer analysiert. 


Für die Analyse wird das Amplitudenverhältnis a der Sollströme mit 
aufgeprägten Schwingungen Is.ı zu den Istströmen Jų bei der Fre- 
quenz der aufgeprägten Ordnung gesucht (ähnlich Kapitel 3.1.3.3 und 
Kapitel 4.4.3). Um das richtige Verhältnis zu bestimmen, wird die Dif- 
ferenz zwischen dem Istwert der aufgeprägten Schwingung Ist und ei- 
ner vorherigen Referenzmessung [pes ohne aufgeprägter Schwingung 
genutzt. 


Tet -1 

g= | Ist Ref | (6.14) 
soul 

6 = Alıst — ref) — (son) (6.15) 


Grundsätzlich zeigen sich hierbei keine erkennbaren Abweichungen im 
Amplitudenverhältnis und der Phasenverschiebung durch die Variation 
der aufgeprägten Amplitude. Dies zeigt beispielhaft das Ergebnis aus 
Abbildung 6.17. 


134 


6.3 Anwendung der Methodik am Prüfstand 


oe 

gs JA 

ER 2A 

E idie 

<E 9 | | | | 4A 
> 50 100 150 200 0-5 A 

2 oF HA 

> Bp —-7A 
5 

32 -100 F 8A 
$= = 

= E 9A 


-200 fi i fi LS 
50 100 150 200 
Ordnung 


Abbildung 6.17: Übertragungsverhalten PI-Regler 


Für weitere Auswertungen wird ein Durchschnittswert aus den 
Amplitudenverhältnissen und Phasenlagen bei verschieden aufgepräg- 
ten Amplituden ermittelt. Folgende Stromregelverfahren wurden dabei 
untersucht: 


e Konventionelles PI-Stromregelverfahren 
Ein häufig verwendetes Stromregelverfahren ist das Verwenden 
eines PI-Reglers. Es besteht aus jeweils einem PI-Regler für die 
geregelten Ströme mit Entkopplung der Regelgrößen und einer 
Vorsteuerung, um der induzierten Spannung entgegen zu wir- 
ken. Die Regelparameter wurden nach dem Betragsoptimum be- 
stimmt. 


e Deadbeat Stromregelverfahren mit linearen Parametern 
[Nus10] 
Um ein besseres dynamisches Verhalten zu erzielen, besteht die 
Möglichkeit, den Regler aus dem Zustandsraummodell der Stre- 
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cke abzuleiten. Hierbei können die Regler auf Deadbeat-Verhal- 
ten eingestellt werden. Dadurch wird theoretisch im nächsten 
Zeitschritt die geänderte Regelgröße eingestellt, was im geteste- 
ten System mit einer Taktfrequenz von fr = 10 kHz einer Zeit 
von Tr = 200 us entspricht. Der integrale Anteil, welcher die 
statische Genauigkeit einregelt, wird mit einer Zeitkonstanten 
von T4 = 40 ms eingestellt. 


Deadbeat Stromregelverfahren mit nachgefiihrten Parame- 
tern [Gem15] 

Moderne elektrische Maschinen besitzen ein oft nichtlineares 
Verhalten, vor allem im Sattigungsbereich. Aus diesem Grund 
können die beiden zuvor dargestellten Stromregelverfahren 
nicht die Ströme in allen Betriebspunkten optimal einregeln. 
Durch den Einsatz eines adaptiven Ansatzes, welcher diese 
Nichtlinearitäten berücksichtigt, kann ein besseres Verhalten 
des Reglers erwartet werden. Die Entkopplung des getesteten 
Verfahrens ist ähnlich dem des PI-Reglers. Allerdings sind die 
Parameter abhängig vom Betriebspunkt des Motors. Da der 
grundsätzliche Aufbau des Reglers ähnlich dem des linearen An- 
satzes ist, werden ähnliche Einstellungen für die Zeitkonstanten 
der Regelung mit T4 = 200 us vorgenommen. 


Abbildung 6.18 zeigt die verschiedenen Strukturen der untersuchten 
Stromregelverfahren. 
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Konventionelles PI-Stromregelverfahren 


Deadbeat-Regler mit linearen Parametern 
ge 
m? 


Deadbeat-Regler mit nachgeführten Parametern 


Abbildung 6.18: Struktur verschiedener Stromregelverfahren nach 
[Gem15] 
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Konventionelles PI-Stromregelverfahren 
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Deadbeat-Regler mit nachgefiihrten Parametern 
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Abbildung 6.19: Vergleich der Regelverfahren 
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Ein Vergleich der Stromregelverfahren in Bezug zu aufgeprägten 
Schwingungen in Abbildung 6.19 zeigt, dass mit dem PI-Regler nur 
schwer sehr hohe Frequenzen aufgeprägt werden können, da das 
Amplitudenverhältnis bei steigender Frequenz, also bei höheren Ord- 
nungen und Drehzahlen, stark abnimmt. Die Deadbeat-Regler hingegen 
können durchaus auch bei sehr hohen Frequenzen Oberschwingungen 
in den Ist-Strömen aufprägen. Somit besteht die Möglichkeit, über den 
gesamten Betriebsbereich Schwingungen aufzuprägen und somit das 
vibroakustische Verhalten zu beeinflussen. 


Ein negativer Effekt, der bei der Aufprägung von Harmonischen in 
Strömen entsteht, sind Nebenordnungen, die sich vor allem bei hohen 
Amplituden zeigen und das vibroakustische Verhalten auch negativ be- 
einflussen können (Abbildung 6.20). 


PI-Stromregelverfahren Regler mit nachgeführten Parametern 


170. Ordnung mit I, (50 Nm 100 rpm) 170. Ordnung mit 1, (50 Nm 100 rpm) 


Differenz I 
a 


o 50 100 150 200 250 300 0 50 100 150 200 250 300 
Ordnung Ordnung 


Differenz 1, [A] 


0 50 100 150 200 250 300 0 50 100 150 200 250 300 
Ordnung Ordnung 


Abbildung 6.20: Entstehung von Nebenordnungen 


In Abbildung 6.20 sind die entstehenden Nebenordnungen in einem Be- 
triebspunkt zu erkennen. Dabei wird die Differenz aus Referenzmes- 
sung und Messung mit aufgeprägter Schwingung gebildet. Die eigent- 
lich aufgeprägte Ordnung, in diesem Fall die 170. Ordnung, nicht dar- 
gestellt, sodass nur die unerwünschten Nebenordnungen, die aufgrund 
des Aufpräge-Verfahrens entstehen dargestellt werden. Erklärt werden 
kann dieses Phänomen, indem die Strangströme betrachtet werden. Da 
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der Regler die aufgeprägten Schwingungen nicht vollständig ausregelt, 
sondern lediglich auf die Abweichungen reagiert, ist es nicht möglich, 
dass dieser den Sollverlauf exakt einstellt. Die Abweichung zwischen 
Ist- und Sollströmen stellt sich im Frequenzbereich vor allem als Ne- 
benordnungen dar. 


24 — Sollverlauf 
— Istverlauf 25 


bk 1 


0 100 200 300 140 160 180 200 
elektrischer Winkel [°] mechanische Ordnung 


Abbildung 6.21: Aufprägung der 170. Ordnung mit 3 A beim Be- 
triebspunkt 100rpm/50Nm mit dem PI-Stromregelverfahren 


Grund für das Entstehen von Nebenordnungen sind fehlende oder ab- 
weichende Informationen über das physikalische Verhalten der Ma- 
schine. So ist zum Beispiel die eigentlich für eine Stromänderung in den 
Strangströmen notwendige Spannung wiederum vom Rotorwinkel und 
dessen Winkeländerung abhängig. Diese Informationen werden im ein- 
fachen Grundwellenregler nicht gesondert rückgeführet und können so- 
mit nicht berücksichtigt werden, wodurch es zu diesen Abweichungen 
kommet. Ebenso können Fehlstellungen und Toleranzen sowie Strom- 
und Winkelsensorfehler zu Nebenordnungen führen. Vergleicht man 
die beiden dargestellten Regler miteinander, so ist zu erkennen, dass bei 
beiden vor allem die Nebenordnungen der + 20. Ordnung im J, Strom 
ausgeprägt sind. Die Ausprägungen im lq Strom sind weniger stark und 
auch weniger kritisch aufgrund der schwach ausgeprägten Übertra- 
gungspfade und werden aus diesem Grund nicht weiter betrachtet. 
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Eine weitere Möglichkeit dieses Phänomen über einen größeren Be- 
reich dazustellen, ist die Isolierung der Amplitude der + 20. Nebenord- 
nung Ay+29 und Ay_20. Um nur eine Größe darzustellen, werden hier 
die arithmetischen Mittelwerte der Amplitude Ay429 und Ay_2o pro 
aufgeprägte Amplitude A, genutzt und der Störfaktor S daraus berech- 


net. 


_ Ans+zo + An-20 


(6.16) 
2: A, 


PI-Stromregelverfahren Regler mit nachgeführten Parametern 


Störfaktor 
Störfaktor 


500 200 
150 
100 


Drehzahl [rpm 0 50 Drehzahl [rpm 0 
P! Ordnung p Ordnung 


Abbildung 6.22: Vergleich des Störfaktors bei verschiedenen Regel- 
verfahren bei 50 Nm 


Die entstehenden Nebenordnungen sind beim Regler mit nachgeführten 
Parametern deutlich geringer als beim konventionellen PI-Stromregel- 
verfahren. Zudem fällt auf, dass bei höheren Drehzahlen bei beiden 
Stromreglern punktuell höhere Nebenordnungen auftreten. 
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6.3.5 Stromoberschwingungen 


Kein anwendbares Stromregelverfahren kann Oberschwingungen in 
den Strömen ganz ausregeln. Es gibt verschiedene Methoden um diese 
zu reduzieren (Kapitel 6: Einleitung), doch aufgrund von Ungenauig- 
keiten in der Modellbildung können diese nie vollständig vermieden 
werden. So zeigen sich in Abbildung 6.23 die entstehenden Stromober- 
schwingungen bei verschiedenen Regelverfahren ohne aufgeprägte 
Schwingungen. 


Spektrum für verschiedene Drehzahlen Spektrum für verschiedene Drehzahlen 
PI-Regler Zustandsregler mit nachgeführten Parametern 


100 600 


400 
20 200 200 
Ordnung Drehzahl [rpm] Ordnung Drehzahl [rpm] 


Abbildung 6.23: Vergleich der Stromoberschwingungen 


Wie schon beschrieben, sind bereits Möglichkeiten zur Reduktion von 
Oberschwingungen in den Strömen bekannt. Allerdings besteht noch 
kein Verfahren oder Methode, um deren Einfluss auf das vibroakusti- 
sche Verhalten experimentell zu ermitteln und somit das Potential für 
eine Verminderung von störenden Schwingungen aufgrund von Ober- 
schwingungen in den Strömen aufzudecken. In den folgenden Absätzen 
wird das während der Arbeit entwickelte Verfahren vorgestellt, um 
diese Potentiale aufzuzeigen und somit auch eine klare Trennung zwi- 
schen den von den Oberschwingungen in den Strömen verursachten 
Vibrationen und den rein aus dem magnetischen Design stammenden 
Vibrationen herzustellen. 
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Hierzu wird eine vereinfachte lineare Modellvorstellung genutzt. 


n LIND] — 
u — cao + y 


r 


Abbildung 6.24: Modellvorstellung zur Trennung der Quellen (I4 und 
Ty) 
q 


Dabei beschreiben G(s) und Gjq(s) die Ubertragungspfade, welche 
die Ist-Ersatzströme J, und l4 durchlaufen, um am Ausgang als Größe 
y gemessen zu werden. Die Störgröße r wirkt in diesem vereinfachten 
Modell direkt auf den Ausgang. Die Sollvorgaben der Ströme werden 
mit verschiedenen Schwingungen überlagert, sodass diese wie in For- 
mel (6.1) und (6.2) beschrieben werden können. 


Die Beschreibung des Systems aus Abbildung 6.24 im Laplace-Bereich 
lautet somit: 


14(5) * Gig(s) + la(s) * Gig(s) + r(s) = y(s) (6.17) 


Betrachtet man dabei eine Kreisfrequenz w = 2r vn, so kann der 
Zusammenhang umgeschrieben werden zu: 


Iq|Giqle*?at 4 + 14|Gqlet?at sera) + A,elPr 


i (6.18) 
= Aye'?y 
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Daraus erhält man ein Gleichungssystem. 


CAUA COS Øq + la |Gial COS Pa + Ar cos p) + 


j . 5 (6.19) 
(Iq |Giq| sin @g + lalGıal sin pa + A, sin Pr) 
lg |Grig| Sin Pg + lalGral sin Pa + A; sin 
Py = tan-? al ial Pq al ial Pa r Pr (6.20) 
IgGl COS Pq + Iq|Giq| COS Pa + A, COS Q, 


Die Unbekannten des Gleichungssystems GAR Gig |Gial, 4Gia, Ar, 
gy, können nach 3 Messungen gelöst werden. Für genauere Ergebnisse 
wird mit Hilfe der kleinsten-Quadrate-Methode die Differenz zwischen 
den erwarteten Werten (Ay), E (yy) und den gemessenen Werten Ay, 


Py minimiert. 


YA, - EA? 


min (6.21) 
1G1ql4G1ql1Gral4G1a-ArPr Yo = E(@y))? 


Abbildung 6.25 zeigt beispielhaft die Aufteilung der verschiedenen Ur- 
sachen auf das vibroakustische Verhalten des Systems mit Hilfe der 
wirkenden Kraft auf den Dämpfer bei verwendeten PI-Stromregelver- 
fahren. Als verwendete Messungen wurden die Kompensationsmessun- 
gen aus 6.3.3 verwendet, welche am Akustikrollenprüfstand mit dem 
Fahrwerksadapter erstellt wurden. Die dargestellten Punkte zeigen die 
Messungen mit aufgeprägten Schwingungen in beiden Steuergrößen, 
also I4 und /,. Die Richtigkeit des Modells wird durch die kleine Norm 
der Residuen sichergestellt. 
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47. Ordnung am Dämpfer 80. Ordnung am Dämpfer 
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Abbildung 6.25: Trennung der Quellen an der Dämpferkraft (J, Ia) 


Hierbei ist zu erkennen, dass kein Potential bei der Verbesserung der 
Stromregelverfahren in der 47. Ordnung zu erkennen ist. Dies ist nur 
logisch, da für diese Ordnung keine Stromoberschwingungen zu erken- 
nen sind. Die Abbildung zeigt weiter die Verstärkungsfaktoren für die 
Gesamtübertragungswege [Gal und |G;q|. Somit stellt die graue Kurve 
das berechnete vibroakustische Verhalten bei einer Ansteuerung ohne 
Stromoberschwingungen und somit das Potential für einen sehr guten 
Regler, wie er beispielsweise in [Karl6, Ricl5c, Pell2b] beschrieben 
ist, dar. Beispielhaft dargestellt sind dazu die Auswirkungen am Dämp- 
fer der 47. und 80. Ordnung. 


Es fällt auf, dass der Verstärkungsfaktor |G,4| der Übertragungswege 
für aufgeprägte Schwingungen vom Ersatzstrom /, sehr gering ist. Aus 
diesem Grund wurde für nachfolgende Untersuchungen ein vereinfach- 
tes Modell angenommen, sodass das vibroakustische Verhalten nur 
noch durch den Ersatzstrom I, beeinflusst werden kann. Der Vorteil bei 
dieser Vereinfachung liegt in einer höheren Auflösung der zur Verfü- 
gung stehenden Messdaten, sodass diese einen wesentlich besseren 
Aufschluss über das Verhalten über den Drehzahlbereich liefern kön- 
nen. 
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+) y 


Abbildung 6.26: Reduziertes Modell 


Somit gilt analog zu der Berechnung aus den Formeln (6.19), (6.20) 
und (6.21): 


(1,|Giq| COS Pq + A, cos p) + 


= (6.22) 
y . f 2 
(1,|Giq| sin q + A, sin Pr) 
I,|Gjq| Sin Pq + A, sin 
Py = tan! (en tiine (6.23) 
14|Giq| COS Pq + A, COS p, 
YA, - Ey)? 
(6.24) 


min 
|Giq|4G1q.Ar.Pr Yo an E(p,))? 


Die Ergebnisse der berechneten Unbekannten zeigen sich in Abbildung 
6.27. Die Abweichungen des Modells werden als die Norm der dadurch 
entstehenden Residuen dargestellt und zeigen eine gute Übereinstim- 
mung mit der Modellvorstellung. 
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47. Ordnung 80. Ordnung 
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Abbildung 6.27: Trennung der Quellen an der Dämpferkraft (/,) 


Aufgrund der vorliegenden Messdaten konnten hier die Intervalle der 
Geschwindigkeiten verkleinert werden, sodass ein genauerer Verlauf 
dargestellt wird. Ebenso kann aus dem Diagramm sehr gut abgeschätzt 
werden, wie hoch die Amplituden zur Kompensation der Wechselkräfte 
sein müssen. Zudem können aus der Abbildung die Verstärkungsfakto- 
ren IG ‘al der Ubertragungswege des [,-Stroms zu den gemessenen Stel- 
len abgelesen werden. Darauf sind als Überhöhungen die verschiedenen 
mechanischen Resonanzen zu erkennen. Weitere Messungen an ande- 
ren Anbindungspunkten und mit geänderten Ordnungen finden sich im 
Anhang (Abbildung 7.9 - Abbildung 7.12). 


Als besonders kritisch in Bezug auf Stromoberschwingungen und deren 
Auswirkungen auf Vibrationen und Drehmoment wird in der Literatur 
vor allem die 6. elektrische Ordnung beschreiben, was im Falle des 
Radnabenantriebs aufgrund seiner 20 Polpaare der 120. Ordnung ent- 
spricht. Grundsätzlich zeigt sich auch hier, dass durch diese Ober- 
schwingungen das vibroakustische Verhalten negativ beeinflusst wird. 
Allerdings muss dies nicht über den gesamten Betriebsbereich der Fall 
sein, wie sich vor allem zwischen den Geschwindigkeiten 15 — 18 km/h 
zeigt. Untersuchungen haben gezeigt, dass auch bei anderen Ordnungen 
dies positive Effekte haben kann. 
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120. Ordnung 
0.4 
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Abbildung 6.28: Trennung der 120. Ordnung 


Mit Hilfe des in der Arbeit vorgestellten Steuerungsansatzes und dem 
dazugehörigen Parameterridentifikationsalgorithmus ist es ebenfalls 
möglich, sich auf die Stromoberschwingungen in den Ersatzströmen zu 
beschränken. Dies würde das aufgezeigte Potential (orangene Fläche) 
in Bezug auf die vibroakustischen Eigenschaften nutzen und hätte zu- 
dem den Vorteil, dass hier auch eine aktive Rückkopplung ohne zusätz- 
lichen Sensoraufwand möglich wäre, da Sensoren für Strangstrommes- 
sung bei der Regelung elektrischer Maschinen notwendig sind. 
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Geschätze Fläche nach der 5. Messung 120. Ordnung Vergleich der 5. Messung mit der Referenzmessung 
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Ordnung 


Beschleunigung [m/s?] 


es ü . 100 150 200 250 300 
Phasenlage [°] 0 Amplitude [A] Ordnung 


Abbildung 6.29: Kompensation I, bei 50 rpm 


Abbildung 6.29 zeigt die Kompensation der harmonischen Ober- 
schwingungen des /,-Stroms und dessen Auswirkungen auf das Dreh- 
moment. So wird die Ordnung im Strom fast vollständig kompensiert 
und im Drehmoment reduziert sich die Welligkeit ebenfalls sehr stark. 
Die teilweise noch vorhandene Welligkeit des Drehmomentes ist rein 
auf das elektromagnetische Design zurückzuführen. 
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6.3.6 Kombinierte Kompensation 


Mit der dargestellten Methodik ist es möglich, auch mehrere Ordnun- 
gen gleichzeitig zu kompensieren, wie Abbildung 6.30 zeigt. Als Prü- 
feinrichtung wurde der Motorenprüfstand des ETI mit Drehmoment- 
messwelle verwendet. Die folgenden Messungen wurden mit dem in 
Kapitel 6.3.4 beschriebenen PI-Regler durchgeführt. 


T T 
-© Referenz 
—©—aufgeprägte Schwingung |- 


= 
w 


Drehmoment [Nm] 
o o 
h N 


0 20 40 60 80 100 120 140 160 
Ordnung 


Abbildung 6.30: Kompensation des Drehmomentes bei 60 rpm mehre- 
rer Ordnungen 


Wie in Abbildung 6.30 zu erkennen ist, wurden an diesem Betriebs- 
punkt mehrere Ordnungen (40., 47. und 120. Ordnung) gleichzeitig 
kompensiert. Dazu wurden die aus der bereits vorgestellten Ermitt- 
lungsmethodik berechneten Parameter verwendet, welche jedoch bei 
einzeln aufgeprägten Ordnungen ermittelt wurden. Grundsätzlich zeigt 
sich eine gute Kompensation der Drehmomentungleichförmigkeiten. 
Allerdings sind auch zusätzlich eingebrachte Störungen, wie etwa die 
60. Ordnung, zu erkennen. Das Übertragungsverhalten der Regelung 
auf einzelne harmonische Oberschwingungen auf den Sollstrom wurde 
bereits in Kapitel 6.3.4 angesprochen. Werden jedoch mehrere Schwin- 
gungen überlagert aufgeprägt, kann je nach Phasenlage zueinander das 
Übertragungsverhalten der Regelung stark oder weniger stark variieren 
und somit auch zusätzliche Störungen eingebracht werden. Dies zeigt 
sich, wenn man den Zeitverlauf der aufgetragenen Sollverläufe im ro- 
tierenden und dreiphasigen statorfesten Koordinatensystem betrachtet 
(Abbildung 6.31). 
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Abbildung 6.31: Ströme zur Kompensation des Drehmoments 


Die zeitliche Stromänderung in den Strangströmen kann stark zu- oder 
abnehmen, je nach Überlagerung der verschiedenen Ordnungen. Für 
eine gleiche Umsetzung der Amplitude und Phasenlage folgt somit ein 
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erhöhter Spannungsbedarf im Vergleich zu der Kompensation mit nur 
einzelnen Ordnungen. Diese Änderung kann mit einem PI-Regler nur 
bei niedrigen Frequenzen der aufgeprägten Schwingungen kompensiert 
werden. Wie bereits beschrieben können mit einer Erweiterung der Re- 
gelung bessere Ergebnisse erzielt werden. 
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6.4 _Betriebsstrategie 


In den vorherigen Kapiteln wurde beschrieben, wie die Parameter für 
die Kompensation gefunden werden können. Diese stellen zunächst die 
vibroakustisch optimalen Parameter für die Sollvorgaben (Amplitude 
und Phasenlage) des PI-Reglers dar. Sie werden in diesem Kapitel ver- 
wendet, um an einem Hochlauf zu demonstrieren, welche Herausforde- 
rungen diese Methodik mit sich bringt und welche Möglichkeiten sie 
bietet. Es wird hier vor allem auf die Kompensation einzelner Ordnun- 
gen über die Drehzahl eingegangen. 


6.4.1 Abhängigkeiten 


Grundsätzlich wurden Effekte, welche vor allem die Stromregelung mit 
sich bringt, in Kapitel 6.3.4 angesprochen. Im Folgenden soll kurz auf 
die Abhängigkeiten eingegangen werden, die für die Kompensations- 
parameter wichtig sind. Im Rahmen des Forschungsvorhabens konnten 
nicht alle angesprochenen Abhängigkeiten, welche die Parameter be- 
einflussen, im Experiment nachgebildet werden. Es ist jedoch klar, dass 
bei einer Änderung der mechanischen Übertragungswege die Parameter 
für die aufgeprägten Schwingungen beeinflusst werden können, sofern 
bei einer Kompensation der Ursache nicht direkt im Motor entgegenge- 
wirkt wird. So ist es zum Beispiel möglich Schwingungen, die aus 
Drehmomentungleichförmigkeiten resultieren, auch wenn sie den Ur- 
sprung im magnetischen Design haben, durch aufgeprägte Schwingun- 
gen im /,-Strom zu kompensieren. Schwingungen, die andere Ursachen 
haben, wie beispielsweise radiale Kräfte auf Rotor und Stator, werden 
vor allem durch das Aufbringen zusätzlicher Drehmomentschwingun- 
gen im späteren Anbindungspunkt kompensiert. Hierbei können Ände- 
rungen im Übertragungsverhalten der mechanischen Komponenten 
eine große Rolle spielen. Gründe für diese Veränderungen sind vor al- 
lem eine Änderung der Steifigkeiten, welche verursacht werden können 
durch: 
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e Alterungsprozesse (z. B. Längslenkerbuchse) 
e Temperaturabhängigkeiten 


Anderung der Beladung 


Luftdruckänderungen im Reifen 
e Abhängigkeiten des Zustandes im Dämpfer (gleitend oder ste- 
hend) 


6.4.2 Betrachtung der Verluste 


Grundsätzlich treten in permanentmagneterregten Synchronmaschinen, 
wie dem Radnabenanrieb, folgende Verluste auf: 


e Mechanische Verluste 
o Reibung Radlager 
o Reibung durch Strömung am Rotor 
e Elektrische Verluste 
o  Kupferverluste (Ohmsche Verluste in den Stator- 
wicklungen) 
o  Eisenverluste 
=  Hystereseverluste 
= Wirbelstromverluste 
= Excessverluste 


Mechanische Verluste spielen beim Vergleich der akustisch optimier- 
ten Ansteuerung im Gegensatz zu einer Ansteuerung ohne aufgeprägte 
Schwingungen keine Rolle. Eine größere Rolle hingegen könnten die 
elektrischen Verluste haben. Da diese jedoch nicht Fokus der Arbeit 
sind, wird nur stichpunktartig auf die Kupferverluste eingegangen. 
Diese können durch folgenden Ausdruck beschrieben werden. 


3 
P= a Rs: (Ig? + Iq”) (6.25) 


154 


6.4 Betriebsstrategie 


Bereits an dieser Formel ist zu erkennen, dass sich bei Oberschwingun- 
gen in den Strangströmen die Verluste aufgrund des quadratischen Zu- 
sammenhangs im Gegensatz zu Strömen ohne Oberschwingungen the- 
oretisch vergrößern. In der Praxis treten jedoch auch ohne gewollt auf- 
geprägte Schwingungen Oberschwingungen in den Strömen auf. Eine 
akustisch optimierte Ansteuerung mittels HCI kann somit bei kleinen 
Amplituden auch dazu führen, dass die quadratische Summe der Ströme 
- (Ia? + Iq”) im Mittel nicht zunehmen muss. Stichprobenartige Mes- 
sungen am Beispiel des Radnabenantriebs zeigen, dass bei aufgepräg- 
ten Schwingungen diese Verluste nicht signifikant zunehmen müssen 
beziehungsweise auf gleichem Niveau liegen. Dazu beispielhaft aufge- 
führt ist in Abbildung 6.32 die quadratische Summe der Ströme, um den 
Unterschied bei aufgeprägten Schwingungen aus dem Beispiel von Ab- 
bildung 6.30 in den ohmschen Verlusten abschätzen zu können. Hierbei 
ist kein messbarer Unterschied festzustellen. 


500 p 
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N 400 | i 
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as | ! Ai | Mi 
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mechanischer Winkel [°] 


Abbildung 6.32: Quadratische Summe von J, und Iq bei aufgeprägten 
Schwingungen 


Hystereseverluste entstehen im Allgemeinen bei der Umpolung der 
weißschen Bezirke in einem ferromagnetischen Material. Um eine tat- 
sächliche Abschätzung treffen zu können, müsste das dafür verantwort- 
liche magnetische Feld in einer FEM Simulation berechnet und analy- 
siert werden. Gleiches gilt für die Wirbelstrom und Excessverluste. 
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Eine einfache Modellierung dieser Effekte ist nicht möglich. Aufgrund 
der für die Kompensation notwendigen geringen Amplituden ist jedoch 
davon auszugehen, dass die durch die aufgeprägten Schwingungen ver- 
ursachten Verluste sehr gering und diese somit dort zu vernachlässigen 
sind. 


6.4.3 Anwendung 


Für die in der Arbeit beschriebene Strategie wird vorerst nur die Kom- 
pensation einer bestimmten Ordnung beim Hochlauf der Drehzahl ge- 
testet. Die Parameter stammen aus den gezeigten Kompensationsmes- 
sungen in Kapitel 6.3.3. Zwischen den Sollvorgaben wird linear inter- 
poliert. Als Prüfeinrichtung wurde der Fahrwerksadapter für das 
Radnabensystem mit der flexiblen Ansteuereinheit am AARPS verwen- 
det. 


Abbildung 6.33 zeigt beispielhaft den Vorgang zur Ermittlung und 
Kompensation der 120. Ordnung am Beschleunigungssensor der 
Längslenkerbuchse in x-Richtung. Als Regelung für die Ströme wurde 
ein PI- Regler verwendet. 


Deutlich erkennbar ist die Kompensation der 120. Ordnung bei höheren 
Drehzahlen im Beschleunigungssignal. Ebenso führen folgende Phäno- 
mene zu einer nicht kompletten Kompensation der störenden Ordnung: 


e Schlechte Interpolation der Parameter zwischen den angepassten 
Parametern aufgrund der Stützstellen 

e Dynamische Vorgänge im Regler beim Durchlaufen der Be- 
triebspunkte 

e Störungen im Winkelmesssignal aufgrund elektromagnetischer 
Einflüsse 
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1. Ermittlung des Kompensationsbereiches 


120. Ordnung 
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kompensierender 
Schwingungen: 
38-53 km/h 
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Abbildung 6.33: Kompensationsvorgang Beschleunigung LLB (x) 
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6 Akustisch optimierter Steuerungsansatz 


6.5 Vergleich verschiedener 
Aufprägemethoden 


Grundsätzlich kann die Methode zur Suche der Parametereinstellung 
(Amplitude und Phasenlage) auch für andere Methoden zur Aufprägung 
von Schwingungen genutzt werden, wie sich bereits in anderen Unter- 
suchungen gezeigt hat. Dazu wurden in dieser Arbeit zwei weitere 
Möglichkeiten und deren Auswirkungen auf die Ströme in einem Ver- 
such am Fraunhofer IISB genauer untersucht. Hierbei wurde die neue 
Version des Radnabenantriebs Gamma genutzt, der an einer Abtriebs- 
maschine angebracht war, sodass ein Gegenmoment erzeugt werden 
konnte. Für die Untersuchung wurde eine feste Drehzahl von 300 rpm 
mit verschiedenen Lasteinstellungen (200 und 400 Nm) getestet. 


6.5.1 Aufgeprägte Schwingung in den 
Spannungen 


Eine andere Möglichkeit, Schwingungen in den Strömen aufzuprägen, 
ist die Addition der aufgeprägten Schwingung in Form einer überlager- 
ten Spannung auf die bereits berechnete Stellgröße aus der Regelung, 
wie sie bereits in [Oks16, Qil5] verwendet wurde. Diese ist nachfol- 


gend dargestellt. 
8 > 
x 
v 
Tu lw ly > 
| | | A + > sin 
— 
>| 
A x | 
q > 
>| 
+ Sellvoreabe Stellgröße aus Regelung ne 
ar Y Ud 
Ua 
L 
Feldorientierte Pa + >) sin 
Regelung >| 
Aq > xX 
| 
[| * Solvorgabe Stellgröße aus Regelung > + —, 
q u 
Ug q 


Abbildung 6.34: Aufgeprägte Schwingungen in der Spannung 
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Allerdings bringt diese Methode, im Gegensatz zu der bereits beschrie- 
ben, auch Nachteile mit sich: 


e Der Einfluss auf das Magnetfeld und somit auch auf die aufge- 
brachten Kräfte entsteht durch eine geänderte Bestromung. So- 
mit ist der Strom die Kenngröße, die den Einfluss unabhängiger 
darstellen kann, auch wenn die gestellte Spannung die Ursache 
der Stromänderung ist. Die Kenngröße Spannung ist also weni- 
ger aussagekräftig. 

e Der Regler versucht, die durch die aufgeprägte Spannung ent- 
stehenden Ströme, wieder auszuregeln. Eine Verbesserung des 
Störverhaltens des Reglers hat somit auch eine Erhöhung der 
Amplitude für die aufgeprägte Spannung zufolge. 

e Starke Abhängigkeit von Drehzahl 


Um den Einfluss der aufgeprägten Spannung zu visualisieren, wurde 
eine ähnliche Berechnungs- und Darstellungsweise wie in den vorheri- 
gen Kapiteln genutzt (Abbildung 6.35). Die orange Fläche zeigt den 
Referenzwert ohne aufgeprägte Schwingungen der 120. mechanischen 
beziehungsweise 6. elektrischen Ordnung. Durch das Aufprägen genau 
dieser Ordnung mit verschiedenen Amplituden und Phasenlagen erge- 
ben sich Bereiche, in denen man eine konstruktive oder destruktive 
Überlagerung erreichen kann. Das System verhält sich linear und die 
beschriebene Methodik zur Parametersuche, auch zur Kompensation 
von mechanischen Schwingungen, könnte somit genutzt werden. Das 
in Abbildung 6.35 gezeigte Spannungsamplitudenverhältnis beschreibt 
lediglich die hier verwendete Spannungsamplitude A, multipliziert mit 
der Drehzahl. 
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120. Ordnung im I, Strom 120. Ordnung im 1, Strom 
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Abbildung 6.35: Einfluss der aufgeprägten Spannungsamplitude auf 
die Ströme (links: 200 Nm, rechts: 400 Nm) 
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6.5 Vergleich verschiedener Aufprägemethoden 


6.5.2 Aufgepragte Ströme in den Phasen U, V, W 


Ebenso ist es möglich, aufgeprägte Schwingungen in Strömen durch 
eine zusätzlich überlagerte feldorientierte Regelung zu erzeugen (ähn- 
lich [Pell2b]). Als Eingang für die überlagerte feldorientierte Regelung 
dient die Differenz der Phasenströme zur perfekten Grundwelle. Diese 
Differenz wird über den mit der Ordnung v multiplizierten Rotorwinkel 
vom stehende ins drehende Koordinatensystem umgewandelt, sodass 
die Phasenlage und Amplitude der aufgeprägten Schwingung als 
Gleichgröße beschrieben und ausgeregelt werden kann. 
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Abbildung 6.36: Aufgeprägte Schwingung in Strom durch tiberlagerte 
feldorientierte Regelung 


Hiermit können vor allem die Strangströme direkt mit der entsprechen- 
den Ordnung v beeinflusst werden. Der Zusammenhang zwischen der 
Aufprägung im statorfesten Koordinatensystem und dem drehenden 
Koordinatensystem wird nachfolgend beschrieben. 
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6 Akustisch optimierter Steuerungsansatz 


Dabei ist die aufgeprägte Schwingung in den Phasenströmen mit der 
Amplitude A, und der Phasenlage p, beschrieben und es gilt für diese 
Ströme: 


Iy = lysin + lw = lusin + Av * sin(v@ + Py) 
, 21 
ly = lysin + Iw = lysin + Ay * sin (ve + Qo, - =) (6.26) 


Art 
Iw = Iwsin + bw = Iwsin + 4 ' sin (vo + Py — =) 


Betrachtet man Iq mit 


la = lasruna + laober (6.27) 


und 


2 21 
lagruna = 3 [st "cos + Iysin * COS (o = —) + Iysin 


cfo- 


ergibt sich daraus 


ldober = Ay* sin(8 (v m 1) + Py) (6.28) 


dementsprechend gilt für /, mit Ig = Iggruna + Iqover 


Igober =A,' cos(0 (v —1)+ Py) (6.29) 
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6.5 Vergleich verschiedener Aufprägemethoden 


Die Aufprägung hat also vor allem Auswirkungen auf die Ströme im 
drehenden Koordinatensystem mit der Ordnung v — 1 . 


Nachteil dieser Methode ist nur die Möglichkeit der direkten Beeinflus- 
sung der Ordnungen im statorfesten Koordinatensystem, welche als 
Folge eine Beeinflussung der Ströme im drehenden Koordinatensystem 
mit der Ordnung v — 1 haben. 


Um hier eine nur in q- oder d-Richtung wirkende Oberschwingung auf- 
zuprägen, wären zwei Regler für die jeweiligen Ordnungen notwendig, 
wie sich bereits in Gleichung (6.7) gezeigt hat. 


Auch hier wird die Möglichkeit zur Aufprägung von Schwingungen mit 
der vorher gezeigten Darstellungsmöglichkeit angewandt (Abbildung 
6.37), um die Linearität des Verfahrens und die Möglichkeit der Funk- 
tion zur Oberschwingungsreduzierung zu zeigen. In orange ist die Re- 
ferenz ohne aufgeprägte Schwingungen gezeigt, in grün die Variation 
der Phase und Amplitude. In violett zeigt sich der Wert bei ausgeschal- 
teter Oberschwingungs-Regelung. 
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Abbildung 6.37: Einfluss der überlagerten Regelung auf die Ströme 
(links: 200 Nm, rechts: 400 Nm) 
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6 Akustisch optimierter Steuerungsansatz 


6.6 Messung am Gesamtfahrzeug 


In Messungen am Gesamtfahrzeug werden die gezeigten Zusammen- 
hänge und vor allem die Wirksamkeit eines vibroakustisch optimierten 
Steuerungsansatzes gezeigt. Hierfür wurde der Akustik-Rollenprüf- 
stand der Schaeffler Engineering GmbH in Werdohl genutzt. Als Ver- 
suchsträger kam ein Ford Fiesta zum Einsatz, auf dessen Hinterachse 
zwei Radnabenantriebe verbaut sind. Im Gegensatz zu der auf den vor- 
herigen Kapiteln vorgestellten Version des Radnabenantriebes „beta“ 
musste hier auf die aktuelle Version „gamma“ zurückgegriffen werden. 
Grund hierfür war dessen Verfügbarkeit im Fahrzeug. Die wesentlichen 
Unterschiede zur Version beta sind die geänderte Statortopologie (60 
Zähne) und die geänderte Rotortopologie (gestaffelter Rotor). 


Grundsätzlich ergeben sich aufgrund der geänderten Motortopologie 
und Konstruktion andere Ordnungen und Resonanzen. Die Messungen 
am Gesamtfahrzeug sollen jedoch auch zeigen, dass die Kompensati- 
onsmethode über die Ansteuerung der elektrischen Maschine generell 
geeignet ist, um akustischen Problemen im elektrifizierten Antriebs- 
strang zu entgegnen. 


Abbildung 6.38: Explosionsdarstellung Radnabenantrieb „gamma“ 
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Allgemeine Ordnung Hauptsächliche 
Berechnung EWD gamma Ursprungskraft 


radiale Kräfte auf Stator 
2p:i 40-i > asym. ÜW des Stators 
> resultierende Axialkraft 
radiale Kräfte auf Rotor 
N'i 60-i > resultierende Axial- 
kräfte 
radiale Kräfte auf Rotor 
N-i+1 60-i+1 > asym. ÜW des Rotors 
+ Rotoreffekt 
tangentiale Kräfte 
(Lastpulsation) 
> Übertragungspfad Sta- 
tor stark ausgeprägt 


kgV(m, 2p) -i kgV (60,40) - i 


tangentiale Krafte 
(Nutrastmoment) 

> Ubertragungspfad Sta- 
tor stark ausgeprägt 


Tabelle 6.1: Entstehende Ordnungen beim EWD „gamma“ 


kgV(N,2p)-i  kgV(60,40)-i 


Es ist erkennbar, dass bei dieser Topologie die erste Harmonische der 
Lastpulsation und des Nutrastmoment auf die 120. Ordnung fällt, was 
bei diesem System auch die hautsächlich störende Anregungsordnung 
ist. Ein weiterer Anregungsmechanismus in diesem System sind Axial- 
kräfte, welche aufgrund der nicht symmetrischen Staffelung des Rotors 
entstehen. So entstehen an den Rändern des Rotors sowie zwischen den 
Staffelungen Kräfte mit der gleichen Ordnungszahl der Drehmomen- 
tungleichförmigkeiten. Ebenso können diese aufgeteilt werden in 
stromabhängige und nicht stromabhängige Effekte. 
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6.6.1 Versuchsaufbau 


Im Fahrzeuginnenraum wurden auf Kopfhöhe drei Mikrofone ange- 
bracht, um die vom Fahrer wahrnehmbaren Geräusche messen zu kön- 
nen. Da bei diesem Rollenprüfstand die linke und rechte Rolle an der 
Hinterachse gekoppelt sind, wurde eine Vorrichtung mit einem Luftbal- 
ken angefertigt (siehe Abbildung 7.4), um Messungen beim Fahrzeug 
mit nur einem drehenden Radnabenantrieb vornehmen zu können. So- 
mit ist es möglich, im Fahrzeuginnenraum die Herkunft der akustischen 
Auffälligkeiten dem jeweiligen Radnabenantrieb zuzuordnen. 


Abbildung 6.39: Akustischer Rollenprüfstand der Schaeffler Enginee- 
ring GmbH 


Eine genauere Beschreibung des Versuchsaufbaus ist in Kapitel A.3 zu 
finden. 
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6.6.2  Versuchsdurchführung und Ergebnisse 


6.6.2.1 Kompensation Radnabenantrieb einzeln 


Bei den Versuchen konnte ähnlich zu den Versuchen am AARP vorge- 
gangen werden. Aufgrund des Aufbaus und der fehlenden direkten 
Rückkopplungsmöglichkeit in der Ansteuerungseinheit für den Radna- 
benantrieb im Gesamtfahrzeug wurde hier ein leicht abgeänderter Ab- 
lauf zur Bestimmung der akustisch optimalen Amplitude und Phasen- 
lage verwendet. So wurde kein iteratives Verfahren verwendet, welches 
sich an die optimale Phasenlage und Amplitude annähert, sondern es 
wurden vorab Amplituden für eine Kompensation geschätzt und diese 
dann mit um 120° versetzten Phasenlagen nacheinander bei einer Dreh- 
zahlrampe des zu untersuchenden Bereichs aufgebracht. Extrahiert man 
nun aus diesen Messungen die zu untersuchende Ordnung, können aus 
den Messungen die vibroakustisch kompensierend wirkenden Parame- 
terwerte (Amplitude und Phasenlage) der aufgeprägten Schwingung 
mit Hilfe der Gleichungen (6.9) bis (6.15) berechnet werden. Voraus- 
setzung ist hier die Linearität des Systems, wobei hier aufgrund der vo- 
rausgegangen Versuche und Simulationen keine starken Nichtlinearitä- 
ten zu erwarten sind. Abbildung 6.40 zeigt den Prüfablauf für diese 
Vorgehensweise. 
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Abbildung 6.40: Ablauf Parameterberechnung im Gesamtfahrzeug 


Vorteile des verwendeten Ablaufs: 


e Durch Drehzahlrampe (alternativ Momentrampe) Berechnung 
der akustisch optimierten Parameter fiir alle Drehzahlen 

e Wesentlich schnellere Möglichkeit für das System die Einstell- 
parameter zu finden 

e Einfache Möglichkeit zur Adaption der Methode bei Systemen, 
die nicht ohne Weiteres mit einer rückführenden Größe arbeiten 
können 


Nachteile des verwendeten Ablaufs: 


e Keine rückführende Größe und somit keine direkte Überprü- 
fung, ob berechnete Amplitude und Phasenlage richtig ist, was 
bei Nichtlinearitäten, Störungen oder Fehlmessungen auftreten 
kann 

e Keine konstante Drehzahl für die Fensterung, wodurch Unge- 
nauigkeiten bei der Berechnung entstehen können 
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Für diese Untersuchung wurde vor allem die 120. Ordnung und die 40. 
Ordnung untersucht. Bei dieser Motortopologie ist die 120. Ordnung 
dem Drehmomentrippel und die 40. Ordnung den Radialkräften zuzu- 
ordnen (siehe Tabelle 6.1). 


Geschätze Fläche bei 40 rpm Geschätze Fläche bei 180 rpm 
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Abbildung 6.41: Veranschaulichung der Gleichung 


So werden für unterschiedliche Drehzahlen die Amplituden und Pha- 
senlagen der einzelnen Ordnungen dargestellt. Mit dem Bestimmtheits- 
maß R? wird ein Wert bestimmt, welcher eine Aussage über die Güte 
des Modells trifft. Es berechnet sich aus 


R2=1 ee 5 unerklärte Variation 
u vn -M Gesamtvariation (6.30) 
mit 


= 


yı=yte 


wobei n die Anzahl der Beobachtungen, e; der Fehler zwischen Mo- 
dellwert y; und Messung ¥, und y der Mittelwert der Messungen ist. 


y= 


ale 


D 


n 
i=1 
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Abbildung 6.42: Berechnung der Parameter tiber Drehzahl 


Abbildung 6.42 zeigt folgende Systemeigenschaften fiir eine vibroakus- 
tisch optimierte Ansteuerung: 


Schlechtes GiitemaB vereinzelt bei 24, 56, 148 rpm 


Grund hierfür ist einerseits die große Änderung der Amplitude und Pha- 
senlage über die Drehzahl (24, 56 rpm), wodurch Ungenauigkeiten bei 
der Synchronisation der Messungen zu diesem Effekt geführt haben 
kann. Drehzahl 148 rpm hingegen deutet auf eine hohe notwendige 
Amplitude zur Kompensation hin, wobei Nichtlinearitäten im System 
eine größere Rolle spielen würden. Aufgrund des niedrigen Schall- 
druckpegels in diesem Bereich kann dieser Eintrag jedoch vernachläs- 
sigt werden. 


Stark variierende Amplituden- und Phasenlagen im Bereich 24, 56, 
95, 115,125, 145, 160 rpm 


Grundsätzlich ist davon auszugehen, dass über die Drehzahl hinweg 
hauptsächlich direkt das Drehmomentrippel kompensiert wird. Dieser 
Bereich hingegen deutet darauf hin, dass auch andere Anregungsme- 
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chanismen in ihrer Resonanz in den Vordergrund treten können. Auf- 
geprägte Schwingungen in den Strömen führen hier zwar auch zu einer 
Kompensation der Schwingung an der Messstelle, allerdings nicht 
durch direkte Kompensation des ursprünglich störenden Kraft oder des 
Drehmomentrippels, sondern durch destruktive Überlagerung verschie- 
dener Anregungsmechanismen. Weitere Untersuchungen haben ge- 
zeigt, dass hier der Anregungsmechanismus Axialkräfte sind, welche 
durch die nicht symmetrische Staffelung beziehungsweise Schrägung 
der Rotormagnete auftreten. 


Aus den berechneten Parametern kann eine nutzbare Kennlinie erzeugt 
werden, indem die Werte unter einem bestimmten Gütemaß nicht be- 
rücksichtigt werden. 
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Abbildung 6.43: Veranschaulichung Punktauswahl über Gütemaß 


Dieser Ablauf kann für alle Momente wiederholt werden, um ein Kenn- 
feld zu generieren, welches für das jeweilige Rad nutzbar ist. In Abbil- 
dung 6.44 sind diese für den linken und den rechten Radnabenantrieb 
dargestellt. Es sind in der Amplitude leichte Unterschiede zwischen den 
Radnabenantrieben erkennbar. Gründe hierfür können leicht unter- 
schiedliche Übertragungswege sowie eine unterschiedliche Anregung 
in den Radnabenantrieben sein. 
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6 Akustisch optimierter Steuerungsansatz 
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Abbildung 6.44: Kennfelder für linken E-WD (links) und rechten E- 
WD (rechts) 


Die gezeigte Vorgehensweise konnte neben der 120. Ordnung auch er- 
folgreich auf die 40. Ordnung im Gesamtfahrzeug angewandt werden, 
wie Abbildung 6.46 zeigt. 
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Abbildung 6.45: Schalldruck rechter E-WD ohne (links) und mit 
(rechts) HCI bei 400 Nm — 120. Ordnung 
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6.6 Messung am Gesamtfahrzeug 
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Abbildung 6.46: Schalldruck rechter E-WD mit/ohne HCI bei 400 rpm 
— 40. Ordnung 


6.6.2.2 Kompensation gesamt 


Mit Hilfe der entstandenen Kennfelder für den linken und den rechten 
Antrieb kann nun die Kompensation für den Schalldruck im Fahrzeug- 
innenraum bestimmt werden. Aufgrund von subjektiven Eindrücken, 
die während der Testreihen im Fahrzeug entstanden, wurden in den 
Versuchen am Gesamtfahrzeug mit beiden Antrieben ausschließlich die 
120. Ordnung kompensiert, da diese den wesentlichen Beitrag für 
vibroakustische Störungen liefert. Die 40. Ordnung hingegen war sub- 
jektiv kaum wahrnehmbar. Die Messungen mit nur einem Antrieb zei- 
gen jedoch, dass auch diese Ordnung durch aufgeprägte Schwingungen 
und der hier vorgestellten Methode sehr gut kompensierbar ist. Die 
Wirksamkeit der Methode konnte somit auch im Gesamtfahrzeug für 
das System Radnabenantrieb „gamma“ gezeigt werden. Aufgrund des 
notwendigen Implementierungsaufwandes in der integrierten Steuer- 
einheit für aufgeprägte Schwingungen konnten nicht mehrere Ord- 
nungskompensationen in einem Betriebspunkt gleichzeitig durchge- 
führt werden. 
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6 Akustisch optimierter Steuerungsansatz 
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Vergleich der 120. Ordnung bei Drehzahlrampe mit und ohne HCI bei 400 Nm 
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Abbildung 6.47: Vergleich mit beiden EWDs 
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Abbildung 6.48: Vergleich mit beiden EWDs ohne (links) und mit 


174 


(rechts) HCI 


7 Zusammenfassung und 
Ausblick 


In dieser Arbeit wurden die grundlegenden Entstehungsmechanismen 
der vibroakustischen Störungen aufgezeigt und mit Hilfe eines multi- 
physikalischen Modells am Beispiel des Radnabenantriebs nachgebil- 
det. Darauf aufbauend wurden mögliche Gegenmaßnahmen abgeleitet 
und ein vibroakustisch optimierter Steuerungsansatz beschrieben, un- 
tersucht und an einem speziell für diese Arbeit entwickelten Prüfaufbau 
sowie am Gesamtfahrzeug angewandt. Mit der im Zentrum der Arbeit 
stehenden neu entwickelten Suchmethodik zur Parameterfindung kann 
dieser Ansatz zur Kompensation elektromagnetischer Anregung in 
elektrischen Maschinen mit nur wenigen Messungen auf andere elekt- 
rische Maschinen übertragen werden. 


7.1 Zusammenfassung 


Die in der Arbeit beschriebenen Inhalte und Zusammenhänge können 
in zwei Themengebiete zusammengefasst werden. 


Thema: Modellbildung und Simulation 


Auf Basis einer 2D-FEM Simulation des magnetischen Feldes wurden 
die wirkenden Kräfte im Luftspalt genauer untersucht. Ebenso konnten 
die Übertragungspfade der Oberschwingungen im Strom zu den wir- 
kenden Kräften beschrieben werden. Mit Hilfe dieser grundlegenden 
Untersuchungen wurden bekannte Simulationsmethoden miteinander 
verglichen und daraus eine neue Möglichkeit zur schnellen Berechnung 
der magnetischen Kräfte mit Berücksichtigung von Stromoberschwin- 
gungen entwickelt (Kraftrekonstruktion). 
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7 Zusammenfassung und Ausblick 


Ebenso wurden mechanische FEM-Modelle der verschiedenen Kompo- 
nenten (Fahrwerk, Stator und Rotor) aufgebaut. Hier konnten die ein- 
zelnen Komponenten sowie das Gesamtsystem auf ihr Schwingungs- 
verhalten untersucht und mit Ergebnissen aus den Experimenten vergli- 
chen werden. 


Aus dem mechanischen und dem magnetischen Modell wurde ein mul- 
tiphysikalisches Gesamtmodell entwickelt, welches Aufschluss über 
die Quellen und Ursachen vibroakustischer Effekte gibt. An diesem 
Modell wurden mögliche Gegenmaßnahmen getestet und bewertet. 


Thema: Kompensation durch akustisch optimierte Ansteuerung 


Der aktuelle Stand der Technik zu Möglichkeiten einer akustisch opti- 
mierten Ansteuerung wurde erläutert. Darauf aufbauend wurde eine 
Steuerungsmethode basierend auf aufgeprägten Schwingungen in den 
Strömen entwickelt und am Gesamtmodell getestet. Hierfür wurde eine 
Berechnungsmöglichkeit zur Bestimmung der vibroakustisch optimier- 
ten Parameter (Amplitude und Phasenlage) aufgezeigt, welche im Ge- 
samtmodell verwendet werden kann. Für den Einsatz dieser Methode 
am Prüfstand ist ein Paramtersuchalgorithmus entwickelt worden, wel- 
cher einfach und effizient auch für andere Maschinen nutzbar ist. Dieser 
ermöglicht, neben der Verbesserung der vibroakustischen Eigenschaf- 
ten, Informationen über die Herkunft der störenden Effekte zu erhalten 
(Stromoberschwingungen oder magnetisches Design). Anschließend 
wurde die vibroakustisch verbesserte Ansteuerungsmethode auf ver- 
schiedenen Stromregelverfahren angewandt und untersucht. Für das 
System des Radnabenantriebs wurde eine Kompensationsstrategie ent- 
wickelt, welche aber auch für andere elektrische Maschinen angewandt 
werden kann. Ebenso wurde ein Nachweis der Wirksamkeit am Ge- 
samtfahrzeug mit geänderter Topologie der elektrischen Maschine er- 
bracht. 
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7.1 Zusammenfassung 


Somit ergibt sich aus dem durch die Arbeit erworbenen Wissen folgen- 
der Nutzen in konkreten Anwendungsbeispielen: 


e Verbesserung der vibroakustischen Eigenschaften durch Anpas- 
sung von Symmetrieeigenschaften in der Struktur und Magnetik 

e Aufdecken des Minderungspotentials der von der Stromrege- 
lung verursachten vibroakustischen Störungen und dessen An- 
wendbarkeit für aufgeprägte Schwingungen 

e Verbesserung des vibroakustischen Verhaltens in kritischen Be- 
triebspunkten des Radnabenantriebs nach dem magnetischen 
und mechanischen Design-Freeze durch akustisch optimierten 
Steuerungsansatz 

e Entwicklung einer universell einsetzbaren Methodik, um das 
vibroakustische Verhalten elektrischer Maschinen zu verbes- 
sern 
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7 Zusammenfassung und Ausblick 


7.2 


Ausblick 


Aus der Arbeit und der darin beschriebenen Modelle und entwickelten 
Methoden und Informationen ergeben sich aber auch neue Herausfor- 


derungen, die nicht vollständig gelöst oder abgearbeitet werden konn- 
ten und somit Potenzial für weitere Forschungsarbeiten liefern. 
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e Verbesserung des Simulationsmodells, um bereits vor der Er- 
stellung erster Prototypen vibroakustisch kritische Mechanis- 
men erkennen und besser einschätzen zu können und auch wei- 
tere mechanische Optimierungspotentiale aufzudecken 


oO 
O 


Materialparameter 
Einfluss von Exzentrizitäten und Fertigungstoleran- 
zen auf das Gesamtsystem 


e Erweiterung des Modells um die Schallabstrahlung bis zum 
wahrnehmbaren Schalldruck in die Fahrzeugkabine 


e Verbesserung des Regelverfahrens zum Aufprägen harmoni- 
scher Schwingungen 


O 


Weniger Beeinflussung bei mehreren Kompensatio- 
nen mit verschiedenen Ordnungen 

Keine oder weniger Nebenordnungen 
Implementierung einer Regelung von Phase und 
Amplitude 


e Erweiterung der Betriebsstrategie 


O 


Kompensationswirkung auf die Fahrzeuginsassen 
(z. B. durch Probandenstudie zur Erfassung der kriti- 
schen Punkte), um Betriebsstrategie zu verbessern 
Entwicklung eines Komfortmodus (nur notwendig, 
wenn sehr hohe Amplituden zur Kompensation nötig 
wären), um den eventuellen erhöhten Energiebedarf 
zu minimieren 


7.2 Ausblick 


Weitere Anwendungsgebiete der modellgestützten Suchmethodik 


Mit der Reduzierung beziehungsweise Kompensation von Beschleuni- 
gungen an der Längslenkerbuchse oder der Dämpfer- und Federkräfte 
kann auf ein verbessertes vibroakustisches Verhalten im Fahrzeugin- 
nenraum geschlossen werden. Da jedoch für unterschiedliche Betriebs- 
punkte unterschiedliche Anregungsordnungen und auch unterschiedli- 
che Übertragungswege in der Fahrgastzelle kritisch sind, ist eine sinn- 
volle Ergänzung zu den bereits gewonnenen Erkenntnissen, wie zum 
Beispiel der Beschleunigungen am Sitz oder Lenkrad, um auch das ge- 
samte vibroakustische Verhalten und deren Beeinflussungsmöglichkei- 
ten innerhalb der Fahrzeugkabine zu untersuchen. Aus diesen Untersu- 
chungen kann sich auch eine Konstruktionsmethode (Mechanik und 
Motortopolgie) entwickeln lassen, bei welcher aufgeprägte Schwingun- 
gen als Kompensationsmöglichkeit vorgesehen sind, um somit Zielkon- 
flikten (zum Beispiel zwischen Akustik und Wirkungsgrad) zu entgeg- 
nen. 


Weiterhin kann diese Methode genutzt werden, um an einem System 
eine Transferpfadanalyse oder Modalanalyse durchzuführen. Grund- 
sätzlich werden diese Analysen mit Hilfe von externen Anregungsquel- 
len, wie einem Shaker oder Impulshammer durchgeführt. Der Vorteil 
bei der Anregung durch die elektrische Maschine ist jedoch, dass diese 
das System in Bereichen anregen kann, die sehr wichtig, aber auch für 
externe Anregungsquellen schwer oder gar nicht erreichbar sind. Au- 
Berdem kann dadurch das System im laufenden Betrieb analysiert wer- 
den, was einen weiteren Vorteil darstellt. 


Ebenso ergibt sich die Möglichkeit, durch gezieltes Aufbringen von 
Schwingungen das Fahrzeug besonders sportlich wirken zu lassen, wie 
zum Beispiel bei stark beschleunigter Fahrt. Während des Projektes hat 
sich auch gezeigt, dass auch der direkt vom Motor abgestrahlte Schall 
beeinflusst werden kann. Somit kann auch die von außen hörbare Akus- 
tik durch das Aufprägen von Schwingungen positiv beeinflusst werden. 
Ein weiteres Beispiel für die Verwendung aufgeprägter Schwingungen 
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7 Zusammenfassung und Ausblick 


ist, akustische Signale nach außen direkt vom Motorgehäuse abzuge- 
ben. Dies ist besonders sinnvoll, um die seit 2019 in Europa bestehen- 
den gesetzlichen Bestimmungen zur Abgabe von Geräuschen zum Fuß- 
gängerschutz bei niedrigen Geschwindigkeiten zu erfüllen, was mit auf- 
geprägten Schwingungen in den Strömen auch ohne zusätzlichen 
Aktuator möglich sein kann. 


Weiter ergibt sich die Möglichkeit, auch die spürbaren Effekte, welche 
nicht aus dem elektrischen Radnabenantrieb stammen, durch aktive 
Kompensation zu minimieren. Vor allem bei periodischen Signalen, die 
direkt mit der Motordrehzahl zusammenhängen, wären weiterführende 
Untersuchungen sehr interessant. So kommen bei Elektrofahrzeugen 
mit zentralem Elektromotor häufig Getriebestufen zum Einsatz, um die 
Vorteile eines hochdrehenden Elektromotors zu nutzen. Durch die Ein- 
griffe der Verzahnung können hier Schwingungen entstehen, welchen 
mit aufgeprägten Schwingungen entgegengewirkt werden könnte. 
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A Anhang 


Flexible Ansteuereinheit 


Um geeignete und optimierte Ansteuerverfahren am Prüfstand testen zu 
können, wurde für diese Arbeit eine Ansteuereinheit aufgebaut, mit der 
auf alle Regelparameter Einfluss genommen werden kann. Für die Re- 
gelung wurde auf ein vom KIT am Elektrotechnischen Institut entwi- 
ckeltes Signalprozessorsystem zurückgegriffen. Die Regelung wird in 
Simulink entworfen und kann über den Simulink Coder auf dem DSP 
ausgeführt werden. Mit einem Bussysstem kann der DSP Daten von 
A/D Wandlerkarten lesen, um Ströme, Winkel oder andere physikali- 
sche Werte zu messen und in der Regelung zu verarbeiten. Als Leis- 
tungselektronik kommt das HybridKit 2 von Infineon zum Einsatz, wel- 
ches über den Signalprozessor gesteuert wird. 


Sollwertvorgabe über PC 


DIGITALER SIGNAL PROZESSOR ng E 

Implementierung der Regelung = Zi m. 

Berechnung des Winkelsignales Ti 

MODULATORKARTE A/D-WANDLERKARTE RESOLVERKARTE 

Generierung der PWM Strangströme/ Erzeugen der Erregerspannung 
Signale Inverter Zwischenkreisspannung Filterung/AD-Wandlung der Werte 
INVERTER 


HybridKit Infineon 


RESOLVER 


Abbildung 7.1: Aufbau der Ansteuereinheit 
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Anhang 


Prüfstandsadapter für Fahrwerk 


Ebenfalls wurde im Zuge dieser Arbeit ein Prüfstandsadapter für den 
Allrad-Rollenprüfstand des Instituts für Fahrzeugsystemtechnik am 
KIT konstruiert und aufgebaut, mit dem die Möglichkeit geschaffen 
wurde, die bei einem Fahrzeug auf die Karosserie wirkenden Wechsel- 
kräfte genauer zu untersuchen. Genutzt wird der Adapter mit dem aus 
dem Versuchsträger stammenden Fahrwerk. Der Adapter wurde sehr 
steif ausgelegt, sodass dieser wenig Einfluss auf das vibroakustische 
Verhalten des Fahrwerks und des Radnabenantriebs hat. Durch die 
Möglichkeit Sand am Rücken der Adapterplatte in die dafür vorgesehe- 
nen Rohre und damit Dämpfung einzubringen, wird der Einfluss weiter 
minimiert. Zusätzlich besteht die Möglichkeit weitere versteifende Ele- 
mente in den Adapter einzubringen, um bei der Überlagerung von Re- 


sonanzen zwischen Eigenschwingungen des Adapters und des Prüflings 
unterscheiden zu können. 


Abbildung 7.2: Aufbau des Prüfstandsadapters 
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Anhang 


Die auf die Karosserie beziehungsweise auf den Prüfstandsadapter wir- 
kenden Kräfte können über eine Kraftmessdose, die hinter dem Dämp- 
fer oder der Feder angebracht wird, gemessen werden. Weiterhin gibt 
es die Möglichkeit, über Beschleunigungssensoren an der Längslenker- 
buchse, Feder oder Dämpfer, die Wechselkräfte abzuschätzen. 


Motorenprüffeld 


Eine weitere Möglichkeit den Radnabenantrieb auf sein vibro-akusti- 
sches Verhalten zu untersuchen ist das Motorenprüffeld des Elektro- 
technischen Instituts (Abbildung 7.3). Dort besteht die Möglichkeit 
durch eine angeschlossene Abtriebsmaschine das Verhalten auch ohne 
den Einfluss des Fahrwerkes oder Reifens zu betrachten. Insbesondere 
war dort die Möglichkeit gegeben durch die verbaute Messwelle (3) 
auch das Drehmoment zu untersuchen. Allerdings ist man aufgrund der 
Rotationsmassen und Torsionssteifigkeit auf kleine Frequenzen be- 
schränkt, da bei höheren Frequenzen das Drehmomentrippel entkoppelt 
wird. Zusätzlich wurde der Aufbau mit Beschleunigungssensoren be- 
klebt, um auch diese am Statorgehäuse (1) und vor allem an der Stator- 
platten (2) genauer untersuchen zu können. 


Abbildung 7.3: Aufbau am Motorenprüffeld 
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Anhang 


Gesamtfahrzeug am Rollenprüfstand 


Zum Test der vibroakustisch optimierten Ansteuerungsmethode wurde 
der Versuchsträger, ein Ford Fiesta mit dem Radnabenantrieb der Ver- 
sion „gamma“, am akustisch verkleideten Rollenprüfstand der Schaeff- 
ler Engineering GmbH in Werdohl getestet. Um die Einflüsse der ein- 
zelnen Antriebe testen zu können, wurde hierfür eine Vorrichtung mit 
Luftbalken konstruiert und aufgebaut, sodass der Versuchsträger mit 
nur einem Radnabentrieb untersucht werden kann. Für die Messungen 


des vom Fahrer wahrgenommenen Geräusches wurden im Fahrzeugin- 
nenraum auf Kopfhöhe drei Mikrofone installiert, welche den Luft- 
schall aufnehmen. 


j 
>` 
a 


Abbildung 7.4: Abbildung des Versuchsaufbaus am Gesamtfahrzeug 
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Anhang 


Messergebnisse: Kompensationspotential 
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Abbildung 7.5: Kompensationspotential an LLBx 
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Abbildung 7.6: Kompensationspotential am Dämpfer 
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Beschleunigung [m/s?] 
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Abbildung 7.8: Kompensationspotential an Statorplatte und Statorge- 
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Messergebnisse: Trennung der Quellen 
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Abbildung 7.9: Trennung der Quellen am Dämpfer 
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